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In  searching  for  a  thesis  topic  I  found  that  there  is 
a  great  interest  within  the  Air  Force  Flight  Dynamics  Lab¬ 
oratory  at  Wright-Patterson  AFB  in  the  concept  of  self¬ 
repairing  flight  control.  This  concept  not  only  holds  prom¬ 
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graduate  of  the  Air  Force  Institute  of  Technology) .  I  would 
especially  like  to  thank  my  advisor.  Dr.  John  D'Azzo,  and  my 
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Abstract 

Two  linearized  models  containing  coupled  aircraft 
equations  are  developed  for  the  AFTI/F-16.  The  first  is  a 
model  of  the  healthy  aircraft  with  all  control  surfaces 
intact,  and  the  second  is  a  model  of  the  aircraft  with  a 
free-floating  right  horizontal  tail  and  all  other  surfaces 
operational . 

The  multivariable  design  technique  of  Professor  Brian 
Porter  and  the  computer  program  MULTI  are  first  used  to 
design  control  laws  for  the  healthy  model.  The  control  laws 
are  tailored  to  perform  seven  maneuvers  at  four  flight  con¬ 
ditions.  Maximum  maneuvers  are  commanded  to  yield  maximum 
control  surface  deflections. 

The  same  control  law  designs  are  then  applied  to  the 
model  with  a  failed  right  horizontal  tail,  and  the  perfor¬ 
mance  is  evaluated.  Some  maneuvers  require  modifications  to 
the  designs  or  lowered  maximum  maneuver  requirements  to 
avoid  overshooting  the  deflection  limits  of  the  operational 
control  surfaces. 

Simulation  responses  are  presented  for  both  the  healthy 
and  failure  aircraft  models.  Generally,  when  the  right  hor¬ 
izontal  tail  fails,  the  left  horizontal  tail  assumes  primary 
pitch  control  and  the  flaperons  take  over  complete  roll 
control.  The  flaperons,  rudder,  and  canards  deflect  to 
counter  the  rolling  and  yawing  moments  produced  by  the  left 
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horizontal  tail  deflection. 


The  design  process  used  in  the  study  is  applicable  to 
other  control  surface  failures.  Recommendations  are  pre¬ 
sented  for  future  study. 
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MULTIVARIABLE  CONTROL  LAW  DESIGN  FOR  THE  AFTI/F-16 


WITH  A  FAILED  CONTROL  SURFACE 


I.  Introduction 


1-1 .  Background 

In  recent  years  much  attention  has  been  directed 
towards  the  design  of  flight  control  laws  which  maintain 
aircraft  stability  after  a  control  surface  failure.  This 
topic  is  of  special  interest  to  the  Air  Force  because  its 
aircraft  may  be  reguired  to  fly  in  hostile  environments 
where  ground  fire,  air-to-air  missiles,  or  other  hazards  may 
damage  control  surfaces.  More  specifically,  the  Air  Force 
Flight  Dynamics  Laboratory  at  Wright-Patterson  Air  Force 
Base  has  begun  a  program  which  stresses  the  need  for  control 
laws  that  can  handle  surface  failures  in  future  aircraft 
designs . 

The  concept  behind  the  flight  Dynamics  Laboratory's 
program  is  called  "self-repairing  flight  control".  This 
concept  promises  benefits  to  the  Air  Force  by  improving 
safety  for  the  pilots  and  by  increasing  aircraft  availabil¬ 
ity  in  wartime.  With  self-repairing  control  an  aircraft 
will  be  able  to  continue  its  primary  mission,  or  take  on  a 
secondary  mission,  even  with  damaged  control  surfaces  or 
failed  actuators.  Self-repairing  control  will  greatly 
reduce  the  need  for  field  maintenance,  and  allow  deferment 


of  repair  of  minor  damage.  Self-repairing  therefore 
increases  aircraft  availability  and  increases  the  possible 
number  of  sorties  for  each  plane.  Additionally,  the  concept 
will  reduce  the  cost  of  aircraft  systems  by  reducing  the 
amount  of  hardware  redundancy  needed  for  safe  flight 
(13:141-146) . 

As  can  be  expected,  self-repairing  flight  control  is  a 
broad  concept  that  involves  several  engineering  disciplines. 
Much  work  must  be  done  in  order  to  develop  a  self-repairing 
flight  control  system  that  performs  the  following  necessary 
tasks.  First,  control  surface  failures  must  be  detected  and 
isolated  in  some  way  so  that  the  system  knows  which  surface 
has  failed.  Second,  a  way  must  be  found  to  estimate  the 
extent  of  the  damage  or  failure.  For  example,  damage  to  a 
surface  can  range  from  a  small  hole  to  complete  destruction 
of  the  surface.  Severe  damage  may  threaten  the  safety  of 
the  pilot  while  slight  damage  may  allow  continuation  of  the 
mission  but  with  degraded  aircraft  performance.  Third,  the 
flight  control  system  must  do  something  to  minimize  the 
effects  of  the  damage.  The  system  must  at  least  keep  the 
aircraft  stable  so  that  the  pilot  can  fly  back  to  friendly 
territory.  Finally,  the  system  must  somehow  inform  the 
pilot  of  what  has  happened  and  what  is  the  new  status  of  the 
system  (13 : 147) . 

Before  too  much  time,  effort,  and  money  are  spent  on 
the  first,  second,  and  last  tasks  listed  above,  it  is  impor¬ 
tant  to  ascertain  that  the  aircraft  is  capable  of  handling 
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control  surface  failures  in  the  first  place.  Advanced 
aircraft  designs,  such  as  the  AFTI/F-16,  have  Control 
Configured  Vehicle  (CCV)  capabilities  which  enhance  combat 
effectiveness  by  allowing  the  aircraft  to  perform  unconven¬ 
tional  maneuvers  such  as  fuselage  pointing  and  translation. 
The  maneuvers  are  accomplished  by  the  use  of  additional 
control  surfaces,  like  canards,  and  by  independent  operation 
of  left  and  right  elevators  and  ailerons.  In  this  way 
surfaces  which  previously  could  only  affect  longitudinal 
movements  can  now  affect  lateral  motions,  and  vice  versa. 
With  these  CCV  capabilities,  other  control  surfaces  can 
apply  the  forces  and  moments  needed  for  stability  when  one 
surface  is  missing  or  damaged. 

While  CCV  capabilties  offer  outstanding  payoffs  in 
combat  effectiveness  and  survivability,  design  of  flight 
control  laws  for  such  aircraft  becomes  more  difficult 
because  of  a  large  number  of  required  feedback  loops  and  the 
multiple  input/output  nature  of  the  problem.  Design  is 
further  complicated  when  considering  the  loss  of  a  control 
surface  because  of  cross-coupling  between  the  longitudinal 
and  lateral  aircraft  equations.  The  multivariable,  high- 
gain,  error  actuated  design  technique  of  Professor  Brian 
Porter  from  the  University  of  Salford,  England,  is  of 
interest  to  the  Air  Force  Flight  Dynamics  Laboratory  as  a 
technique  for  handling  CCV  designs.  Professor  Porter's 
design  method  is  outlined  in  Appendix  A. 


1-2.  Problem 


Loss  of  a  control  surface  may  cause  instability  and 
degraded  performance  in  modern  air  superiority  aircraft 
unless  the  flight  control  system  is  designed  to  handle  such 
a  change  in  the  aircraft's  configuration.  The  purpose  of 
this  research  is  to  present  a  technique  for  evaluating  the 
performance  of  an  aircraft  after  a  control  surface  failure 
has  occurred.  Control  laws  are  developed  for  a  healthy 
aircraft  model  using  the  Porter  technique,  and  then  new 
control  laws  are  designed,  also  using  the  Porter  technique, 
to  provide  maximum  performance  with  a  particular  surface 
failure . 

1-3 .  Scope 

The  control  surface  failure  considered  in  this  study  is 
a  free-floating  right  horizontal  tail  half.  The  failed 
surface  is  assumed  free-floating  so  that  it  does  not  create 
forces  and  moments  on  the  aircraft.  Although  this  investi¬ 
gation  is  limited  to  one  control  surface  failure,  the 
modeling  and  design  procedures  used  can  be  applied  to  other 
failed  surfaces. 

Aircraft  performance  is  evaluated  at  four  flight 
conditions:  subsonic  (0.6  Mach  at  30,000  feet),  transonic 

(0.9  Mach  at  20,000  feet),  supersonic  (1.6  Mach  at  30,000 
feet),  and  landing  approach  (0.3  Mach  at  30  feet).  These 
flight  conditions  are  chosen  because  they  sample  a  variety 
of  points  in  the  aircraft's  flight  envelope.  Seven 


maneuvers  are  considered  in  this  study  for  the  evaluation  of 
aircraft  performance.  The  maneuvers  are:  g-command,  pitch 
pointing,  longitudinal  translation,  roll,  sideforce,  yaw 
pointing,  and  lateral  translation.  These  maneuvers  are  des¬ 
cribed  in  detail  in  Chapter  IV. 

1-4 .  Assumptions 

For  most  preliminary  designs  using  linearized  pertur¬ 
bation  equations  the  longitudinal  and  lateral  equations  of 
motion  are  considered  to  be  decoupled.  However,  in  this 
study  the  lateral  and  longitudinal  equations  are  coupled 
from  the  start  because  an  unsymmetrical  control  surface 
failure  causes  cross-coupling  between  the  two  sets  of 
equations . 

The  engine  thrust  model,  discussed  in  Chapter  II,  is 
assumed  to  be  a  first  order  lag  and  is  included  only  to  hold 
the  forward  velocity  constant  during  maneuvering.  Although 
two  vertical  canards  are  present  on  the  AFTI/F-16,  the 
canards  are  considered  to  be  one  surface  in  order  to  simpli¬ 
fy  the  designs.  The  remaining  assumptions  are  similar  to 
those  used  in  a  recent  thesis  by  Mr.  Finley  Barfield  (4:6) 
and  are  listed  below. 

*  The  aircraft  is  a  rigid  body  with  constant  mass. 

*  The  atmosphere  is  fixed  with  respect  to  the  earth. 

*  The  earth's  surface  is  an  inertial  reference  frame. 

*  Linearization  about  an  operating  condition  is  valid 
for  preliminary  aircraft  models. 


*  Aerodynamics  are  constant  for  Mach  ani  altitude. 

The  constant  mass  assumption  is  valid  because  the 
system  responses  occur  over  a  relatively  short  period  of 
time.  The  rigid  body  assumption  is  valid  as  long  as  the 
resulting  design  bandwidths  sufficiently  limit  noise  due  to 
structural  interactions.  The  other  assumptions  greatly 
simplify  the  design  and  are  widely  accepted  in  preliminary 
design  efforts.  Further  investigations  can  take  into  ac¬ 
count  the  effects  of  noise  and  wind  gusts,  but  these  effects 
are  not  considered  in  this  study. 

1-5 .  Approach 

The  AFTI/F-16  is  chosen  for  this  study  because  some 
control  laws  for  the  healthy  aircraft  have  been  developed  in 
Barfield's  thesis  (4),  and  because  most  of  the  needed  infor¬ 
mation  on  the  aircraft  is  readily  available  from  the  Flight 
Dynamics  Laboritory  (1).  Two  linear  models,  one  for  the 
healthy  aircraft  and  one  for  the  aircraft  with  a  failed 
right  horizontal  tail,  are  developed  for  the  AFTI/F-16  and 
placed  in  state  space  form. 

Control  laws  are  designed,  using  the  Porter  technique 
and  the  computer  program  MULTI,  for  the  healthy  aircraft 
model  at  the  four  flight  conditions.  The  designs  are 
tailored  to  the  maneuvers  performed  at  each  flight  condi¬ 
tion.  The  designs  are  then  applied  to  the  failure  model, 
and  the  aircraft  performance  is  evaluated.  If  necessary, 
the  control  laws  are  redesigned  to  achieve  optimal 


performance  and  to  keep  surface  deflections  within  the 
appropriate  limits. 

1-6 .  Overview 

Appendix  A  summarizes  the  theory  of  the  multivariable 
design  technique  of  Professor  Porter,  and  Appendix  B  lists 
the  equations  used  in  developing  the  aircraft  models.  The 
actual  models  used  for  each  flight  condition  are  given  in 
Appendix  C  which  also  presents  data  files  used  with  the 
computer  program  MULTI. 

Chapter  II  begins  with  a  description  of  the  AFTI/F-16 
aircraft  and  developes  the  healthy  and  failed  models  used  in 
this  study.  Actuator  and  thrust  models  are  also  discussed 
in  Chapter  II.  Chapter  III  describes  in  detail  the  process 
used  in  designing  the  control  laws  for  this  study,  and 
Chapter  IV  describes  the  maneuvers  considered. 

The  simulation  responses  for  the  transonic  flight 
condition  are  presented  in  Chapter  V,  and  the  simulation 
responses  for  the  remaining  flight  conditions  are  given  in 
Appendix  D.  Chapter  VI  concludes  the  study  with  a  general 
discussion  of  results  and  recommendations  for  future  study. 


II.  Aircraft  Description  and  Models 

II-l.  AFTI/F-16  Description 

The  AFTI  (Advanced  Fighter  Technology  Integration) /F-16 
aircraft,  shown  in  Figure  2-1,  is  an  F-16A  air  superiority 
fighter  modified  to  be  a  testbed  for  evaluating  new  tech¬ 
nologies  (3:1).  The  modifications  include:  addition  of 
two  vertical  canards  mounted  on  the  engine  inlet;  modifica¬ 
tions  to  the  functions  of  existing  control  surfaces  by 
allowing  independent  motion  of  the  trailing  edge  flaps  and 
independent  motion  of  the  horizontal  tail  halves;  and 
addition  of  a  redundant,  digital  fly-by-wire  flight  control 
system.  The  AFTI/F-16  is  not  an  operational  fighter,  but 
the  technologies  tested  on  it  will  be  employed  in  the  Air 
Force's  future  Advanced  Tactical  Fighter  (ATF)  aircraft. 

By  design  the  unaugmented  aircraft  is  statically  un¬ 
stable  in  the  longitudinal  axis  for  subsonic  flight.  This 


Figure  2-1:  The  AFTI/F-16 
Source:  (1) 


is  because  the  center  of  gravity  is  located  behind  the 
aerodynamic  center  of  the  aircraft.  The  instability,  mani¬ 
fested  by  an  unstable  short  period  root,  allows  the  aircraft 
to  withstand  higher  load  factors  and  reduces  drag.  Also, 
the  Dutch  roll  mode  of  the  aircraft  is  lightly  damped.  It 
is  the  primary  purpose  of  the  flight  control  system  to  sta¬ 
bilize  the  aircraft  longitudinally  and  improve  the  Dutch 
roll  damping  (4:29-31). 

Like  other  fighters  in  the  Air  Force  inventory,  the 
AFTI/F-16  can  perform  conventional  maneuvers,  such  as  pitch¬ 
ing  longitudinally,  rolling  laterally,  and  turning  with  zero 
sideslip  (coordinated  turn) .  AFTI/F-16  pilots  can  also 
command  maneuvers  that  require  decoupling  of  the  aircraft's 
forces  and  moments  which  cannot  be  commanded  in  conventional 
fighters,  like  pitch-pointing,  yaw  pointing,  and  lateral  and 
longitudinal  translation.  Conventional  and  decoupled 
maneuvers  can  be  blended  to  give  the  AFTI/F-16  a  versatile 
inventory  of  maneuvers. 

The  control  surfaces  used  in  maneuvering  the  aircraft 
are  shown  in  Figure  2-2.  The  surfaces  that  this  thesis  is 
concerned  with  are  shown  in  heavy  black  and  include  the  left 
and  right  horizontal  tail  halves,  left  and  right  flaperons, 
the  rudder,  and  the  left  and  right  canards.  The  horizontal 
tail  halves  can  deflect  symmetrically  as  elevators  to  pitch 
the  aircraft,  or  deflect  asymmetrically  to  augment  rolling. 
Similarly,  the  flaperons  also  have  a  dual  function.  They 
can  deflect  symmetrically  as  flaps  or  asymmetrically  as 
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ailerons  for  primary  roll  control  (4:29).  For  the  purposes 
of  this  study,  the  canards  are  considered  to  be  one  surface 
that  provides  sideforces  on  the  aircraft.  Actually,  the 
canards  can  also  be  used  in  a  snowplow  configuration  as  a 
speed  brake.  The  rudder  is  used  for  yawing  as  on  conven¬ 
tional  aircraft.  Figure  2-2  defines  the  positive  direction 
of  surface  deflections. 

The  AFTI/F-16's  maneuvers  are  tailored  to  primary 
mission  tasks  which  are  divided  into  four  major  configura¬ 
tions  of  the  flight  control  system.  The  first  configuration 
is  the  normal  mode  which  is  used  for  takeoff,  cruise,  land¬ 
ing,  and  air  refueling.  The  second,  air-to-air  gunnery 
mode,  lets  the  pilot  command  precise  pointing  to  targets  and 
evasive  maneuvering.  The  air-to-surface  gunnery  mode,  the 
third  configuration,  provides  pointing  for  strafing  ground 
targets  as  well  as  evasive  maneuvering  for  greater  surviv¬ 
ability.  The  last  configuration  of  the  flight  control 
system,  air-to-surface  bombing  mode,  allows  the  pilot  to 
precisely  control  the  aircraft's  velocity  vector  for  accu¬ 
rate  bombing  (3:4).  This  study  investigates  the  design  of 
control  laws  for  the  air-to-air  gunnery  mode. 

The  many  advanced  technologies  incorporated  into  the 
AFTI/F-16  set  it  apart  from  conventional  fighters.  Lessons 
learned  through  the  use  of  the  AFTI/F-16  testbed  will  aid 
the  Air  Force  in  developing  its  next  generation  of  fighter 


aircraft. 
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1 1  —  2 .  Aircraft  Models 

Reference  4  is  the  main  source  of  information  for  the 
development  of  the  aircraft  model  for  this  study.  The 
aircraft  is  modeled  by  a  set  of  first  order  differential 
equations  in  the  state  space  form 

x  =  Ax  +  Bu  (2-1) 

The  state  equations  are  derived  from  the  aircraft  equations 
of  motion  which  consist  of  forces  at  the  center  of  gravity 
(CG)  and  moments.  The  detailed  derivations  of  the  state 
perturbation  equations  are  given  in  Reference  4. 

Equation  (2-2)  is  the  state  space  model  for  the  healthy 
aircraft,  where  the  primed  terms  are  dimensionalized  deriva¬ 
tives  in  the  body  axis.  The  equations  used  for  calculating 
the  primed  derivatives  are  listed  in  Appendix  B.  The  body 
axis  reference  frame  is  chosen  for  the  design  because  all 
accelerations  and  rates  measured  for  feedback  signals  are 
sensed  in  the  body  axis,  and  because  pilots  desire  to  com¬ 
mand  body  accelerations  and  rates  that  they  can  feel  (6:10). 
Figure  2-3  defines  the  body  axis  as  well  as  positive 
directions  of  forces,  moments,  and  angles. 

For  the  healthy  aircraft  model  of  Equation  (2-2)  the 
surfaces  are  defined  as  follows: 

5et  is  the  total  elevator  where  both  horizontal  tail 
halves  deflect  symmetrically 

6ffc  is  the  total  flap  where  both  flaperons  deflect 
symmetrically 
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Figure  2-3:  Positive  Directions  of  Forces,  Moments,  and 
Angles  in  the  Body  Reference  Frame 


Source:  (1) 


5a  is  the  combination  of  differential  flaperon  (6DF) 
and  differential  tail  (6DT)  deflections 

5r  is  the  rudder  deflection 

5c  is  the  total  symmetrical  canard  deflection 

The  aileron  deflection,  6a,  demonstrates  that  the  dif¬ 
ferential  tail  augments  the  flaperons  in  rolling  the 
aircraft.  The  equation  governing  the  aileron  deflection  is 

a  =  6  DF  +  0.25  (6DT)  (2-3) 

where  0.25  is  an  arbitrary  constant  corresponding  to  the 
reciprocal  of  a  gain  on  the  differential  horizontal  tail 
actuators.  If  the  constant  is  decreased  the  actuator  gain 
increases,  which  increases  the  contribution  of  the  differ¬ 
ential  tail  in  rolling  the  aircraft.  The  reason  for  model¬ 
ing  the  aileron  in  this  way,  and  not  splitting  the  flaperons 
and  tail,  is  due  to  a  lack  of  necessary  feedback  signals. 
This  is  discussed  in  detail  in  the  next  chapter. 

Equation  (2-4)  is  the  model  for  the  aircraft  with  a 
failed  right  horizontal  tail  half.  In  this  model  the  flap¬ 
erons  and  tail  are  split  with  the  control  surfaces  defined 
as  follows: 

6e.  is  the  left  elevator  (left  horizontal  tail  half) 
deflection 

6fr  is  the  riqht  flaperon  deflection 

6f^  is  the  left  flaperon  deflection 

5r  is  the  rudder  deflection 

6c  is  the  total  symmetrical  canard  deflection 


A  simple  model  for  engine  thrust  is  included  in  both 
the  healthy  aircraft  and  failure  models  to  prevent  the 
aircraft  velocity  from  bleeding  off  during  maneuvers.  It  is 
assumed  that  the  thrust  is  aligned  with  the  x-axis  so  that 
no  moments  are  created,  and  that  the  engine  response  is 
modeled  by  a  simple  first  order  lag,  1/ (s  +  1). 

Table  2-1  lists  thrust  data  for  each  of  the  flight 
conditions  considered  in  this  thesis.  The  dimensionalized 
derivative,  X^T,  is  calculated  by  subtracting  the  thrust 
required  to  trim  the  aircraft  from  the  maximum  thrust  avail¬ 
able  at  a  given  flight  condition,  and  then  dividing  by  the 
mass  of  the  aircraft.  For  example,  at  0.9  Mach,  20,000  feet 

X’  =  ll6-5°° - ,^932-61  ESHBds  ,  19  254  £eet/SeC<)nd2 

ST  652.733  slugs 

The  change  in  thrust,  ST,  is  measured  in  pounds  of  thrust. 

TABLE  2-1 
Thrust  Data 


Flight  Condition 

(Mach/feet) 

Maximum  Thrust 

(pounds) 

Trim  Thrust 

(pounds) 

Y  1 

XST 

2 

(feet/sec  ) 

0.3 

/  30 

19,000 

4,449.7 

22.291 

0.6 

/  30,000 

9,000 

2,110.7 

10.555 

0.9 

/  20,000 

16,500 

3,932.6 

19.254 

1.6 

/  30,000 

20,500 

15,780.8 

7.230 

Source 


II-3.  Output  Vector 

An  output  vector,  y,  must  be  added  to  the  models  of 
Equations  (2-2)  and  (2-4)  in  order  to  command  maneuvers. 

The  output  vector  used  in  Reference  4  is  also  chosen  for 
this  study  because  the  vector  elements  are  naturally  com¬ 
manded  by  the  pilot  when  performing  maneuvers,  and  because 
all  the  outputs  can  either  be  accurately  measured  or  easily 
calculated  from  other  measurements,  which  is  necessary  for 
feedback  signals.  The  three  longitudinal  outputs  commanded 
are  forward  velocity,  pitch  rate,  and  normal  acceleration  at 
the  pilot's  station,  and  the  three  lateral  outputs  are  yaw 
rate,  roll  rate,  and  lateral  acceleration  at  the  pilot's 
station.  Six  outputs  are  required  since  there  are  six 
control  element  inputs  in  the  aircraft  models.  Forward 
velocity  is  chosen  as  an  output  because  it  is  directly 
affected  by  thrust. 

The  output  vector  is  added  to  the  state  space  aircraft 
models  with  an  output  matrix,  C,  and  a  feed-forward  matrix, 
D.  The  aircraft  models  then  take  the  form 

x  =  Ax  +  Bu 

y  =  Cx  +  Du  (2-5) 

Because  forward  velocity,  pitch  rate,  roll  rate,  and  yaw 
rate  are  already  state  variables,  they  are  easily  added  to 
the  models  by  placing  ones  in  the  appropriate  positions  in 
the  C  matrix.  The  normal  and  lateral  accelerations,  how¬ 
ever,  are  linear  combinations  of  state  equations  and  require 


some  derivation.  Reference  4  is  the  source  of  information 


for  deriving  the  acceleration  equations. 

An  expression  for  the  acceleration  along  the  z  axis  at 
the  center  of  gravity  is  found  from  the  aircraft  force 
equation 

F  =  m(W  +  pV  -  qU  -  gcos0cos$)  (2-6) 

ZCG 


thus 


JCG 


Az  =  — - —  =  W  +  pV  -  qU  -  gcos0  cos<f> 
CG 


:  2-7  ] 


Equation  (2-7)  is  written  as  a  perturbation  equation 


Azcg  ■  V  +  V  *  V  *  V  *  (Z5et)5et  *  (Z5ft,5£t  (2-SI 


where 


Z  =  c 
a  m  z 


(2-9) 


Zd  "  fUi  Cz*  *  0 
a 


(2-10) 


Zq  '  2Um  Cz 


:  2-ii ) 


Z  =  c 


(2-12) 


k  - 

9 

--  - 

Z  =  Q 

6ft  m  z6ft 

(2-14) 

n 

The  two  control  surface  derivatives,  Z.„  and  Ztf  , 

5et  6ft 

are  for 

- 

the  healthy  aircraft  model.  For  the  failed  right  elevator 

i 

model  the  control  surface  derivatives  are  Z$e  ,  Z ^ 

e£  rr 

Z„  ,  where 
s£t 

,  and 

r 

L* 

Z.  =  0.5  Z. 

5et  Set 

(2-15) 

- 

Z5fr  =  Z6f^  =  °*5  Z5ft 

(2-16) 

Equation  (2-8)  gives  the  acceleration  in  the  z 

direc- 

• 

tion  at  the  center  of  gravity,  but  the  normal  acceleration 

at  the  pilot's  station  is  desired.  Normal  acceleration  is 

-  ‘ 

in  the  minus  z  direction 

s 

(• 

A  =  -  A 

nCG  ZCG 

(2-17) 

u 

and  acceleration  at  the  pilot's  station  is  determined  from 

»•  i 

•  ■ 

A  =  An  "  (1x) (pr  _  3)  -  (*  ) (rq  +  p) 

p  CG  Y 

" 

+  dz)  (P2  +  q2) 

(2-18) 

where  l  x,  l and  4z  are  distances  in  the  body  axis 

from  the 

>; 

center  of  gravity  to  the  pilot's  station.  Assuming 

small 

values  for  i  and  i  ,  and  a  small  contribution  from 

y  z 

lateral 

- 

terms.  Equation  (2-18)  reduces  to 

-  « 

-■ .  '  .  « 

. 

• 

A  =  A  +  (l  )  q 

p  CG 

(2-19) 

^  i 
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r  f 
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A*  .  4m*  1*-.*  rn-*  m-Z  .--V.,.*-  m-J - - - - t-t.  ^  - 

or,  from  Equation  (2-17) 


A 


n. 


Ux)q  - 


'CG 


(2-20) 


Substituting  Equation  (2-8)  and  the  state  equation  for 
q  into  Equation  (2-20),  the  normal  acceleration  at  the 
pilot's  station  is  given  by 


A  =  (£  M')e  +  (£  M'  -  Z  )u  +  (Z  M'  -  Z  )a 
np  x  ®  x  u  u  *x  o  a 


+  (Z  M*  -  Z  )q  +  (Z  Mi  -  Z ,  )6e,_ 

x  q  q  ^  x  fiet  5et  u  t 


*  (*xM«  -  ZS£t)5ft  *  ,ZxMic»{C 


(2-21) 


where  the  primed  terms  are  defined  in  Appendix  B,  and  z  is 
equal  to  13.95  feet  for  the  AFTI/F-16.  Similarly,  for  the 
failed  right  elevator  model 


n 


•  'W*  *  (tx%  -  vu  *  <*xm;  -  z„>» 


+  (£  M’  -  Z  )q  +  (Z  Ml  -  Z.  >56,, 
x  q  q  M  x  {ej  5 e^  °  l 


*  UxMSfr  -  ZS£rIS£r  *  (ixMSft~  Z«f|lIS£l 


*  <»xM'c)5c 


(2-22) 


The  equation  for  the  lateral  acceleration  at  the 
pilot's  station  is  developed  in  a  similar  manner  where 


V  - 

c 


1 


\  -  A  -  (t  )  (pq  -  r)  -  (I  )  (p2  ♦  r2) 

Yp  yCG  Y 


+  (iz)  (qr  -  p)  (2-23) 


can  be  reduced  to 


Q  (• 


A  =  A  +  (£  )  r 

yp  yCG  X 


(2-24) 


For  the  healthy  aircraft  model  the  lateral  acceleration  at 
the  pilot's  station  is  given  by 

A  =  (i  Ni  +  YJ3  +  (£  N'  +  Y)p  +  (JIN'  +  Y)r 
y  x  $  8  x  p  p  ^  x  r  r 


+  <*xN6a  +  Y6a)5a  +  (*xN6r  +  Y6r)5r 


+  <*xN6c  +  Y6c)5c 


(2-25) 


and  for  the  failed  right  elevator 


(Z  N '  +  Y  ) 0  +  (l  N 1  +  Y  ) p  +  (4  N'  +  Y  ) r 
x  g  g  x  p  p  c  x  r  r 


*  'Vie,  *  Y5et»5et  *  <V«fr  *  Y«r»S£r 

*  (‘xHif.  *  Y5f,)5f£  *  “xNir  *  YSrISr 

l  l 

+  (AN'  +  Y.  )6  c 
x  6c  6c 


(2-26) 


where  the  primed  terms  are  defined  in  Appendix  B,  and  the  Y 
derivatives  are  given  by 


Yfl  =  c 
8  m  yn 
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(2-27) 
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6  r 


6  c 


as  c 

m  y6: 


as  c 
m  y5 1 


Y.  =  -  0.5  (0.25)  Y 
6e^  5  DT 


Y6fr  “  °‘5  Y6DF 


Y6f£  0,5  Y5DF 


(2-28) 

(2-29) 

(2-30) 

(2-31) 

(2-32) 

(2-33) 

(2-34) 

(2-35) 


Equations  (2-21)  and  (2-25)  are  used  to  form  the  state 
space  model  of  the  output  equation  for  the  healthy  aircraft 
given  by  Equation  (2-36) .  The  output  equation  for  the  fail¬ 
ed  right  elevator  model  is  given  by  Equation  (2-37)  which 
contains  Equations  (2-22)  and  (2-26). 

The  output  state  equations  include  a  feed-forward 
matrix,  D,  but  the  design  technique  does  not  allow  for  the 
feed-forward  matrix.  The  D  matrix  is  eliminated  by  includ¬ 
ing  actuator  states  in  the  state  space  models.  Reference  2 
gives  a  fourth-order  actuator  model  for  all  of  the  control 
surfaces  (2:  A-19 ,  A-20) : 

23 


_ (20.2)  (144.8)  (71.4)  _ 

(s  +  20.2)  ( s  +  144.8)  [ s2  +  2 ( 0 . 736 )  ( 71 . 4 ) s  +  71. 42]  (2-38) 

For  this  study.  Equation  (2-38)  is  approximated  by  a  second- 
order  model: 


(20) (72)  (2-39) 

(s  +  20)  (s  +  72) 

As  shown  in  Figure  2-4,  the  frequency  response  of  the 
second-order  model  agrees  closely  with  that  of  the  fourth- 
order  model  for  frequencies  less  than  50  radians/second. 

In  order  to  add  just  one  state  to  the  model  for  each 
actuator,  the  simplified  actuator  model  is  split  up  with  the 
20/ (s  +  20)  portion  included  in  the  state  space  aircraft 
model  and  the  72/ (s  +  72)  portion  added  to  the  model  exter¬ 
nally  with  the  program  MULTI.  The  new  actuator  inputs  are 
marked  with  a  subscript,  i,  for  "intermediate."  Equation 
(2-40)  forms  the  new  state  space  model  for  the  healthy  air¬ 
craft.  The  new  actuator  inputs  and  an  intermediate  thrust 
input  have  been  added  to  the  original  models  in  order  to 
eliminate  the  feed-forward  matrix.  Equation  (2-41)  forms 
the  new  model  for  the  aircraft  with  a  failed  right  elevator. 

The  next  chapter  discusses  in  detail  the  design  proce¬ 


dure  using  the  two  aircraft  models  developed  here.  Problems 
encountered  in  the  design  process  are  also  discussed. 
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III. 


Design  Process 


The  aircraft  models  used  in  this  study  are  developed  in 
Chapter  II.  Equation  (2-40)  is  the  model  for  the  healthy 
aircraft,  and  Equation  (2-41)  is  the  model  for  the  aircraft 
with  a  failed  right  elevator.  This  chapter  describes  how 
control  laws  are  designed  for  the  two  models  using  the  com¬ 
puter  program,  MULTI.  Appendix  A  outlines  the  theory  behind 
the  design  process,  while  this  chapter  describes  in  detail 
the  practical  aspects  of  the  design. 

III-l.  MULTI 

Because  the  multivariable  design  method  of  Professor 
Porter  involves  the  manipulation  of  matrices  and  vectors,  a 
computer  program  can  be  a  great  aid  to  the  designer.  MULTI 
is  an  interactive  program  made  up  of  approximately  70 
numbered  options  that  allow  the  user  to  perform  a  digital, 
multivariable  design  (8:2).  The  program  is  currently 
running  on  a  CDC  CYBER  computer  with  the  NOS  operating 
system.  Writing  of  MULTI  was  begun  by  Douglas  Porter  in  1981 
as  part  of  his  thesis  research  at  the  Air  Force  Institute  of 
Technology  (10) .  Several  other  AFIT  students  have  made 
modifications  and  additions  to  the  program,  and  improvements 
are  still  being  made. 

The  latest  additions  to  MULTI  were  written  by  Major 
Terry  Courtheyn  as  part  of  his  AFIT  thesis  requirements  (5). 
His  additions  include  the  formation  of  an  open-loop  transfer 
function  matrix  and  a  conversion  routine  to  change 


simulation  responses  from  radians  to  degress.  The  unit 
conversion  routine  is  used  extensively  for  the  designs  in 
this  thesis. 

Reference  8  is  the  MULTI  User's  Manual  which  describes 
the  use  of  MULTI  and  all  of  its  options;  however,  the 
designer  should  be  aware  of  subtle  details  that  are  not 
discussed  in  the  manual.  These  details  are  discussed  later 
in  this  chapter  when  needed. 

III-2.  Plant  Zeros 

For  this  study,  the  plant  consists  of  either  one  of  the 
aircraft  models  (Equations  (2-40)  and  (2-41)).  Before  the 
plant  is  entered  into  MULTI  the  system  zeros  are  checked  for 
the  presence  of  decoupling  zeros  and  the  locations  of 
transmission  zeros.  The  check  is  made  by  using  another 
interactive  computer  program  called  ZERO  (7:85-136).  if  no 
decoupling  zeros  are  present,  the  system  is  completely  con¬ 
trollable  and  observable  as  required  by  the  design  method. 

As  discussed  in  Appendix  A,  the  transmission  zeros  are  not 
altered  by  output  feedback  and  must  be  located  in  a  stable 
region  of  the  complex  plane. 

For  the  plants  of  this  study  the  program  ZERO  indicates 
that  there  are  no  decoupling  zeros  present?  therefore,  the 
plants  are  completely  controllable  and  observable.  However, 
the  system  has  four  transmission  zeros  located  at  the 
origin.  Because  the  presence  of  four  transmission  zeros  at 
the  origin  implies  the  system  is  unstable  regardless  of 


output  feedback,  more  discussion  is  required  on  this  point. 

As  pointed  out  by  Barfield  (4:83-85),  the  four  trans¬ 
mission  zeros  are  located  at  the  origin  because  of  assump¬ 
tions  made  in  forming  the  aircraft  models.  First  it  is 
assumed  that 

q  =  0  (3-1) 

and 

p  «  $  (3-2) 

The  two  integrations  implied  by  these  equations  account  for 

two  of  the  transmission  zeros  at  the  origin.  If  q  and  p 

were  commanded  with  step  functions,  8  and  $  would  ramp  to 

infinity,  thus  some  of  the  system  responses  are  unbounded 

for  bounded  inputs.  However,  for  any  practical  maneuver,  q 

and  p  are  commanded  with  functions  that  resemble  pulses 

more  than  steps  so  that  0  and  $  reach  steady-state  values. 

The  remaining  two  transmission  zeros  are  accounted  for 

by  two  other  assumptions  made  in  the  model  development.  The 

acceleration  along  the  z  axis  at  the  center  of  gravity  is 

given  by  the  same  equation  (Equation  (2-8) )  as  the  a  state 

equation.  The  normal  acceleration  at  the  pilot's  station  is 

therefore  equal  to  a  linear  combination  of  a  and  q.  As  with 

the  q  and  p  inputs  above,  a  bounded  command  of  AR  can  yield 

P 

unbounded  responses  for  the  a  and  q  states.  Similarly,  the 

lateral  acceleration  at  the  pilot's  station  is  a  linear  com- 

•  • 

bination  of  8  and  r,  so  if  A  is  commanded  by  a  step,  the 


8  and  r  responses  ramp  to  infinity.  In  practice,  though, 
acceleration  commands  resemble  pulses  which  yield  steady- 
state  values  for  the  other  system  responses. 

I I 1-3.  Entering  Actuator  Models  into  MULTI 

As  discussed  in  Chapter  II,  a  second-order  model  is 
used  for  the  control  surface  actuators.  Part  of  the  second- 
order  model  is  included  in  Equations  (2-40)  and  (2-41)  in 
order  to  avoid  having  a  feed-forward  matrix  in  the  aircraft 
models.  The  remaining  portion  of  the  actuator  model  is 
implemented  using  MULTI  Option  #4  (8:6)  and  is  external  to 
the  plant. 

MULTI  Option  #4  requires  an  external  model  for  every 
control  input.  Because  the  first-order  thrust  model  is 
already  included  in  the  plant,  a  transfer  function  of 
1000/  (s  +  1000)  is  used  for  the  external  thrust  model  for 
Option  #4.  This  external  thrust  model  approximates  unity. 

The  controllable  and  observable  plant  information  and 
external  actuator  models  are  stored  in  data  files  for  use  by 
MULTI  Option  #9  (8:7).  The  data  files  for  this  study  are 
listed  in  Appendix  C. 

III-4.  Measurement  Matrix 

The  plants  used  in  this  thesis  are  irregular  as  de¬ 
scribed  in  Appendix  A.  It  is  therefore  necessary  to  feed 
back  the  derivatives  of  certain  states  through  a  measurement 
matrix,  M.  The  measurement  matrix  should  be  chosen  such 
that  its  nonzero  elements  are  as  sparce  as  possible  and  so 


that  FB,  when  formed,  has  full  rank.  The  measurement  matrix 
is  also  required  to  yield  stable  transmission  zeros  for  the 
system. 

As  mentioned  in  Chapter  II  it  would  have  been  desirable 
to  model  the  healthy  aircraft  with  left  and  right  differen¬ 
tial  tails  and  left  and  right  flaperons  rather  than  the 
combinations  used  in  Equation  (2-40) .  It  would  have  then 
been  possible  to  model  a  failed  surface  by  eliminating  that 
surface  from  the  state  equations.  Modeling  in  this  way 
would  greatly  simplify  the  interperation  of  control  surface 
responses  when  comparing  healthy  and  failed  aircraft 
simulations . 

The  problem  in  trying  to  model  the  healthy  aircraft 
with  split  tails  and  flaperons  is  that  there  are  not  enough 
independent  feedback  sources  in  the  plant.  The  healthy 
model  would  have  required  seven  inputs  and  seven  outputs, 
and  CB  would  have  had  a  rank  deficiency  of  five.  The 
measurement  matrix  would  then  be  required  to  feed  back  five 
independent  derivatives  of  states.  Because  the  lateral 
acceleration,  An  ,  is  a  linear  combination  of  the  a  and  q 

p. 

equations,  a  and  q  cannot  be  fed  back  together  through  M. 

Similarly,  because  A  is  included  in  the  output  vector,  3 

yp 

and  r  cannot  be  fed  back  together.  Also,  because  8  and  $ 

•  • 

are  dependent  on  q  and  p  respectively,  8  and  $  cannot  be  fed 
back  through  M.  With  these  limitations,  only  four  independ¬ 
ent  derivatives  of  states  are  available  to  be  fed  back 
through  M  instead  of  the  required  five. 


Modeling  the  healthy  aircraft  as  in  Equation  (2-40) 
reduces  the  number  of  inputs  and  outputs  by  one,  and  the 
rank  deficiency  of  CB  is  four.  Keeping  in  mind  the  limita¬ 
tions  described  above,  the  four  derivatives  fed  back  through 

•  •  •  • 

M  are  u,  q,  p,  and  r.  These  derivatives  are  easily  measured 
or  calculated  from  aircraft  sensor  output  signals.  The  same 
derivatives  are  also  fed  back  for  the  failed  right  elevator 
model.  The  measurement  matrix  used  with  both  the  healthy 
and  failed  aircraft  models  is 


The  choices  of  the  nonzero  values  are  somewhat  arbitrary, 
but  the  values  are  less  than  one  in  order  to  avoid  large 
gains  in  Kq  and  which  may  cause  instabilities.  Also,  the 
locations  of  the  transmission  zeros  due  to  the  measurement 
matrix  are  approximately  equal  to  minus  the  reciprocal  of 
the  matrix  elements.  From  Equation  (3-3) ,  two  transmission 
zeros  are  added  near  -10,  and  two  more  are  added  near  -4  on 
the  real  axis. 

III-5.  Design  Parameters 

Appendix  A  describes  the  formation  of  the  gain 
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matrices,  KQ  and  K^,  for  an  irregular  plant.  After  the 
measurement  matrix,  M,  and  the  design  parameters,  <* ,  z,  and 
e,  have  been  entered  into  MULTI,  the  program  calculates  the 
gain  matrices.  There  are  no  rules  that  tell  how  to  choose 
any  of  the  design  parameters.  The  design  of  control  laws 
reduces  to  a  trial  and  error  process;  however,  the  design 
steps  discussed  below  help  speed  up  the  process. 

Because  the  longitudinal  and  lateral  equations  of  the 
healthy  aircraft  model  (Equation  (2-40) )  are  coupled  only  by 
the  M^c  term,  which  always  has  a  small  value,  the  equations 
are  considered  to  be  decoupled  for  purposes  of  choosing 
design  parameters.  The  designs  of  Reference  4  were  achieved 
for  separate  longitudinal  and  lateral  aircraft  models. 

Since  the  healthy  aircraft  model  in  this  thesis  is  essen¬ 
tially  decoupled,  the  design  parameters  used  in  Reference  4 
are  used  as  starting  points  for  the  designs  of  this  thesis. 
When  designing  a  control  law  for  a  longitudinal  maneuver, 
the  o  values  (elements  of  the  diagonal  weighting  matrix,  Z_) 
corresponding  to  lateral  quantities  are  held  constant,  while 
a  values  corresponding  to  longitudinal  quantities  are  varied 
with  successive  trials.  Similarly,  when  designing  for  lat¬ 
eral  maneuvers,  longitudinal  a's  are  held  constant  while 
lateral  o'  s  are  varied. 

For  preliminary  design  trails  the  ratio  of  proportional 
to  integral  control,  a,  and  the  weighting  matrix  multiplier, 
e,  are  set  to  one  while  the  diagonal  weighting  matrix 
elements  are  varied.  The  o  values  are  adjusted  so  that  the 
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system  is  stable  and  the  corresponding  elements  of  the 
output  vector,  y,  closely  follow  the  input  commands. 

If  satisfactory  responses  cannot  be  achieved  by  chang¬ 
ing  the  elements  of  alone,  a  is  varied  while  holding  the 
other  design  parameters  constant.  Increasing  the  value  of  a 
tends  to  shorten  the  settling  time  of  the  system  responses. 
Originally,  the  a  parameter  in  MULTI  was  implimented  as 
the  reciprocal  of  the  definition  in  Appendix  A,  which  meant 
that  the  proportional  gains,  Kg,  were  changed  while  the 
integral  gains,  ,  were  held  constant.  It  is  desired, 
however,  to  hold  Kg  constant  while  varying  ;  therefore, 
the  MULTI  program  has  been  modified  to  do  this. 

During  the  design  process  the  weighting  matrix  multi¬ 
plier,  e,  is  useful  for  finding  the  maximum  gain  values  that 
still  permit  stable  operation.  It  is  normally  set  back  to 
one,  with  the  a's  adjusted  accordingly,  for  the  final 
designs . 

III-6.  Design  Simulation 

Each  design  iteration  is  checked  with  the  simulation 
routine  in  MULTI.  Before  a  simulation  can  be  run  an  input 
command  vector,  v,  must  be  supplied  to  MULTI.  The  elements 
of  the  input  command  vector  are  the  same  as  those  of  the 
output  vector,  therefore, 

v  =  (  u,  An  ,  q,  A  ,  p,  r }T  (3-2) 

P  YP 

The  command  inputs  can  be  either  ramp  or  step  commands  but 


cannot  be  left  floating.  Because  each  element  of  v  must  be 
commanded,  some  calculations  are  necessary  when  performing 
certain  maneuvers.  The  maneuvers  considered  in  this  study 
are  described  in  the  next  chapter  along  with  the  development 
of  the  corresponding  command  vectors. 

It  is  extremely  important  that  any  control  law  design 
results  in  stable  system  responses.  It  would  have  been 
faster  and  less  costly  in  computer  time  to  check  the  closed- 
loop  transfer  functions  of  the  system  responses  for  stabil¬ 
ity  before  running  a  simulation;  however,  the  factoring 
routine  in  MULTI  Option  #6  (8:6)  becomes  error  prone  when 
transfer  functions  are  greater  than  fifth  order.  Large 
systems,  like  the  ones  In  this  study,  may  in  fact  be  stable 
even  though  Option  #6  indicates  poles  in  the  far  right  half 
plane.  It  is  suggested  that  the  problem  be  corrected  in  a 
future  research  effort. 

After  running  a  simulation  on  MULTI,  stability  is  first 
checked  by  examining  the  figures  of  merit  of  the  output  re¬ 
sponses.  Large  peak  values  in  the  figures  of  merit  usually 
indicate  an  unstable  system.  Attempting  to  plot  an  unstable 
response  with  a  MULTI  plotting  routine  may  lead  to  a  fatal 
error  and  loss  of  data  due  to  an  arithmetic  overflow.  If 
the  peak  values  are  close  to  the  commanded  steady  state 
values,  the  system  may  be  stable.  The  time  responses  of  the 
system’s  states  and  outputs  are  then  plotted  for  the  time 
interval  desired  and  checked  for  diverging  ramps  or  oscilla¬ 
tions.  It  may  be  necessary  to  increase  the  time  interval 


and  rerun  the  simulation  to  ensure  that  instabilities  do  not 
become  evident  after  a  longer  period  of  time. 

III-7.  Final  Designs 

The  trial  and  error  design  process  ends  when  desirable 
system  responses  are  achieved.  Desirable  responses  are 
those  the  that  reach  a  steady  state  value  quickly  without  an 
excessive  overshoot  and  are  free  of  oscillations.  Control 
surface  responses  must  remain  within  the  rate  and  maximum 
deflection  limits  shown  in  Table  3-1.  To  avoid  exceeding 
the  rate  limits  the  input  commands  are  ramped  to  steady 
state  values  rather  than  stepped.  The  duration  of  the  ramp 
is  determined  by  examining  the  control  surface  responses. 

The  ramp  time  is  increased  until  all  the  surface  responses 
are  within  their  respective  rate  limits. 

When  an  acceptable  control  law  is  finally  achieved  for 
a  maneuver  using  the  healthy  aircraft  model,  the  same  design 
parameters  are  tried  for  the  same  maneuver  and  flight  condi¬ 
tion  using  the  failed  right  elevator  model  of  Equation 
(2-41).  New  gain  matrices  are  calculated  because  the 
control  surfaces  are  represented  differently  in  the  failure 
model  than  in  the  healthy  aircraft  model.  Changing  the 
gains  implies  that  the  aircraft  must  be  able  to  detect  a 
failure  and  then  reconfigure  the  flight  control  system 
accordingly.  If  the  new  gains  result  in  unacceptable  re¬ 
sponses  for  the  failed  system,  the  design  parameters  are 
varied  until  acceptable  responses  are  finally  achieved.  If 


Table  3-1 

Surface  Position  and  Rate  Limits 


Surface 

Position  Limit 
(degrees) 

No  Load  Rate  Limit 
(degrees /second) 

Elevators 

±  25 

60 

Flaperons 

±  20 

52 

Rudder 

±  30 

120 

Canards 

±  27 

100 

Source:  (4:93) 


control  surface  limits  are  exceeded,  it  may  be  necessary  to 
redefine  the  maximum  maneuvers  that  the  aircraft  can  perform 
with  the  failed  surface. 

The  next  chapter  describes  in  detail  the  maneuvers 
considered  in  this  study. 


IV.  Maneuvers 


Chapter  III  discusses  the  design  method  used  in 
developing  control  laws  for  the  aircraft.  Each  control  law 
is  tailored  for  a  specific  maneuver  at  a  specific  flight 
condition.  This  chapter  describes  the  maneuvers  considered 
in  this  study  and  discusses  the  development  of  a  command 
vector,  v,  for  each  maneuver. 


IV— 1 .  Longitudinal  Maneuvers 

The  three  longitudinal  maneuvers  considered  in  this 
thesis  are:  g-command  (constant-g  pullup) ,  pitch  pointing, 
and  longitudinal  translation. 

The  g-command,  or  constant-g  pullup,  maneuver  is  used 
to  climb  from  wings  level  flight.  The  pilot  commands  a 
constant  normal  acceleration  that  he  feels  with  his  body, 
and  the  aircraft  pitches  and  climbs  while  generating  the 

constant  acceleration.  The  flight  path  angle,  y,  changes  as  *1^*- 

■*--  •  * 

the  pitch  angle  and  angle  of  attack  change,  and  is  governed 

by  the  equation  -v  v' 

Y  =  0  -  a  (4-1)  ..... 

The  command  vector,  v,  for  the  g-command  maneuver  in-  ; 

eludes  a  step  command  on  A  in  order  to  achieve  a  constant 

n 

P  -  .  ■ 

normal  acceleration  at  the  pilot's  station.  Because  none  of  .  v/.' 

the  elements  of  v  can  be  left  floating,  the  pitch  rate,  q, 

must  also  be  commanded  with  a  step.  The  magnitude  of  the 

step  is  developed  as  follows  (4:107).  For  trimmed,  wings 
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level  flight  Equation  (2-7)  can  be  written  as 


A  =  au  -  qU  (4-2) 

ZCG 

Substituting  Equation  (4-2)  into  Equations  (2-17)  and  (2-19) 


•  • 

An  =  qU  -  OU  +  (*  )q 

P 


(4-3) 


For  constant  acceleration  and  pitch  rate  commands,  the  value 
of  a  reaches  a  constant  and  thus  a  and  q  equal  zero; 
therefore 


An  (g)  =  [U  (ft/sec)] 
P 


q  (radian/sec)  (4-4) 


32.2  ft/sec 

The  magnitude  of  the  step  command  on  q  is  then  given  by 


32.2  _  .  .  c. 

q  "  U  An  (4-5) 

P 

A  step  command  on  An  and  q  implies  mathematically  that 

P 

the  commands  are  applied  for  infinite  time.  For  a  practical 
maneuver,  the  commands  are  cut  off  after  a  finite  time; 
therefore,  the  simulation  responses  are  observed  for  only  a 
few  seconds. 

As  with  the  rest  of  the  maneuvers  in  this  study,  the 
change  in  forward  velocity,  u,  is  commanded  to  zero.  Also, 
when  performing  a  longitudinal  maneuver,  the  lateral 
elements  of  v  are  commanded  to  zero. 

Although  not  attempted  in  this  study,  the  g-command 
maneuver  can  be  combined  with  the  roll  maneuver,  discussed 
later,  to  perform  a  coordinated  turn.  The  pilot  still 
commands  a  constant  normal  acceleration,  but  the  direction 


of  the  acceleration  is  at  an  angle  equal  to  the  roll  angle 
with  respect  to  a  flat  earth  reference  frame. 

The  pitch  pointing  maneuver  is  used  to  precisely  aim 
the  fuselage  of  the  aircraft  longitudinally  at  a  target. 

The  maneuver  gives  the  pilot  a  great  advantage  in  air-to-air 
combat.  When  performing  a  pitch  pointing  maneuver,  the  air¬ 
craft  pitches  without  changing  its  flight  path  angle.  From 
Equation  (4-1)  the  angle  of  attack  must  be  equal  to  the 
pitch  angle.  Also,  for  a  pitch  pointing  maneuver,  no  normal 
acceleration  is  developed. 

It  is  desirable  to  command  pitch  angle  to  demonstrate  a 
pitch  pointing  maneuver;  however,  pitch  rate,  q,  is  an 
element  of  the  command  vector  and  8  is  not.  The  pitch  rate 
is  therefore  commanded  with  a  pulse.  The  integral,  or  area, 
of  the  pulse  determines  the  pitch  angle  commanded.  All 
other  elements  of  v  are  commanded  to  zero. 

Like  the  g-command  maneuver,  pitch  pointing  can  also  be 
combined  with  a  roll  maneuver.  This  combination  is  valuable 
for  air-to-air  combat,  but  is  not  attempted  in  this  thesis. 
It  is  suggested,  though,  that  the  combination  of  maneuvers 
be  investigated  in  future  research. 

The  longitudinal  translation  maneuver  involves  the 
development  of  a  constant  normal  velocity  without  changing 
the  aircraft's  pitch  angle.  The  flight  path  angle  does 
change  and,  for  wings  level  flight,  is  equal  to  the  change 
in  angle  of  attack.  The  longitudinal  translation  maneuver 
is  useful  in  air-to-air  combat  and  air  refueling. 


In  order  to  command  a  constant  normal  velocity,  the 

normal  acceleration,  An  ,  is  commanded  with  a  pulse.  The 

P 

integral  of  the  pulse  determines  the  steady  state  value  of 
the  normal  velocity.  The  other  elements  of  v  are  commanded 
to  zero. 

IV-2 .  Lateral  Maneuvers 

Four  lateral  maneuvers  are  considered  in  this  study. 
They  are:  roll,  sideforce  (flat  turn),  yaw  pointing,  and 
lateral  translation. 

The  roll  maneuver  involves  developing  a  roll  rate  about 
the  aircraft's  velocity  vector  without  accelerating  lateral¬ 
ly  or  normally  and  without  changing  the  fight  path  angle. 

The  reason  for  rolling  about  the  velocity  vector  and  not  the 
aircraft's  x  axis  is  due  to  the  modeling  of  the  aircraft  in 
the  body  reference  frame  (4:36).  At  high  angles  of  attack, 
rolling  about  the  x  axis  results  in  large  sideslip  angles. 

In  order  to  properly  command  a  roll  about  the  velocity 
vector,  two  problems  arise.  The  first  problem  involves 
commanding  sideslip,  6  ,  to  zero,  when  it  is  not  an  element 
of  the  command  vector,  and,  at  the  same  time,  commanding  yaw 
rate,  r,  to  the  proper  value.  As  developed  in  Reference  4, 
an  estimate  for  sideslip  rate  is  derived  from  the  y  axis 
force  equation 

F  =  m(V  +  Ur  -  pW)  -  gm  cos0sin$  (4-6) 

yCG 

Rearranging  Equation  (4-6)  and  dividing  by  mass,  m. 


U(B  +  r  -  po  -  3.  cos8sin<t>) 


Solving  for  B  and  inserting  Equation  (2-24) 


■  y  (  )  r 

3  =  — ^  ^ —  -  r  +  po  +  3.  cos9  sin4> 


(4-8) 


For  constant  roll  rates,  r  is  zero.  The  equation  that 
approximates  sideslip  rate  is  finally  written  by  neglecting 
the  gravity  terms  in  Equation  (4-8) 


-  r  +  par 


(4-9) 


The  lateral  acceleration  in  Equation  (4-9)  is  commanded  to 
zero;  therefore,  for  the  roll  maneuver  Equation  (4-9) 
becomes 


B  =  paT  -  r 


(4-io : 


To  overcome  the  problem  of  commanding  B  to  zero,  the  command 
and  output  vectors  for  the  aircraft  models  are  changed  to 


X  =  X  -  {u,  An  ,  q,  Av  ,  p,  (r  -  paT)} 


!  4-11 ) 


The  output  matrix,  C,  in  Equation  (2-40)  then  becomes  the 
matrix  in  Equation  (4-12) .  A  similar  change  is  made  for  the 
output  matrix  in  Equation  (2-41).  By  then  commanding 


(r  -  paT)  to  zero,  f} ,  and  therefore  B,  is  commanded  to  zero 
also.  The  roll  rate,  p,  is  commanded  with  a  pulse  so  that 
the  steady-state  roll  angle  is  equal  to  the  integral  of  the 
pulse . 


The  second  problem  stems  from  neglecting  the  gravity 
terms  in  Equation  (4-8) .  If  the  steady-state  roll  angle  is 
less  than  ten  degrees,  the  approximation  is  valid,  but  for 
large  roll  angles  a  change  in  the  model  is  required.  As  the 
aircraft  rolls,  the  direction  of  gravity  changes  with 
respect  to  the  aircraft  body  axis.  In  order  to  maintain  a 
constant  altitude,  the  normal  and  lateral  accelerations  must 
combine  to  cancel  gravity.  Unfortunately,  to  command  the 
accelerations  properly  requires  sinusoidal  inputs  to  the 
command  vector.  At  this  time,  MULTI  does  not  allow  for 
sinusoidal  commands.  To  overcome  this  problem,  the  term 
in  the  aircraft  models  (Equations  (2-40)  and  (2-41) )  is 
replaced  by  zero.  By  doing  this,  the  lateral  equations  in 
the  models  do  not  include  a  roll  angle  term;  therefore,  the 


orientation  of  the  gravity  vector  is  not  taken  into  account 
when  rolling  the  aircraft.  With  this  change,  acceptable 
roll  maneuver  responses  are  achieved. 

The  sideforce  (flat  turn)  maneuver  is  used  for  changing 
the  heading  angle  of  the  aircraft  without  developing  a  roll 
angle  or  sideslip.  To  command  the  maneuver,  the  lateral  ac¬ 
celeration  is  commanded  with  a  step  function.  The  yaw  rate, 
r,  is  also  commanded  with  a  step  function,  but  a  with  magni¬ 
tude  determined  by  the  acceleration  command  and  aircraft 
velocity.  The  equation  governing  the  r  command  magnitude  is 
developed  from  Equation  (4-9) .  Assuming  that  no  sideslip  or 
roll  angle  is  generated  with  the  maneuver,  3  and  p  are  set 
to  zero  and  Equation  (4-9)  becomes 

r  =  nr1  Ay  (4-13) 

p 

where  the  proper  unit  conversion  is  included  as  in  Equation 
(4-4)  .  The  command  vector  for  the  sideforce  maneuver  con¬ 
sists  of  step  commands  on  lateral  acceleration  and  yaw  rate, 
and  the  remaining  vector  elements  are  commanded  to  zero. 

The  yaw  pointing  maneuver  is  similar  to  pitch  pointing, 
but  in  the  lateral  direction.  Yaw  pointing  is  used  to  aim 
the  fuselage  laterally  at  a  target  without  rolling  or  accel¬ 
erating  laterally.  Equation  (4-9)  is  used  in  forming  the 
command  vector  for  this  maneuver.  By  commanding  the  lateral 
acceleration  and  roll  rate  to  zero.  Equation  (4-9)  reduces 
to 


/  -•  -■ 
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A  desired  sideslip  angle,  or  pointing  angle  for  this 
maneuver,  can  be  developed  by  commanding  r  with  a  pulse  of 
the  opposite  sign.  The  command  vector  for  the  yaw  pointing 


maneuver  therefore  consists  of  a  pulse  on  r  and  zero  com¬ 
mands  on  the  other  vector  elements.  The  area  of  the  pulse 
on  r  determines  the  pointing  angle  commanded. 

Lateral  translation  is  similar  to  longitudinal  trans¬ 
lation,  but  in  the  lateral  direction.  When  commanding  this 
maneuver,  it  is  desired  to  develop  a  constant  lateral 
velocity  without  rolling  or  yawing.  To  perform  the  maneu¬ 
ver,  the  lateral  acceleration  is  commanded  with  a  pulse. 

The  integral  of  the  pulse  approximates  the  lateral  velocity 
developed.  The  other  elements  of  the  command  vector  are  set 
to  zero. 

Table  4-1  summarizes  the  command  vectors  for  the 
maneuvers  considered  in  this  study. 

The  next  chapter  discusses  the  simulation  responses  of 
the  maneuvers  described  in  this  chapter  at  the  four  flight 


conditions . 


Table  4-1 
Maneuver  Commands 


Maneuvers 

COMMAND 

Longitudinal 
u  An 

P 

VECTOR  ELEMENTS 

Lateral 
q  Ay 

P 

p 

r 

G-Command 
(Constant-G 
Pul lup) 

0 

step 

32.2  , 
U  A 

n  0 

P 

0 

0 

Pitch 

Pointing 

0 

0 

pulse 

0 

0 

0 

Longitudinal 

Translation 

0 

pulse 

0 

0 

0 

0 

Roll  About 

Velocity 

Vector 

0 

0 

0 

0 

pulse 

* 

0 

Sideforce 
(Flat  Turn) 

0 

0 

0 

step 

0 

32.2  , 

0  Ay 

P 

Yaw  Pointing 

0 

0 

0 

0 

0 

pulse 

Lateral 

Translation 

0 

0 

0 

pulse 

0 

0 

* 

Actually  commanding 

r  -  p«T 

for  roll 

maneuver. 

.  V-;  k-:  irv  r 
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V.  Design  Responses 

This  chapter  presents  the  simulation  responses  for  the 
control  laws  developed  in  this  thesis.  Chapter  III  des¬ 
cribes  in  detail  the  process  used  in  designing  the  control 
laws,  and  Chapter  IV  discusses  the  maneuvers  considered. 
Appendix  C  lists  the  data  used  in  the  aircraft  models  for 
each  flight  condition.  The  four  flight  conditions  consider¬ 
ed  in  this  study  are:  transonic  flight  (0.9  Mach  at  20,000 
feet),  subsonic  flight  (0.6  Mach  at  30,000  feet),  supersonic 
flight  (1.6  Mach  at  30,000  feet),  and  landing  approach  (0.3 
Mach  at  30  feet) .  This  chapter  lists  the  resonses  for  the 
transonic  flight  condition  only.  The  simulation  responses 
for  the  remaining  three  flight  conditions  are  presented  in 
Appendix  D. 

The  main  concern  for  this  study  is  that  the  control 
surface  responses  are  within  the  rate  and  angular  deflection 
limits  listed  in  Table  3-1  for  both  the  healthy  and  failed 
right  elevator  models.  The  command  vector  elements,  listed 
in  Table  4-1,  are  therefore  commanded  with  magnitudes  that 
yield  maximum  control  surface  deflections  in  order  to  deter¬ 
mine  the  maximum  capabilities  of  the  aircraft.  In  some 
cases  large  roll  and  pitch  angles  are  generated  which  con¬ 
tradict  the  small  angle  assumptions  made  when  developing  the 
aircraft  perturbation  equations.  Although  this  invalidates 
the  linear  aircraft  models,  a  general  idea  of  the  control 
surface  deflections  needed  to  perform  the  maneuvers  is 
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obtained.  The  commanded  magnitudes  are  decreased  from  the 
maximum  values  if  a  large  angle  of  attack  or  sideslip  is 
generated. 

For  each  of  the  flight  conditions,  when  simulating  a 
longitudinal  maneuver  the  elements  of  the  diagonal  weighting 
matrix  I  corresponding  to  lateral  quantities  are  set  to  the 
values  used  in  the  sideforce  maneuver  design.  Similarly, 
when  simulating  a  lateral  maneuver  the  elements  of  £  corre¬ 
sponding  to  longitudinal  quantities  are  set  to  the  values 
used  in  the  g-command  maneuver  design.  This  is  done  because 
all  elements  of  £  affect  the  gain  matrices  KQ  and  K1  and 
because  the  sideforce  and  g-command  maneuvers  are  conven¬ 
tional  maneuvers. 

This  chapter  is  divided  into  seven  sections  correspond 
ing  to  each  maneuver  at  the  transonic  flight  condition. 

Each  section  includes  tables  that  list  the  design  parameters 
and  simulation  response  plots  generated  with  the  program 
MULTI.  Appendix  D  is  organized  in  a  similar  manner. 

Air-to-air  combat  often  begins  in  the  transonic  region 

of  the  flight  envelope.  Because  of  this,  the  transonic 

flight  condition  (0.9  Mach  at  20,000  feet)  is  an  important 

condition  for  demonstrating  the  maneuverability  of  the 

aircraft.  Because  of  a  relatively  high  dynamic  pressure 
_  2 

(q  =  552  lbs/ft  )  at  this  flight  condition,  the  control  sur¬ 
faces  are  extremely  effective  in  maneuvering  the  aircraft, 
and  designs  at  this  flight  condition  are  fairly  easy. 


V-l.  G-Command  Maneuver 

Tables  5-1  and  5-2  list  the  design  parameters  for  the 
g-command  maneuver  for  the  healthy  and  failed  right  elevator 
models.  As  shown  in  the  tables,  the  same  design  parameters, 
a,  E,  and  e  ,  are  used  for  both  aircraft  models. 

Figures  5-1  through  5-7  show  the  simulation  responses 
for  the  g-command  maneuver.  For  this  maneuver  the  healthy 
aircraft  longitudinal  responses  change  very  little  when  the 
right  horizontal  tail  fails.  Also,  the  healthy  aircraft 
lateral  responses  are  zero  and  are  very  small  with  the 
failed  model.  For  these  reasons  Figures  5-1  through  5-3 
give  simulation  responses  for  both  the  healthy  and  failed 
aircraft  models. 

The  maximum  normal  acceleration  commanded  (9g)  is 
limited  by  the  flaperon  deflection,  as  shown  in  Figures  5-4 
and  5-5.  A  comparison  of  Figures  5-4  and  5-5  shows  that, 
for  the  failure  case,  the  left  elevator  deflects  approxi¬ 
mately  twice  the  amount  required  when  both  tail  halves  are 
working.  Figure  5-5  indicates  that  the  flaperon  deflection 
is  slightly  asymmetric  for  the  failure  case.  This  asymmetry 
compensates  for  the  rolling  moment  created  by  the  left 
elevator  deflection.  Also,  Figure  5-7  shows  that  the  rudder 
and  canard  compensate  for  the  adverse  yaw  caused  by  a 
difference  in  drag  between  the  left  and  right  halves  of  the 


aircraft 


Table  5-1 


G-Command:  Healthy  Model,  0.9  Mach  at  20 


Sampling  Time:  T  =  0.02  second 


a  =  1.0 


c  =  1.0 


E  =  diag  (2.0,  0.1,  2.35,  1.0,  2.3,  1.0} 


«0 


0.0 

-0 . 2267E-03 

-0.3534E-01 

-0.1370E-02 

0 . 5821E-04 

0.0 

0 . 8425E-03 

-0.5018E-01 

0 . 8348E-03 

-0.3546E-04 

0.0 

0.0 

0.0 

0 . 6581E-02 

-0 . 701BE-02 

0.0 

0.0 

o 

o 

0.2548E-01 

0.1858E-02 

0.0 

0.0 

0.0 

0 . 2794E-01 

-0.1187E-02 

0. 1039E+01 

-0 . 1555E-02 

0 . 1944E+00 

0.8464E-03 

-0.3595E-04 

*1  =  «0 


Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =0.0 

A  =  9.0  g  (step) 

P 

q  =  0.3105  radian/ second  (step) 


000  Feet 


0 . 1407E-02 
-0 . 8572E-03 
-0.1200E-01 
-0.S830E-01 
-0.2870E-01 
-0.8691E-03 


Table  5-2 


G-Command:  Failed  Model,  0.9  Mach  at  20,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  1.0 
c  =  1.0 


Z  =  d 

liag 

{2.0,  0.1,  2 

.35,  1.0, 

2.3,  1. 

0} 

0.0 

-0.4529E-03 

-0.7062E-01 

-0.2750E-02 

0 . 1285E-C3 

C . 2833E-02 

0.0 

0.7775E-03 

-0.6031E-01 

0.8910E-02 

-0.9049E-02 

-0.1590E-01 

K  _  = 

0.0 

0.9072E-03 

-0.4008E-01 

-0.7235E-02 

0. 8970E-02 

0.1417E-01 

-0 

0.0 

-0.4018E-04 

-0 . 6264E-02 

0 . 2640E-01 

0. 6243E-03 

-0.6018E-01 

0.0 

-0.4012E-05 

-0.6256E-03 

0.2B04E-01 

-0.1310E-02 

-0 . 2888E-01 

0 . 1039E+01  -0.1554E-02 

0.1946E+00 

0.8576E-03 

-0 . 4007E-04 

-0. 8835E-03 

Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =0.0 

A  =  9.0  g  (step) 

np 

q  =  0.3105  radian/second  (step) 


ALPHA  (OfcC)  t  LATERAL  STATES 

.CO  2.00  4.00  6.  CO  6.00  tQ.OO 


Figure  5-3:  0.9M  Healthy/Failed  G-Command  —  Pitch 

Normal  Acceleration,  and  Lateral  Acceleration 


Figure  5-4:  0.9M  Healthy  G-Command  —  Elevator  and  Flap 

Deflections 


Figure  5-5:  0.9M  Failed  G-Command  —  Left  Elevator  and 

Flaperon  Deflections 


V-2.  Pitch  Pointing  Maneuver 


Tables  5-3  and  5-4  give  the  design  parameters  for  the 
pitch  pointing  maneuver.  As  with  the  g-command  maneuver  the 
design  parameters  a,  and  e  are  the  same  for  the  healthy 
and  failed  models. 

The  simulation  responses  for  the  pitch  pointing  maneu¬ 
ver  are  given  in  Figures  5-8  through  5-14.  For  the  same 
reasons  discussed  in  the  g-command  maneuver  section,  Figures 
5-8  through  5-10  give  simulation  responses  for  both  the 
healthy  and  failed  aircraft  models.  The  failure  model 
requires  that  the  left  elevator  must  deflect  approximately 
twice  the  angle  of  the  healthy  elevator  deflection.  Also, 
for  the  failure  model,  the  flaperons  deflect  asymmetrically 
to  counter  the  roll  induced  by  the  failed  right  horizontal 
tail,  and  the  rudder  and  canard  deflect  to  compensate  for 
adverse  yaw. 


RD-A151  908  MULTIVARIABLE  CONTROL  LAN  DESIGN  FOR  THE  AFTI/F-16  WITH 
A  FAILED  CONTROL  SURFACE(U)  AIR  FORCE  INST  OF  TECH 
MRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENG I.  .  R  A  ESLINGER 
UNCLASSIFIED  DEC  84  AFIT/GE/ENG/84D-28  F/G  1/2 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  196^  A 


Table  5-3 


Pitch  Pointing:  Healthy  Model,  0.9  Mach  at  20,000  Feet 

Sampling  Time:  T  =  0.02  second 

o  =  2.0 
e  =  1.0 

E  =  diag  (2.0,  3.0,  2.35,  1.0,  2.3,  1.0} 


0.0 

-0 . 6800E-02 

-0 . 3534E-01 

-0 . 1370E-02. 

0. 5821E-04 

0. 1407E-02 

0.0 

0.2527E-01 

-0.5018E-01 

0 . 8348E-03 

-0.3546E-04 

-0 . 8572E-03 

0.0 

0.0 

0.0 

0.6581E-02 

-0.7018E-02 

-0.1200E-01 

0.0 

0.0 

0.0 

0. 2548E-01 

0.1858E-02 

-0.5830E-01 

0.0 

0.0 

0.0 

0.2794E-01 

-0.11B7E-02 

-0 . 2870E-01 

0 . 1039E*01 

-0 . 4666E-01 

0. 1944E*00 

0 . 8464E-03 

-0.3595E-04 

-0.8691E-Q3 

0.0 

-0.1360E-01 

-0.7068E-01 

-0.2741E-02 

0 . 1164E-03 

0.2814E-02 

0.0 

0.5055E-01 

-0 . 1004E+00 

0.1670E-02 

-0 . 7092E-04 

-0.1714E-02 

0.0 

0.0 

0.0 

0.1316E-01 

-0.1404E-01 

-0. 2400E-01 

0.0 

0.0 

0.0 

0.5095E-01 

0.3715E-02 

-0. II66E+00 

0.0 

0.0 

0.0 

0.5589E-01 

-0.2374E-02 

-0.5739E-01 

0 . 2077E+01 

-0.9331E-01 

0.3888E»00 

0 . 1693E-02 

-0.7190E-04 

-0.1738E-02 

Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =  0.03491  radian/second  (1  second  pulse) 

A  =  0.0 

p  =0.0 

r  =0.0 
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Table  5-4 


V  VT.^V| 


Pitch  Pointing:  Failed  Model,  0.9  Mach  at  20,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  2.0 


c  =  1.0 

E  =  diag  {2.0,  3.0,  2.35,  1.0,  2.3,  1.0} 


0.0  -0.13S9E-02  -0.7062E-01  -0.2750E-02  0.1265E-03  0.2833E-02 
0.0  0.2332E-01  -0.6031E-01  0.8910E-02  -0.9049E-02  -0.1590E-01 
0.0  0.2722E-01  -0.4008E-01  -0.7235E-02  0.8970E-02  0.1417E-01 
0.0  -0 . 1205E-02  -0.6264E-02  0.2640E-01  0.6243E-03  -0.6018E-01 
0.0  -0.1204E-03  -0.6256E-03  0.2804E-01  -0.I31OE-02  -0.2B88E-01 
0. 1039E+01  -0.4662E-01  0.1946E*00  0.8S76E-03  -0.4007E-04  -0.8835E-03 


0.0  -0.2718E-01  -0.1412E+00  -0.5499E-02  0.2569E-03  0.5665E-02 
0.0  0.4665E-01  -0.1206E*00  0.1782E-01  -0.1810E-01  -0.3179E-01 
0.0  0 . S443E-0I  -0.8017E-01  -0.1447E-01  0.1794E-01  0.2834E-01 
0.0  -0.24UE-02  -0.1253E-0I  0.5280E-01  0.1249E -02  -0.1204E*00 
0.0  -0.2407E-03  -0.1251E-02  0.5607E-01  -0.2620E-02  -0.5777E-01 
0.2077E+01  -0.9324E-01  0.3892E»00  0.1715E-02  -0.8013E-04  -0.1767E-02 


Input  Ramp  Time:  0.2  second 


Command  Vector: 


u  =0.0 


Sr 

q  =  0.03491  radian/second  (1  second  pulse) 


A  =  0.0 

p  =0.0 

r  =0.0 


*  V*  .  O/dA.V.  .V.  v-  ■  .  -  .  -V -  .  -  .  '  .  •  «  .  •  .  •*.  -  .  ••  «  .  -  .  •  .  «“  **.  '  .  .*. 


kyy. 


-3.  Longitudinal  Translation  Maneuver 


Tables  5-5  and  5-6  give  the  design  parameters  for  the 
longitudinal  translation  maneuver.  Once  again  the  design 
parameters  a,  E_,  and  c  are  the  same  for  the  healthy  and 
failed  models. 

The  simulation  responses  for  the  longitudinal  transla¬ 
tion  maneuver  are  given  in  Figures  5-15  through  5-21.  As 
discussed  in  the  g-command  maneuver  section  Figures  5-15 
through  5-17  give  simulation  responses  for  both  the  healthy 
and  failed  aircraft  models.  The  maximum  longitudinal 
translation  command  of  0.9g  for  one  second,  which  yields  a 
normal  velocity  of  approximately  29  feet/second,  is  limited 
by  the  flaperon  deflection. 

As  with  the  previous  two  longitudinal  maneuvers  the 
failure  model  requires  that  the  left  elevator  must  deflect 
approximately  twice  the  angle  of  the  healthy  elevator 
deflection.  Also,  for  the  failure  model,  the  flaperons 
deflect  asymmetrically  to  counter  the  roll  induced  by  the 
failed  right  horizontal  tail,  and  the  rudder  and  canard 
deflect  to  compensate  for  adverse  yaw. 


Table  5-5 


Longitudinal  Translation:  Healthy  Model, 
0.9  Mach  at  20,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  1.0 
e  =  1.0 

I  =  diag  (2.0,  2.0,  2.35,  1.0,  2.3,  1.0} 


0.0 

-0 . 4533E-02 

-0 . 3534E-01 

-0.1370E-02 

0.5821E-04 

0.1407E-02 

0.0 

0 . 16B5E-01 

-0 . 5018E-01 

0. 8348E-03 

-0 . 3546E-04 

-0.8572E-03 

0.0 

0.0 

0.0 

0.6581E-02 

-0 . 7018E-02 

-0 . 1200E-0 1 

0.0 

0.0 

0.0 

0.2548E-01 

0 . 1858E-02 

-0 . 5830E-01 

0.0 

0.0 

0.0 

0.2794E-01 

-0 . 1187E-02 

-0.2870E-01 

0 . 1039E+01 

-0 . 3110E-01 

0 • 1944E+00 

0.8464E-03 

-0.3595E-04 

-0.8691E-03 

*1 


*0 


Input  Ramp  Time:  0.2  second 
Command  Vector: 
u  =0.0 

A  =  0.9  g  (1  second  pulse) 


Table  5-6 


Longitudinal  Translation:  Failed  Model, 
0.9  Mach  at  20,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  1.0 
e  =  1.0 


diag 

(2. 

CM 

O 

• 

<N 

O 

.35,  1.0 

,  2.3,  1. 

0} 

0.0 

-0.9059E-02 

-0.7062E-01 

-0.2750E-02 

0.1285E-03 

0.28  3  3E-02 

0.0 

0.1555E-01 

-0.6031E-01 

0.8910E-02 

-0.9049E-02 

-0. 1590E-01 

0.0 

0.1814E-01 

-0.4008E-01 

-0.7235E-02 

0.8970E-02 

0.141 7E-0 1 

0.0 

-0.8036E-03 

-0.6264E-02 

0.2640E-01 

0.6243E-03 

-0 . 6018E-01 

0.0 

-0.802SE-04 

-0.6256E-03 

0 . 2B04E-01 

-0.1310E-02 

-0 . 2888E-01 

0. 1039E+01 

-0.3I08E-01 

0 . 1946E*00 

0.8576E-03 

-0.4007E-04 

-0 . 8835E-03 

*1  '  «0 


Input  Ramp  Time:  0.2  second 
Command  Vector: 
u 


0.0 

0.9  g  (1  second  pulse) 


Figure  5-17:  0.9M  Healthy/Failed  Longitudinal  Translation 

—  Pitch  Rate,  Normal  Acceleration,  and  Lateral 
Acceleration 


r  r,  * -  •  •  "  “  V  -..  J 


V-4.  Roll  About  Velocity  Vector  Maneuver 

The  roll  maneuver  requires  some  modifications  to  the 
healthy  and  failure  models.  These  modifications  are  dis¬ 
cussed  in  Chapter  IV  and  are  included  in  Tables  5-7  and  5-8, 
which  give  the  design  parameters  a,  £_,  and  e  for  the  roll 
maneuver.  The  parameters  are  the  same  for  the  healthy  and 
failed  models. 

The  simulation  responses  for  the  roll  maneuver  are 
given  in  Figures  5-22  through  5-28.  For  this  maneuver  the 
longitudinal  responses  are  very  small  for  both  the  healthy 
and  failure  models.  Also,  the  healthy  aircraft  lateral 
responses  change  very  little  when  the  right  horizontal  tail 
fails.  For  these  reasons  Figures  5-22  through  5-24  give 
simulation  responses  for  both  the  healthy  and  failed 
aircraft  models. 

The  roll  rate  commanded  is  limited  to  308  degrees  per 
second  for  three  seconds  in  order  to  avoid  generating  too 
large  a  sideslip  angle.  The  large  roll  rate  command  creates 
a  steady-state  roll  angle  of  approximately  1000  degrees 
which  corresponds  to  rolling  the  aircraft  2.8  times.  This 
large  roll  angle  contradicts  the  small  angle  assumptions 
made  when  developing  the  linear  aircraft  models;  however, 
approximations  of  the  surface  deflections  required  to  obtain 
the  commanded  roll  rate  are  obtained.  This  maneuver  should 
be  simulated  with  a  full  nonlinear  model  in  the  future. 

The  aileron  deflection  for  the  healthy  aircraft  (Figure 
5-27)  is  less  than  the  flaperon  deflections  for  the  failure 
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case  (Figure  5-26)  because  the  aileron  deflection  includes  a 
contribution  from  the  horizontal  tail.  When  the  right  hori¬ 
zontal  tail  fails  the  flaperons  must  deflect  more  to  achieve 
the  commanded  roll  rate.  The  rudder  and  canard  deflections 
are  unaffected  by  the  failure  of  the  right  horizontal  tail. 


Table  5-7 


Roll  About  Velocity  Vector: 

Healthy 

Model , 

0. 

9  Mach  at 

20,000 

Feet 

Sampling 

Time:  T  =  0 

.02  second 

a  = 

1.0 

e  = 

1.0 

|M 

li 

diag 

(2.0,  1.0,  2 

.35,  1.0, 

H ‘ 

• 

O 

** 

H* 

* 

0} 

0.0 

-0.2267E-03 

-0.3534E-01  - 

■0.1370E-02 

0.2531E-04 

0.1407E-02 

0.0 

0.8425E-03 

-0.5018E-01 

0 . 8348E-03 

-0.1542E-04 

-0 . 8572E-03 

0.0 

0.0 

0.0 

0.6581E-02 

-0.3051E-02 

-0 . 1200E-01 

*0 

0.0 

0.0 

0.0 

0.2548E-01 

0 . 8077E-03 

-0.5830E-01 

0.0 

0.0 

0.0 

0.2794E-01 

-0.S161E-03 

-0.2870E-01 

0 . 1039E+01  -0.1555E-02 

0 . 1944E+00 

0. 8464E-03 

-0.1563E-04 

-0.8691E-03 

*1  =  «0 

Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =0.0 

A  =  0.0 

p  =  5.376  radians/second  (3  second  pulse) 


Table  5-8 


Roll  About  Velocity  Vector:  Failed  Model, 
0.9  Mach  at  20,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  1.0 
e  =  1.0 

Z  =  diag  {2.0,  1.0,  2.35,  1.0,  1.0,  1.0} 


-0.4529E-03 
0 . 7775E-03 
0 . 9072E-03 
-0.4018E-04 
-0.4012E-05 
-0.1554E-02 


-0.7062E-01 
-Q.6031E-01 
-0.4008E-01 
-0.6264E-02 
-0.6256E-03 
0 . 1946E+00 


-0 . 2750E-02 
0 . 8910E-02 
-0.7235E-02 
0 . 2640E-01 
0.2B04E-01 
0.8576E-03 


0.55B5E-04 

-0.3934E-02 

0.3900E-02 

0.2715E-03 

-0.5695E-03 

-0.1742E-04 


Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =0.0 

An  =  0.0 

P 

q  =  0.0 

A  =  0.0 

p  =  5.376  radians/second  (3  second  pulse) 
r-aT  =0.0 
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Figure  5-22:  0.9M  Healthy/Failed  Roll  —  Sideslip  Angle 

Pitch  Angle,  Angle  of  Attack,  and  Forward  Velocity 


r,rE.  SECONDS 


Figure  5-23:  0.9M  Healthy/Failed  Roll  —  Roll  Angle,  Roll 

Rate,  and  Yaw  Rate 
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Figure  5-25:  0.9M  Healthy  Roll  —  Elevator  and  Flap 

Deflections  and  Thrust 


V-5.  Sidef orce  (Flat  Turn)  Maneuver 

Tables  5-9  and  5-10  give  the  design  parameters  a  , 
and  e  for  both  the  healthy  and  failed  cases  of  the  sideforce 
maneuver.  The  same  parameters  are  used  for  both  cases. 

Figures  5-29  through  5-35  present  the  simulation 
responses  for  a  0.8g  sideforce  command.  The  maximum  value 
commanded  is  limited  by  the  canard  deflection.  As  with  the 
roll  maneuver  the  healthy  and  failed  longitudinal  responses 
are  very  small,  and  the  lateral  responses  change  very  little 
when  the  right  horizontal  tail  fails;  therefore,  Figures 
5-29  through  5-31  give  simulation  responses  for  both  the 
healthy  and  failed  models. 

For  this  maneuver  the  horizontal  tail  is  mainly  used  to 
counter  the  longitudinal  pitching  moment  produced  by  the  ca¬ 
nard  deflection  although  it  also  contributes  to  roll  control 
for  the  healthy  aircraft  case.  When  the  right  horizontal 
tail  fails  the  left  horizontal  tail  deflects  approximately 
twice  the  deflection  of  the  healthy  elevators.  The  flap- 
erons  deflect  slightly  asymmetrically  to  compensate  for  the 
roll  created  by  the  failed  right  horizontal  tail.  The 
rudder  and  canard  deflections  change  only  slightly  for  the 


failure  case. 


Table  5-9 

Sideforce  Maneuver:  Healthy  Model, 
0.9  Mach  at  20,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  1.0 
e  =  l.O 


diag 

{2. 

0,  1.0,  2 

.35,  1.0 

,  2.3,  1. 

0} 

o 

Li 

-0 . 2267E-03 

-0 . 3534E-01 

-0 . 1370E-02 

0.5821E-04 

0.1407E-02 

o 

o 

0.B425E-03 

-0.5018E-01 

0 . 8  348E-03 

-0.3546E-04 

-0.8572E-03 

o 

o 

0.0 

0.0 

0.6581E-02 

-0.7018E-02 

-0.1200E-01 

0.0 

o 

o 

0.0 

0 . 2548E-01 

0.1858E-02 

-0.5830E-01 

0.0 

0.0 

0.0 

0 . 2794E-01 

-0.11B7E-02 

-0.2870E-01 

jo. 1039E+01 

-0 . 1555E-02 

0. 1944E+00 

0 . 8464E-03 

-0.3595E-04 

-0. 8691E-03 

*1  =  ^0 

Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =0.0 

A  =  0.8  g  (step) 

YP 

p  =0.0 

r  =  0.02760  radian/second  (step) 
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Table  5-10 


Sideforce  Maneuver:  Failed  Model, 
0.9  Mach  at  20,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  l.O 
e  =  l.O 

Z  =  diag  (2.0,  1.0,  2.35,  1.0,  2.3,  1.0} 


Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

A  =0.0 

P 

q  =0.0 

A  =  0.8  g  (step) 
yP 

p  =0.0 

=  0.02760  radian/second  (step) 


r 


c.cc 


•  .cc 


6 .  cc 


°  cc 


t.CC  o.CC  *.CC  i  cc 
_ rrr£.  ScCCNGS 

Figure  5-29:  0.9M  Healthy/Failed  Sideforce  —  Sideslip 

Angle,  Pitch  Angle,  Angle  of  Attack,  and  Forward  Velocity 


Figure  5-30:  0.9M  Healthy/Failed  Sideforce  —  Roll  Angle 

and  Roll  Rate 


Figure  5-31:  0.9M  Healthy/Failed  Sideforce  —  Lateral 

Acceleration  and  Yaw  Rate 
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V-6.  Yaw  Pointing  Maneuver 

The  design  parameters  <x ,  E,  and  c  for  the  yaw  pointing 
maneuver  are  given  in  Tables  5-11  and  5-12  for  both  the 
healthy  and  failed  models.  The  parameters  are  the  same  for 
both  cases. 

Figures  5-36  through  5-42  present  the  simulation 
responses  for  a  four  degree  yaw  pointing  command.  The 
maximum  angle  commanded  for  this  maneuver  is  limited  by  the 
ability  to  hold  the  roll  angle  to  a  small  value.  As  with 
the  roll  and  sideforce  maneuvers.  Figures  5-36  through  5-38 
give  simulation  responses  for  both  the  healthy  and  failed 
aircraft  models. 

For  the  yaw  pointing  maneuver  the  horizontal  tail  is 
used  for  both  pitch  and  roll  control.  When  the  right  hori¬ 
zontal  tail  fails  the  flaperons  take  over  the  roll  control 
and  deflect  asymmetrically  to  counter  the  rolling  moment  of 
the  left  horizontal  tail,  which  deflects  for  pitch  control. 
The  rudder  and  canard  responses  are  not  affected  much  by  the 
failure  for  this  maneuver. 


Table  5-11 


Yaw  Pointing:  Healthy  Model,  0.9  Mach  at  20,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  6.0 


e  =  1.0 


I  =  diag  {2.0,  1.0,  2.35,  3.0,  3.25,  3.0) 


— 0 


0.0 

-0.2267E-03 

-0.3534E-01 

-0.411 IE-02 

O.0225E-O4 

0 . 4222E-02 

0.0 

0.8425E-03 

-0 . 501 8E-01 

0 . 2504E-02 

-0 . 5011E-04 

-0.2572E-02 

o 

o 

0.0 

o 

o 

0.1974E-01 

-0.9916E-02 

-0 . 3600E-0 1 

o 

o 

o 

o 

0.0 

0 . 7643E-01 

0 . 2625E-02 

-0 . 1 749E-*  00 

o 

o 

o 

o 

o 

o 

0 . 8383E-01 

-0 . 1677E-02 

-0 . 8609E-0 1 

0 . 1039E+01 

-0.1555E-02 

0 . 1944E+00 

0 . 2539E-02 

-0.5080E-04 

-0.2607E-02 

K 


1 


0.0 

-0.1360E-02 

-0.2120E*00 

-0.2467E-01 

0.4935E-01 

0.2533E-01 

0.0 

0.5055E-02 

-0 . 3011E+00 

0.1503E-01 

-0.3006E-03 

-0.1543E-01 

o 

O 

o 

o 

0.0 

0.11B5E*00 

-0.5950E-01 

-0 . 2160E*00 

0.0 

0.0 

0.0 

0 . 4586E*00 

0. 1575E-01 

-0.1049E*01 

0.0 

0.0 

0.0 

0 . 5030E+00 

-0.1006E-01 

-0.5165E*00 

0.6232E+01 

•0.9331E-02 

0 . 1166E+01 

0. 1523E-01 

-0. 3048E-03 

-0. 1564E-01 

Input  Ramp  Time:  0.2  second 
Command  Vector: 
u  =  0.0 

An  =  0.0 


P 

p  =  0.0 

r  =  -0.05818  radian/second  (1  second  pulse) 


Table  5-12 


Yaw  Pointing:  Failed  Model,  0.9  Mach  at  20,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  6.0 
e  =  1.0 

Z  =  diag  {2.0,  1.0,  2.35,  3.0,  3.25,  3.0} 


-0 . 2718E-02 
0 . 4665E-02 
0.5443E-02 
-0 . 2411E-03 
-0.2407E-04 
-0.9324E-02 


-0.4237E*00 
-0.3619E+00 
-0.2405E*00 
-0.3759E-01 
-0.3754E-02 
0. 2169E+01 


-0.4949E-0X 
0.1604E*00 
-0.1302E*00 
0.4752E*00 
0. 5046E+00 
0.1544E-01 


Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =  0.0 

A  =  0.0 

np 

q  =  0.0 


p  =  0.0 

r  =  -0.05818  radian/second  (1  second  pulse) 


Table  5-13 


Lateral  Translation:  Healthy  Model, 
0.9  Mach  at  20,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  3.0 


diag 

{2.0 

,  l.o,  2 

.35,  2.5, 

3.0,  3. 

25} 

___ 

0.0 

-0.2267E-03 

-0.3534E-01 

-0 . 3426E-02 

0.7593E-04 

0 . 4573E-02 

0.0 

0 . S425E-03 

-0.5018E-01 

0 . 2087E-02 

-0 . 4625E-Q4 

-0 . 2786E-02 

0.0 

0.0 

0.0 

0 . 1645E-01 

-0 . 9 154E-02 

-0.3900E-01 

0.0 

0.0 

0.0 

0.6369E-01 

0.2423E-02 

-0 . 1B95E+00 

0.0 

0.0 

0.0 

0 . 6986E-01 

-0 . 1548E-02 

-0.9326E-01 

I  0 . 1039E+O1 

-0.1555E-02 

0. 1944E+00 

0. 2116E-02 

-0 . 4689E-04 

-0. 2825E-02 

i  '• 


I 


"~0.0 

-0 . 6800E-03 

-0.1060E.00 

-0.1028E-01 

0.2278E-03 

0. 1372E-01™ 

0.0 

0.2527E-02 

-0.1505E-00 

0 . 6261E-02 

-0.1388E-03 

-0.8358E-02 

0.0 

0.0 

0.0 

0 . 4936E-01 

-0.2746E-01 

-0 . 1170E*00 

0.0 

0.0 

0.0 

0. 1911E-00 

0 . 7269E-02 

-0. 5685E-00 

•  •'  * 

0.0 

0.0 

0.0 

0.2096E*00 

-0.4645E-02 

-0.2798E.00 

•  4 

.  4 

0. 3116E+01 

-0 . 4666E-02 

0 . 5832E+00 

0 . 6348E-02 

-0.1407E-03 

-0 . 8474E-02 

.  - 

— — 

'  *-  > 

Input  Ramp  Time:  0.3  second 
Command  Vector: 


u  =0.0 


A  =0.0 
nP 

q  =0.0 


P 


0.6  g  (1  second  pulse) 


0.0 


r  =0.0 
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Table  5-14 


Lateral  Translation:  Failed  Model, 
0.9  Mach  at  20,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  3.0 

e  =  1.0 

Z  =  diag  (2.0,  1.0,  2.35,  2.5,  3.0,  3.25} 


o.o 

0.0 

0.0 

0.0 

0.0 

0.1039E*01 


-0.4529E-03 
0.7775E-03 
0 . 9072E-03 
-0 . 4018E-04 
-0.4012E-05 
-0.1554E-02 


-0 . 7062E-01 
-0.6031E-01 
-0.4008E-01 
-0.6264E-02 
-0 . 6256E-03 
0 . 1946E+00 


-0.6874E-02 
0.2227E-01 
-0.1809E-01 
0.6600E-01 
0 . 7009E-01 
0. 2144E-02 


0 . 1676E-03 
-0.1180E-01 
0.1170E-01 
0.8144E-03 
-0.1708E-02 
-0.5226E-04 


0.9206E-02 
-0 . 5166E-01 
0 . 4605E-01 
-0.1956E-00 
-0.9387E-01 
-0.2871E-02 


“I 


— - 

0.0 

-0.1359E-02 

-0 . 2119E*00 

-0.2062E-01 

0 . 5027E-03 

0 . 2  762E-01 

0.0 

0 . 2332E-02 

-0 . 1809E+00 

0.6682E-01 

-0.3541E-01 

-0.1550E*00 

0.0 

0.2722E-02 

-0 . 1203E+00 

-0.5426E-01 

0.3510E-01 

0 . 1382E*00 

0.0 

-0.1205E-03 

-0.1879E-01 

0. 1980E*00 

0.2443E-02 

-0.5868E*00 

0.0 

-0.1204E-04 

-0 . 1877E-02 

0. 2103E+00 

-0 . 5125E-02 

-0. 2816E*00 

0.3116E+01 

-0.4662E-02 

0. 5838E+01 

0.6432E-02 

-0.1568E-03 

-0.8614E-02 

Input  Ramp  Time:  0.3  second 

Command  Vector: 

u  =0.0 

A  =  0.0 
n 

P 

q  =0.0 

A  =  0.6  g  (1  second  pulse) 

yp 

p  =0.0 
r  =0.0 
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Figure  5-43:  0.9M  Healthy/Failed  Lateral  Translation  — 

Sideslip  Angle,  Pitch  Angle,  Angle  of  Attack,  and  Forward 
Velocity 


Figure  5-44:  0.9M  Healthy/Failed  Lateral  Translation 

Roll  Angle,  Roll  Rate,  and  Yaw  Rate 


Figure  5-45:  0.9M  Healthy/Failed  Lateral  Translation 

Normal  Acceleration  and  Lateral  Acceleration 
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Appendix  D  lists  the  remainder  of  the  simulation  re¬ 
sponses  generated  in  this  study.  As  shown  in  this  chapter 
the  transonic  flight  condition  simulation  responses  show 
that  the  aircraft  is  still  able  to  perform  the  commanded 
maneuvers  when  the  right  horizontal  tail  is  failed. 

Appendix  D  presents  some  cases  where  the  aircraft  perform¬ 
ance  is  degraded  when  the  right  horizontal  tail  fails. 

The  next  chapter  discusses  the  results  of  this  study  in 
general  terms  and  presents  recomendations  for  further  study. 


Conclusions  and  Recommendations 


In  this  study  control  laws  are  developed  for  the 
AFTI/F-16  using  the  Porter  method.  Acceptable  designs  are 
achieved  for  a  model  of  the  healthy  aircraft  with  all 
control  surfaces  operational  and  for  a  model  with  the  right 
horizontal  tail  free-floating.  Designs  are  accomplished  for 
seven  maneuvers  at  four  key  flight  conditions  within  the 
aircraft's  flight  envelope.  The  control  laws  consist  of  the 
gain  matrixes  and  Kq  which  can  be  implemented  in  a  digit¬ 
al  flight  control  computer.  Each  control  law  is  tailored  to 
a  specific  maneuver  at  a  particular  flight  condition.  The 
maneuvers  include  both  conventional  and  CCV  maneuvers. 


VI-1.  Design  Summary 


Control  laws  are  first  developed  for  the  healthy  air¬ 
craft  model  with  the  aid  of  the  computer  program  MULTI, 
which  proves  invaluable  when  utilizing  the  Porter  method. 

The  transonic  and  supersonic  flight  conditions  (0.9  Mach  at 
20,000  and  1.6  Mach  at  30,000)  are  the  simplest  cases  to 
design  because  of  the  high  dynamic  pressure  at  these  flight 
conditions . 

Maximum  inputs  are  commanded  for  each  maneuver  in  order 
to  estimate  the  control  surface  deflections  required  for 
maximum  performance.  Because  commanding  maximum  maneuvers 
invalidates  the  small  angle  assumptions  made  in  deriving  the 
linearized  perturbation  models,  some  of  the  simulation 
responses  are  open  to  interpertation.  A  nonlinear  simula- 
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tion  of  the  designs  would  eliminate  any  problems  due  to  the 
linearization  assumptions. 

After  acceptable  control  laws  are  achieved  for  the 
healthy  aircraft  model,  the  same  design  parameters  (M,  a,  E, 
and  e )  are  used  to  calulate  gain  matrices  for  the  model  with 
a  failed  right  horizontal  tail.  The  calculation  of  new  gain 
matrices  implies  that  the  system  must  be  able  to  detect  the 
failure  and  reconfigure  accordingly.  The  detection  problem 
is  beyond  the  scope  of  this  study,  but  it  is  assumed  that  a 
failure  detection  scheme  can  be  developed. 

For  most  cases  considered  in  this  thesis,  satisfactory 
simulation  responses  are  achieved  for  the  failure  model 
without  requiring  further  design  trials.  The  four  cases 
that  do  require  changes  in  the  design  parameters  or  input 
command  magnitudes  are:  pitch  pointing  at  supersonic  speed; 
longitudinal  translation  at  supersonic  speed;  roll  at 
supersonic  speed;  and  roll  at  landing  approach.  In  each  of 
these  four  cases  the  failure  of  the  right  horizontal  tail 
requires  an  increased  deflection  of  the  remaining  control 
surfaces  to  compensate  for  the  loss.  For  these  cases,  if 
changing  the  design  parameters  does  not  yield  a  satisfactory 
result,  then  the  maximum  maneuver  commanded  is  decreased. 

In  general,  when  the  right  horizontal  tail  fails,  the 
left  horizontal  tail  assumes  primary  pitch  control  and  the 
flaperons  take  over  complete  roll  control.  The  flaperons 
also  compensate  for  unwanted  rolling  and  pitching  moments 
created  by  the  left  horizontal  tail  deflection.  The 
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deflection  of  the  left  horizontal  tail  also  creates  a  yawing 
moment  due  to  the  difference  in  drag  between  the  left  and 
right  halves  of  the  aircraft.  The  rudder  and  canards  com¬ 
pensate  for  this  yawing  moment. 

The  simulation  results  of  this  study  indicate,  for  the 
linearized  models  used  here,  that  the  aircraft  is  capable  of 
maneuvering  with  a  free-floating  right  horizontal  tail.  Al¬ 
though  much  work  must  still  be  accomplished  before  these 
results  are  accepted  in  practice,  the  results  do  imply  that 
it  is  worthwhile  to  continue  research  into  the  effects  of  a 
failed  right  horizontal  tail.  The  procedures  used  in  this 
study  are  applicable  to  other  control  surface  failures,  such 
as  a  failed  flaperon. 


VI-2.  Recommendations  for  Future  Study 

Although  the  computer  program  MULTI  has  proved  to  be  a 
great  aid  in  performing  multivariable  designs  with  the 
Porter  method,  additional  improvements  can  still  be  made  in 
the  program.  As  pointed  out  in  Chapter  III,  the  factoring 
routine  in  MULTI  Option  #6  becomes  error  prone  when  large 
order  transfer  functions  are  calculated.  Correction  of  this 
problem  will  allow  the  user  to  check  stabilty  without  run¬ 
ning  several  time  consuming  simulation  runs. 

Even  if  MULTI  Option  #6  had  the  desired  increased 
factoring  accuracy,  it  is  still  desirable  at  times  to  check 
simulation  responses  over  long  periods  of  time.  Currently 
the  plotting  routine  in  MULTI  does  not  generate  enough  data 


points  for  plotting  long  simulation  times.  If  more  data 
points  are  plotted  then  more  details  can  be  observed  on  the 
response  plots,  which  may  indicate  limit  cycles  or  slowly 
growing  instabilities.  This  change  should  be  readily 
achievable. 

The  addition  of  a  nonlinear  simulation  block  to  MULTI 
is  strongly  recommended.  Such  an  addition  would  allow  for 
more  realistic  simulations  of  the  control  law  designs.  A 
nonlinear  simulaion  on  MULTI  is  feasible  because  the  differ¬ 
ential  equation  solver  subroutine  (called  ODE)  currently 
used  with  MULTI  is  capable  of  handling  nonlinear  equations. 

Other  extentions  of  this  thesis  can  be  studied  with 
MULTI  in  its  current  state.  To  add  more  realism  to  the 
simulations,  a  computational  time  delay  and  sensor  models 
should  be  added  to  the  models  of  this  study.  Other  phenom¬ 
ena,  like  the  effects  of  wind  gusts  and  sensor  noise,  should 
also  be  studied. 

The  results  of  this  thesis  are  a  small  step  toward  the 
goal  of  designing  flight  control  systems  that  can  compensate 
for  battle  damage  or  control  surface  failures.  If,  however, 
the  control  laws  developed  in  this  study  are  simulated  under 
more  realistic  conditions  and  fine  tuned,  the  results  could 
make  a  significant  contribution  to  the  self-repairing  flight 
control  program  undertaken  by  the  Air  Force  Flight  Dynamics 
Laboratory. 
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Appendix  A:  Multivariable  Control  Theory 

This  thesis  utilizes  the  multivariable  design  method  of 
Professor  Brian  Porter  of  the  University  of  Salford,  England 
(9).  The  design  method  uses  output  feedback  with  high-gain 
error-actuated  controllers.  Output  feedback  is  advantagous 
since  state  variables  may  be  difficult  to  measure  while 
system  response  data  are  more  readily  available. 

System  State  Equations 

Porter's  method  works  equally  well  for  either  contin¬ 
uous  or  discrete  systems,  but  it  is  often  easier  to  first 
examine  a  system  in  the  continous  time  domain.  This  is 
because  designs  can  be  more  easily  visualized  in  the  s-plane 
than  in  the  z-plane.  A  continuous  time  system  is  represent¬ 
ed  by  the  state  space  model: 

• 

x  =  Ax  +  Bu 

X  =  Cx  (A-l) 


where 

A  =  continuous  plant  matrix  (n  x  n) 

B  =  continuous  input  control  matrix  (n  x  m) 
C  =  continuous  output  matrix  (£  x  n) 
x  =  state  variable  vector  with  n  states 


u  =  input  vector  with  m  inputs 
Y_  =  output  vector  with  £  outputs 


The  system  inputs  for  an  aircraft  are  the  control  surface 
deflections  or  actuator  input  commands,  and  the  system 
outputs  are  aircraft  responses  affected  by  the  inputs. 

The  method  does  not  allow  for  a  feedforward,  D,  matrix. 
If  such  a  matrix  is  present  in  the  original  state  space 
model,  the  control  inputs  must  be  redefined  as  states  so 
that  the  D  matrix  is  absorbed  into  the  C  matrix.  This  can  be 
accomplished  by  incorporating  the  actuator  dynamics  into  the 
plant  model.  Actuator  inputs  then  become  control  inputs. 

To  employ  Porter's  method,  it  is  desirable  (but  not 
necessary)  to  partition  the  system  state  equations  as 
follows : 


The  equations  are  partitioned  so  that  B 2  and  C 2  are  square 
(m  x  m)  and  (  l  x  & )  matrices,  respectively.  The  method 
requires  that  the  number  of  inputs  to  the  system  equals  the 
number  of  outputs  which  means  m  =  i  ,  and  therefore  the 
dimension  of  B 2  equals  the  dimension  of  C^.  It  is  often 
possible  to  form  the  state  equations  so  that  B^  =  0^. 

For  the  discrete  case  the  system  equations  are  written 
as  follows: 


A-2 


x [ (k+1 ) T]  =  $x(kT)  +  j/u(kT) 
£(kT)  =  rx(kT) 


where 


<J>  =  exp  (AT)  =  discrete  plant  matrix 


'I'  =  f  exp(AT)Bdt  =  discrete  input  control  matrix 

J  0 

r  =  c  =  discrete  output  matrix 


In  the  above  equations  T  is  the  sampling  period,  and  k 
takes  on  integer  values  from  zero  to  plus  infinity. 


System  With  Output  Feedback 


Figure  A-l  shows  the  block  diagram  for  a  continuous 
output  feedback  system,  where  v  is  the  command  input  vector, 
and  £  is  the  desired  output  vector.  The  blocks  for  the 
plant  are  derived  directly  from  the  system  state  equations, 
Equation  (A-l).  The  proportional  plus  integral  controller 
has  three  parameters,  Kq,  and  g,  which  must  be  deter¬ 

mined  by  the  designer.  The  output  signal  of  the  controller, 
u,  is  given  in  the  following  control  law  equation: 


u  =  g(KQe  +/ K x  edt) 


(A-4) 


where 


u  is  the  output  signal  of  the  controller 
e  is  the  error  signal  at  the  input  of  the  controller 
Kq  is  the  proportional  gain  matrix 


is  the  gain  matrix  for  the  integral  term 
g  is  the  scalar  forward  path  gain 

A  measurement  matrix  M  is  included  in  the  system  if  the 
plant  is  irregular.  Regular  and  irregular  plants  are  dis¬ 
cussed  later. 

The  discrete  system  block  diagram,  shown  in  Figure  A-2, 
is  similar  to  the  continuous  system,  but  Equation  (A-4) 
becomes 

u(kT)  =  (1/T) [KQe(kT)  +K1z(kT)]  (A-5) 

where  the  forward  path  gain  g  equals  the  sampling  frequency, 
(1/T).  The  zJkT)  matrix  is  derived  from  the  backward 
difference  equation, 

z[(k+l)T]  =  z(kT)  +  Te(kT)  (A-6) 

The  steps  to  be  taken  next  in  the  design  method  depend 
on  whether  or  not  [CB]  has  full  rank,  i.e.,  does  it  have  an 
inverse.  If  the  matrix  [CB]  has  full  rank,  the  plant  is 
called  "regular"  and  no  measurement  matrix  M  is  needed. 
However,  if  [CB]  does  not  have  full  rank,  the  plant  is 
called  "irregular"  and  M  is  needed  to  form  a  new  matrix  [FB] 
(see  Equations  (A-12)  through  (A-14))  which  does  have  an 
inverse.  This  is  explained  in  more  detail  in  the  next 
sections.  When  the  partitioned  B  matrix  in  Equation  (A-2) 


has  the  form 


Regular  Plant 


For  the  system  to  be  classified  as  "regular"  the  first 
Markov  parameter  [CB]  must  have  full  rank.  If  this  is  true 
the  gain  matricies  can  be  found  from 

Kq  «  [CB] ~1Z  (A-10) 

and 

Kx  =  o[CB]-1I  ( A-ll ) 

where 

a  is  a  constant  which  assigns  the  ratio  of  proportional 
to  integral  control 

I  is  the  diagonal  weighting  matrix 

The  diagonal  weighting  matrix,  Z_  =  diag  (o^,  ...  ,  , 

is  specified  by  the  designer.  Each  (i  =  l,...,l)  deter¬ 
mines  the  weighting  of  the  effect  of  a  particular  error 
signal  on  each  control  input. 


If  the  first  Markov  parameter  [CB]  is  rank  deficient, 
then  the  plant  is  called  "irregular".  In  this  case,  the  C 
matrix  must  be  replaced  by 


F  =  [F1 *  F2]  (A-12) 

where 

Fx  =  [Cx  +  Mul  (A- 13 ) 

F2  =  [C2  +  Mi2]  (A-14) 

The  matrix  M  in  the  above  equations  is  a  measurement  matrix 
which  is  chosen  such  that  the  matrix  [FB]  has  full  rank. 

The  designer  chooses  the  measurement  matrix  so  that  it  is  as 
sparse  as  possible,  thus  the  smallest  number  of  additional 
measurements  are  required.  Reference  12  gives  an  approach 
for  selecting  the  measurement  matrix  to  achieve  optimal 
decoupling.  Once  M  is  formed,  Kq  and  K1  are  computed  by 

Ho  =  ( A-15 ) 

Hi  =  a[FB]  ~h  ( A-16 ) 

which  are  similar  to  Equations  (A-10)  and  (A-ll) . 

For  irregular  plants  the  error  vector  e  is  defined  as 

e  =  v  -  w  (A- 17) 


where 


For  step  inputs  the  values  of  the  rates,  x^ ,  become  zero  in 
the  steady  state  because  they  represent  kinematic  variables. 

Asymptotic  Characteristics 

As  the  gain  factor  of  the  system,  g  (or  1/T  for  the 
discrete  case) ,  approaches  infinity,  the  system  transfer 
function  matrix  G(s)  assumes  the  asymptotic  form 

r  (x)  =  ru)  +  r  (X)  <A-i9) 

where 

r_(X)  is  the  slow  transfer  function  matrix 
f ( X )  is  the  fast  transfer  function  matrix 


The  roots  of  the  asymptotic  closed-loop  transfer  func¬ 
tion  may  be  grouped  into  three  sets:  ,  Z_2 ,  and  Z^.  Table 

A-l  gives  the  equations  for  finding  these  asymptotic  roots. 
Sets  and  Z_2  correspond  to  the  slow  modes  of  the  system, 
where  the  modes  associated  with  the  roots  in  Z^  become 
uncontrollable,  and,  for  regular  plants,  the  modes  associat¬ 
ed  with  the  roots  in  Z_2  become  unobservable  as  the  gain 
increases.  Set  >  the  infinite  roots,  are  associated  with 
the  fast  modes  of  the  system  which  become  dominant  as  the 
gain  increases. 

The  roots  in  set  Z^  correspond  to  the  transmission 
zeros  of  the  system  which  are  not  altered  by  output  feed¬ 
back.  Since  infinite  gain  cannot  be  implemented,  and  is  not 


desirable,  the  closed-loop  roots  of  the  system  tend  to 
migrate  toward  the  transmission  zeros.  This  may  adversely 
affect  the  system  stability  if  the  location  of  these  zeros 
is  in  the  unstable  region.  Reference  11  gives  a  procedure 
for  locating  the  transmission  zeros  of  a  system. 

As  the  gain  increases  the  system  output  responses 
become  increasingly  decoupled.  The  asymptotic  closed-loop 
transfer  function  for  the  continuous  case  has  the  form 


gcr  -i  go  ? 

MX)  =  diag  {  - ,  - ,..., 

X+gc^  x+g(J2 


(A— 20) 


For  the  discrete  case  the  form  is 


o  i  o  2 

MX)  =  diag  {  - ,  - ,..., 

X-l+Oj^  X-l+o2 


°9 

- * - }  ( A-21 ) 


where  the  (i  =  1 are  determined  by  the  weighting 
matrix,  £. 

The  computer  program  MULTI  is  of  great  value  in 


reducing  the  time  required  to  achieve  a  satisfactory  design. 
The  MULTI  User's  Manual  (8)  describes  the  program  and  its 
operation. 


Control  Surface  Derivatives  —  Healthy  Aircraft 
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(B-23) 


X '  =  ^  (C  ), 

6ft  m  x5f  b 


(B-24 ) 
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(B-27 ) 


(3-28) 


~  (cm  ), 

lY  m5c  b 


(B-29) 


Note:  Reference  1  is  the  source  of  information  for  (C  ) .. 

m.  b 
oc 

The  rigid  data  points  (2:4-55  -  4-123)  are  averaged  for 
canard  deflections  of  15  and  25  degrees  and  then  flex- 
iblized  with  altitude  (2:9-17). 
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Stability  to  Body  Axis  Conversions 


Note:  In  the  following  equations,  a  6  by  itself 

indicates  that  the  conversion  is  the  same  for  all  control 
surface  derivatives. 


(CXu)b  =  ’  (CDu  +  2CD)  C°s2aT  "  (CLa+CD)  Sin2(XT 

+  (CD  +CL  +CL)  cosaT  sir,aT  (B-57) 

a  u 


ICX  'b  -  (CL-CD  1  °os2»T  *  ,CL  *2CL>  sil,2“T 

a  a  u 


f  (CT  -Cn  -C„)  cosar„  sina, 


(B-58 ] 


(cn  )b  =  Cn  cos  aT  +  sin  aT 

r  r  p 

+  (Ct  +Cn  )  sinaT  cosoT 
r  P 

(Cn5»b  ■  %  cosoT  *  %  sinaT 

Note:  A  question  was  raised  late  in  the  study  about 

consistency  of  units  in  the  longitudinal  stability  to  body 

axis  conversion  equations  that  involve  the  terms  CL  ,  CD  , 

u  u 

and  C  .  All  of  the  body  axis  aerodynamic  coefficients 
u 

calculated  with  these  equations  should  be  dimensionless; 

however,  the  stabiltiy  axis  aerodynamic  coefficients  CL  , 

u 

,  and  C  listed  in  Appendix  C  have  units  of  per  feet  per 
Du  mu 

second  squared.  These  aerodynamic  coefficients  were  taken 

from  Reference  4  and  are  also  available  from  Reference  1. 

In  order  to  make  the  units  consistent  in  the  body  to 

stability  axis  conversion  equations,  CT  ,  C„  ,  and  C 

L  u  m 

u  u  u 

should  first  be  multiplied  by  the  trim  velocity  before  sub¬ 
stitution  into  the  conversion  equations.  Equations  (B-57)  , 
(B-61) ,  and  { B— 6 6 )  are  particularly  affected  by  the  cor¬ 
rection;  however,  the  primed  dimensional  derivatives  X',  Z', 

u  u 

and  NT,  calculated  with  Equations  ( B— 2 ) ,  (B-6),  and  (B-10), 

remain  small  in  magnitude.  The  designs  and  simulations  in 
this  study  are  therefore  essentially  unaffected  by  the 
correction  and  are  still  valid. 


(B-81) 

( B—  8  2 ) 


The  same  error  exists  in  several  references  and 


Appendix  C:  Aerodynamic  Data 


This  appendix  contains  all  of  the  aerodynamic  data  used 
in  the  thesis.  Reference  4  is  the  source  of  most  of  the 
data;  however,  the  primed  dimensional  derivatives  were 
checked  with  the  equations  in  Appendix  B.  The  thrust  data 
is  given  in  Table  2-1,  and  the  derivative  M£c  is  calculated 
with  Equation  (B-29) . 

The  Tables  in  this  appendix  list  aircraft  aerodynamic 
coefficients  in  the  stability  axis,  primed  dimensional  deri¬ 
vatives  in  the  body  axis,  and  MULTI  data  files  for  each 
flight  condition  considered  in  this  study.  The  aerodynamic 
coefficients  are  entered  into  the  equations  of  Appendix  B  to 
calculate  the  primed  dimensional  derivatives.  The  dimen¬ 
sional  derivatives  are  entered  directly  into  the  aircraft 
models,  given  by  Equations  (2-40)  and  (2-41).  The  units  for 
the  derivatives  are:  radians  for  angles;  radians  per  second 
for  angular  rates;  feet  per  second  for  velocity;  and  g's  for 
accelerations.  MULTI  data  files  are  given  for  easy  verifi¬ 
cation  of  design  responses. 

Table  C-l  lists  aircraft  parameters  that  are  common  to 


all  flight  conditions 


Table  C-2 


Aerodynamic  Coefficients  in  the  Stability  Axis 
for  0.9  Mach  at  20,000  Feet 
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(dynamic  pressure)  =  552.11295  Ibs/ft 
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C„  ( 
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P 


)  = 


C  ( — 4^ — 
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)  = 
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0.023766 

0.000160 

0.002636 


0.000173 

0.000351 

-0.022052 

0.057563 

0.557635 

-0.000088 

0.001708 


Table  C-2  (continued) 
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Table  C-3 


Primed  Dimensional  Derivatives  in  the  Body 
for  0.9  Mach  at  20,000  Feet 


X;  =  -32.1830 

0 

Z'  =  -0.001120 
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Me  = 
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C-5 


Axis 


0.000309 

0.000130 

4.27171 

0.777221 

24.0581 

6.47269 

0.9864 

7.2370 

0.023184 

0.362530 

1.25006 

5.13710 

5.80890 

5.89254 


Table  C-4 


Healthy  Aircraft  Data  File  for  use  with  MULTI  Option  #9 
for  0.9  Mach  at  20,000  Feet 


l 

14  6  6 

oooioooooooooo 

-32.1830  -.012075  38.2906  -30.1376  0000  2.00593  2.31681  000  19.254 
-.00112  -.000022  -1.48446  .994789  0000  -.149227  -.244924  0000 
.000309  -.000130  4.27171  -.777221  0000  -24.0581  -6.47269  0  0  -.9864  0 
00000010000000 

0000  .034486  -.343554  .032636  -.997556  0  0  .0052813  .037032  .026734  0 
00000  -55.2526  -2.80004  .145674  0  0  -63.73245  10.3955  5.53185  0 
00000  7.237  -.023184  -.36253  0  0  -2.534335  -5.8089  5.89254  0 
00000000  -20  00000 
000000000  -20  0000 
0000000000  -20  000 
00000000000  -20  0  0 
000000000000  -20  0 
0000000000000-1 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
20  0  0  0  0  0 
0  20  0  0  0  0 
0  0  20  0  0  0 
0  0  0  20  0  0 
0  0  0  0  20  0 
0  0  0  0  0  1 

01000000000000 

.0001339  .0005813  44.8737  -.1857  0000  -6.09774  4.2943  0  0  -.4273  0 
OOOIOOOOOOOOOO 

00000  -6.8217  -.005034  -.08623  0  0  -.9449  -1.4433  3.3276  0 
00000010000000 
00000001000000 
N 
Y 

111111 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-1000 

1000 

1 

N 


Table  C-5 


Failed  Aircraft  Data  File  for  use  with  MULTI  Option  #9 
for  0.9  Mach  at  20,000  Feet 


l 

14  6  6 

00010000000000 

-32.1830  -.012075  38.2906  -30.1376  0000  1.002965  1.1584  1.1584  0  0  19.254 
-.00112  -.000022  -1.48446  .994789  0000  -.07461  -.1225  -.1225  000 
.000309  -.000130  4.2717  -.777221  0000  -12.0291  -3.2363  -3.2363  0  -.9864  0 
00000010000000 

0000  .03449  -.3436  .03264  -.9976  -.003326  -.000686  .000686  .0370  .0267  0 
00000  -55.2526  -2.80004  .145674  6.3411  -25.5251  25.5251  10.3955  5.5319  0 
00000  7.237  -.023184  -.36253  .642138  -.62503  .62503  -5.8089  5.89254  0 
00000000  -20  00000 
000000000  -20  0000 
0000000000  -20  000 
00000000000  -20  00 
000000000000  -20  0 
0000000000000-1 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
20  0  0  0  0  0 
0  20  0  0  0  0 
0  0  20  0  0  0 
0  0  0  20  0  0 
0  0  0  0  20  0 
0  0  0  0  0  1 

01000000000000 

.0001339  .0005813  44.8737  -.1857  0000  -3.04887  2.14715  2.14715  0  -.4273  0 
00010000000000 

00000  -6.8217  -.005034  -.08623  .1818  -.2906  .2906  -1.4433  3.3276  0 

00000010000000 

00000001000000 

N 

Y 

111111 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-1000 

1000 

1 

N 


Table  C-6 


Aerodynamic  Coefficients  in  the  Stability  Axis 
for  0.6  Mach  at  30,000  Feet 


2 

(dynamic  pressure)  =  158.81  lbs/ft 

(trim  velocity)  =  596.91  ft/sec 

(trim  angle  of  attack)  =  0.08212  radians 


=  0.439013 


C 


D 


ft/ sec 
1  . 


deg 

1 


radian 

1 

radian 


'deg' 


0.000158 

0.073559 

-1.029896 

2.437286 

0.009473 


'D  ft/sec 
u 


)  = 


CDa(deg) 


C[56e  (d^ 
t 


C  ( — i — ) 

LL$f t  deg 

C 

m 

Cmu  ^  f t/sec^ 

C  is^~) 
ma  deg 

C  .  ( — i-r - ) 

ma  radian 

r  ( — = - ) 

m  radian 

q 


0.015850 

0.0 

-0.000106 

0.004356 

-0.747194 

-2.859448 


C 

C 

c 

c 

c 

c 


(— ) 
D6  deg 

t 

(_ L_) 

Yg  deg' 

( — l _ 

y  radian 
P 

,  1 


y  'radian 
1  r 

( — - — 
ySDF  de9 

(  -  A—  ) 

^SDT  deg 


)  = 

)  = 

)  = 


0.044151 

0.000050 

0.008210 


0.000019 

0.001808 

-0.021995 

0.132102 

0.536904 

0.000051 

0.002055 


Table  C-6  (continued) 


Cm  =  -0.010229 

m  deg 

°  t 

C  (— L-)  =  -0.000383 

m5f t  deg 


0.003021 


c 

m5c  deg 

= 

-0.000463 

Cy$c  <deg) 

= 

0.001047 

C^B(deg) 

= 

-0.002209 

%'deg’ 

= 

0.001972 

C  (  1  ) 

'radian' 
P 

= 

-0.243246 

C  ( - - - ) 

n  radian 

P 

= 

-0.013813 

C  (  1  ) 

l  'radian' 

= 

0.071941 

C  ( - i-r - ) 

n  radian 
r 

= 

-0.484330 

cW'd^ 

= 

-0.002141 

CnSDF(deg’ 

= 

0.000035 

Ct5DT,d^’ 

= 

-0.001742 

C  ( 

n6DT  deg 

= 

-0.000940 

C*6r  W 

= 

0.000364 

Sr  <deg) 

= 

-0.001520 

Cs-5C  <deg) 

= 

0.000138 

Sc  (deg) 

= 

0.001121 

Table  C-7 


Primed  Dimensional  Derivatives  in  the  Body  Axis 
for  0.6  Mach  at  30,000  Feet 


=  -32.0915 

=  -0.004425 

M0  =  0.000313 

X'  =  -0.005142 

Z'  =  -0.000109 

M'  =  -0.000337 

u 

u 

u 

X'  =  23.0402 

Z'  =  -0.526422 

M'  =  2.52708 

a 

a 

a 

X'  =  -48.8785 

Z'  =  0.997184 

M'  =  -0.341902 

q 

q 

q 

X£  =  3.17035 

Z  1  =  -0.066156 

6  0 , 

M'  =  -5.86214 

et 

t 

t 

XI,  =  -2.09855 

5ft 

Z'  =  -0.111711 

5ft 

M'  =  -0.211773 
6ft 

X^T  =  10.555 

— 

m6c  =  -°*2657 

Y!  =  0.05376 

<P 

— 

— 

Yg  =  -0.154099 

Lg  =  -19.2246 

Ng  =  2.29583 

Y'  =  0.082387 

L'  =  -0.893601 

N'  =  -0.000888 

P 

P 

P 

Y^  =  -0.998322 

L^.  =  0.318845 

N^.  =  -0.278676 

Y6DF=  °*000357 

L-df=  -17.4468 

N6DF=  -°-268403 

Y’dt=  0.014398 

L5DT=  -13-5832 

N6DT=  -1-50547 

YI  =  0.021165 

0  r 

Li  =  3.92325 

o  r 

n:  =  -1.96651 

5  r 

Y1  =  0.007335 

oc 

Ll  =  0.414519 

6c 

N j  c  =  1.51008 

L0 


Table  C-8 


Healthy  Aircraft  Data  File  for  use  with  MULTI  Option  #9 


for  0.6  Mach  at  30,000  Feet 


l 

14  6  6 

oooioooooooooo 

-32.0915  -.005142  23.0402  -48.8785  0000  3.17035  -2.09855  0  0  10.555 
-.004425  -.000109  -.526422  .0997184  0000  -.066156  -.111711  0000 
.000313  -.000337  2.52708  -.341902  0000  -5.86214  -.211773  0  0  -.2657  0 
00000010000000 

0000  .05376  -.154099  .082387  -.998322  0  0  .0039565  .021165  .007335  0 
00000  -19.2246  -.893601  .318845  0  0  -20.8426  3.92325  .414519  0 
00000  2.29583  -.000888  -.278676  0  0  -.6447705  -1.96651  1.51008  0 
00000000  -20  00000 
000000000  -20  0000 
0000000000  -20  000 
00000000000  -20  0  0 
000000000000  -20  0 
0000000000000-1 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
20  0  0  0  0  0 
0  20  0  0  0  0 
0  0  20  0  0  0 
0  0  0  20  0  0 
0  0  0  0  20  0 
0  0  0  0  0  1 

01000000000000 

.0001356  .001875  10.8534  -.09592  0000  -1.3133  1.9791  0  0  -.1151  0 
OOOIOOOOOOOOOO 

00000  -1.862  .004608  -.08962  0  0  -.2060  -.4596  .7902  0 


1 

14  6 

;  e 

0  0 

0 

10  0  0 

0  0  0  0  0  0 

0 

-32. 

092  -.00514  23.040  -48 

.879  0000  1.58518  -1.04928 

-1.04928  0  0 

10.555 

-.00443  -.000109  -.52642  . 

099718  0000  -.03308  -.05586 

-.05586  0  0 

0 

.000313  -.000337  2.5271  -. 

34190  0000  -2.9311  -.1059  - 

.1059  0  -.2657  0 

0  0 

0 

0  0  0  1 

0  0  0  0  0  0 

0 

0  0 

0 

0  .0538 

-.1541  .0824 

-.9983  -.00178  .0001785  -.0001785  .02117  . 

00734  0 

0  0 

0 

0  0  -19 

.225  -.8936  . 

3188  1.6979  -8.7234  8.7234  3.9233  .41452  0 

0  0 

0 

0  0  2.2958  -.000888 

-.27868  .1882  -.1342  .1342  -1. 

9665  1.51008 

0 

0  0 

0 

0  0  0  0 

0  -20  000 

0  0 

0  0 

0 

0  0  0  0 

0  0  -20  0  0 

0  0 

0  0 

0 

0  0  0  0 

0  0  0  -20  0 

0  0 

0  0 

0 

0  0  0  0 

0000  -20 

0  0 

0  0 

0 

0  0  0  0 

00000  -20  0 

0  0 

0 

0  0  0  0 

0  0  0  0  0  0 

-1 

0  0 

0 

0  0  0 

0  0 

0 

0  0  0 

0  0 

0 

0  0  0 

0  0 

0 

0  0  0 

0  0 

0 

0  0  0 

0  0 

0 

0  0  0 

0  0 

0 

0  0  0 

0  0 

0 

0  0  0 

20  0 

1  0 

10  0  0 

0  20 

0 

0  0  0 

0  0 

20 

0  0  0 

0  0 

0 

20  0  0 

0  0 

0 

0  20  0 

0  0 

0 

0  0  1 

0  1 

0 

0  0  0  0 

0  0  0  0  0  0 

0 

.0001356  .001875  10.8534  - 

.09592  0000  -.65665  .98955 

.9855  0  -.1151  0 

0  0 

0 

10  0  0 

0  0  0  0  0  0 

0 

0  0 

0 

0  0  -1.1 

362  .004608  - 

.08962  .04816  -.05483  .05483  - 

.4596  .7902  0 

0  0 

0 

0  0  0  1 

0  0  0  0  0  0 

0 

0  0 

0 

0  0  0  0 

1  0  0  0  0  0 

0 

N 

Y 

1  1 

1 

111 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

I  -1000 

I  1000 

1 

N 

q  (dynamic  pressure) 


=  1129.3126  lbs/ft 
VT  (trim  velocity)  =  1591.75  ft/sec 
aT  (trim  angle  of  attack)  =  0.02932  radians 


CT  (: 


'L  ft/sec 
u 


CT  (; 


)  = 


'Vdeg 


CT.  (- 


'L‘  radian 

a 


CL  ^radian 

q 

Se,^ 


m 


r  ( — - - 

mu lf t/sec 

Cm  (35?> 

a 


)  - 


)  = 


)  - 


C  •  (: 


'm‘  radian 

a 


C  ( — - - 

m  radian 

q 


)  = 


)  = 


0.060695 

0.000027 

0.066129 

-2.974699 

3.458724 

0.006180 

0.004004 

0.0 

-0.000046 

-0.010705 

0.501981 

-1.889366 


'D  ft/sec 


CD  <di?> 

a 


)  = 


C^5et 

cys'3i?' 


C  — — 
y  radian 
P 

1 

C  ( 


)  = 


y  radian 


)  = 


C  (-5^-) 

ySDF  deg 

r  (_ — ) 

*6dt  deg 


0.046574 

0.000002 

0.001457 


-0.000135 

0.000565 

-0.019192 

0.026463 

0.516864 

0.000362 

0.000281 


Table  C-10  (continued) 


<3ig> 


C  (  ■■  ■) 
m5  deg 

t 


C 

C 

C 

C 

C 

C 

C 

C 


(— — ) 
m5c  deg 

(_L_) 

%  deg 

( — — - ) 

l  radian 
P 

( _ i - ) 

Jtr  radian 

*6 df  <d®g) 

(_1_) 
*6  DT  deg 

(  — — ) 
£5r  deg 

(_i_) 

£6c  deg 


1 

-0.008061 

^deg1 

-0.001432 

— 

0.000050 

C  (^) 

y6c  deg 

-0.001702 

-0.257331 

0.018261 

Cng (deg* 

C  (  1 

~n  'radian 
P 

C  (  1 

^nr 'radian 

-0.000250 

Cn5DF(degl 

-0.000811 

Cn6DT<deg) 

0.000077 

Sr  ‘deg1 

0.000232 

Sc  (deg) 

0.000552 

0.001541 

0.00088 

0.004371 

-0.246367 

-0.000236 

-0.000317 

-0.000349 

0.000973 


Table  C-ll 


Primed  Dimensional  Derivatives  in  the  Body 
for  1.6  Mach  at  30,000  Feet 


X^  =  -32.1861 

Z'  =  -0.000593 

0 

II 

X'  =  -0.030046 
u 

Z'  =  -0.000002 
u 

M'  = 
u 

X'  =  45.1941 

a 

Z'  =  -1.25100 

a 

M’  = 
a 

X'  =  -46.5132 

q 

Z*  =  0.995991 

q 

M'  = 

q 

XI  =  9.40431 

5et 

Z'  =  -0.115337 

6et 

X'  =  -13.3021 

Sft 

Z'  =  -0.075084 

6ft 

■V 

X^T  =  7.230 

— 

M6c 

Y'  =  0.02022 

— 

— 

Y'  =  -0.358564 
p 

Lg  =  -100.032 

Ns  ■ 

Yp  =  0.029374 

L'  =  -2.46287 

P 

»p  ■ 

=  -0.998410 

L^.  =  0.160287 

Nr  = 

Y'df=  0.006763 

L  6DF=  -14*1689 

N '  = 

6DF 

Yj  DT=  °*005250 

L5DT=  -46*6078 

N6DT= 

Y^r  =  0.010313 

L  L  =  4.96942 

or 

N6r  = 

Y’  =  0.028790 

6C 

Li  =  12.0662 

oc 

N6c  = 

Axis 


0.000075 

0.000802 

43.8012 

0.351737 

32.8919 

5.85004 

0.2043 

7.42422 

0.005040 

0.383316 

2.37673 

3.45190 

3.27086 

9.36290 


Table  012 


Healthy  Aircraft  Data  File  for  use  with  MULTI  Option  #9 
for  1.6  Mach  at  30,000  Feet 


l 

14  6  6 

00010000000000 

-32.1861  -.030046  45.1941  -46.5132  0000  9.40431  -13.3021  000  7.230 
-.000593  -.000002  -1.251  .995991  0000  -.115337  -.075084  0000 
-.000075  .000802  -43.8012  -.351737  0000  -32.8919  -5.85004  0  0  .2043  0 
00000010000000 

0000  .02022  -.358564  .029374  -.99841  0  0  .0080755  .010313  .02879  0 
00000  -100.032  -2.46287  .0160287  0  0  -25.82085  4.96942  12.0662  0 
00000  7.42422  -.00504  -.383316  0  0  -3.2397  -3.27086  9.3629  0 
00000700  -20  00000 
000000000  -20  0000 
0000000000  -20  000 
00000000000  -20  00 
000000000000  -20  0 
0000000000000-1 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
20  0  0  0  0  0 
0  20  0  0  0  0 
0  0  20  0  0  0 
0  0  0  20  0  0 
0  0  0  0  20  0 
0  0  0  0  0  1 

01000000000000 

-.00003249  .0004463  42.865  .04579  0000  -8.5483  1.1772  0  0  .08851  0 
00010000000000 

00000  -14.5086  .0004118  -.08747  0  0  -1.004339  -.90723  5.4795  0 

00000010000000 

00000001000000 

N 

Y 

111111 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-1000 

1000 

1 


Table  C-13 


Failed  Aircraft  Data  File  for  use  with  MULTI  Option  #9 


for  1.6  Mach  at  30,000  Feet 


l 

14  6  6 

00010000000000 

-32.186  -.030046  45.1941  -46.5132  0000  4.7022  -6.6511  -6.6511  0  0  7.230 
-.000593  -.000002  -1.251  .995991  0000  -.05767  -.03754  -.03754  000 
-.000075  .000802  -43.8012  -.3517  0000  -16.4460  -2.925  -2.925  0  .2043  0 
00000010000000 

0000  .02022  -.35856  .02937  -.9984  -.000656  .00338  -.00338  .01031  .02879  0 
00000  -100.032  -2.4629  .0160287  5.826  -7.0845  7.0845  4.9694  12.0662  0 
00000  7.42422  -.00504  -.383316  .4315  -1.1884  1.1884  -3.27086  9.3629  0 
00000000  -20  00000 
000000000  -20  0000 
0000000000  -20  000 
00000000000  -20  00 
000000000000  -20  0 
0000000000000-1 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
20  0  0  0  0  0 
0  20  0  0  0  0 
0  0  20  0  0  0 
0  0  0  20  0  0 
0  0  0  0  20  0 
0  0  0  0  0  1 

01000000000000 

-.00003249  .0004463  42.865  .04579  0000  -4.2742  .5886  .5886  0  .08851  0 
00010000000000 

00000  -14.5086  .0004118  -.08747  .1545  -.3477  .3477  -.90723  5.4795  0 
00000010000000 
00000001000000 
N 
Y 

111111 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-1000 

1000 

1 

N 


The  aerodynamic  coefficients  for  the  previous  three 
flight  conditions  appear  in  Reference  4  and  were  originally 
obtained  from  Reference  1;  however,  in  both  these  refer¬ 
ences,  landing  approach  data  is  given  for  0.2  Mach  at  30 
feet  with  the  flaperons  completely  deflected  downward.  The 
trim  angle  of  attack  <*T  for  this  flight  condition  is  14.98 
degrees.  For  the  study  in  this  thesis  it  is  desirable  to 
evaluate  the  aircraft's  performance  with  the  flaperons 
trimmed  at  zero  degrees  and  at  a  higher  speed,  which  yields 
a  smaller  aT>  A  landing  approach  condition  of  0.3  Mach  at 
30  feet  is  therefore  chosen  for  this  thesis. 

It  is  assumed  that  the  aircraft's  equations  of  motion 
are  linear  for  a  small  change  in  velocity;  therefore,  most 
of  the  aerodynamic  coefficients  for  the  0.2  Mach  flight 
condition  are  applicable  to  tht  0.3  Mach  flight  condition. 
The  lift  coefficient  CL  changes  with  the  inverse  of  the 
dynamic  pressure  q,  and  the  drag  coefficient  CD  is  decreased 

by  Cn  times  the  maximum  flaperon  deflection,  because  the 
of 

new  trim  condition  for  the  flaperons  is  at  zero  degrees. 

The  trim  angle  of  attack  is  calculated  from  the  lift 
and  moment  equations; 


+  C  am 
L  T 
a 

(C-l) 

*  Cm  aT 
a 

(C-2 ) 

These  equations  are  solved  simultaneously  for  a™,  yielding 


Table  014 


Aerodynamic  Coefficients  in  the  Stability  Axis 
for  0.3  Mach  at  30  Feet 


q  (dynamic  pressure)  =  133.163  lbs/ft 

VT  (trim  velocity)  =  334.9  ft/sec 

«T  (trim  angle  of  attack)  =  0.14484  radians 


r  ( _ = _ 

'“L  'ft/sec 
u 


Vdeg1 


)  = 


CL*  ^radian 

a 


)  = 


CL  ^radian' 

q 

CL«etlai?l 


'm 


n  i — _ _ 

m  'ft/sec 
u 

c,»  'a??1 

a 


)  = 


Cm*  ^radian' 

a 


Cm  (— TT-— ) 
nig  rcKaicin 


0.526120 

-0.000178 

0.053031 

-0.775892 

1.428590 

0.010306 

0.004004 

0.0 

0.000032 

0.001755 

0.017115 

-6.716320 


'D 


Cr,  (: 


'D  ft/sec 
u 


CD  (deg} 

a 


)  = 


cv‘a*) 

CD«£t,a^' 

cy  'as?1 
B 


Cy  ^radian 
P 


)  = 


^y  ^radian 
1  r 


)  = 


Cy  <35?» 

y<5  DF  y 

C  (-r^-) 
y6  DT  de^ 


0.167204 

0.000063 

0.027636 


0.002519 

0.002519 

-0.020613 

0.214209 

0.640146 

0.001493 

0.003154 


Table  C-14  (continued) 


CmSet(deg) 
%.  <ai5> 

t 

C  (^:) 
m5c  deg 


CfLg  *deg* 


C0  ( 


l  radian 
P 


C„  (: 


l  radian 
r 

C  ( 
lSDT  de9 

Sr 


'he  deg 


-0.010622 

0.001000 

-.0195035 

-0.003015 

-0.201582 

0.224680 

-0.001373 

-0.001658 

0.000033 

0.000074 


<3ig> 


0.003103 


C 

C 

c 

c 

c 

c 

c 

c 


(— L_ ) 

y6c  deg 

ng  (deg) 

( — — - ) 

n  radian 

P 

( — - - ) 

n  radian 


n6DFldeg 
(  -  ^ 

n6DT  deg 

(_i_ 

n6r  deg 

(JL 

n5c  deg 


) 

) 

) 

) 


0.000893 

0.002824 

-0.065499 

-0.509547 

-0.000566 

-0.000769 

-0.001716 

0.000665 


Table  C-15 


Primed  Dimensional  Derivatives  in  the  Body  Axis 
for  0.3  Mach  at  30  Feet 


x0  = 

-31.8627 

Z '  = 

0 

-0.013880 

M' 

0 

-0.000033 

Xu  = 

-0.016994 

Z  1 
u 

-0.000326 

M' 

u 

-0.000364 

X'  = 
a 

-36.9184 

Z  ' 
a 

-0.629959 

M' 

a 

0.833706 

X' 

q 

-48.3011 

Z  ' 

q 

0.995634 

M' 

q 

-0.950374 

x:  = 

5et 

-3.52376 

-0.110589 

-5.10864 

X 1  = 
5ft 

-8.00387 

-0.090118 

0.48071 

X5T 

22.291 

— 

M'  = 
5C 

-0.16371 

% 

0.095141 

— 

— 

►< 

V*  - 

il 

-0.215833 

H  ' 

-23.1367 

N'  = 

8 

2.53712 

Y'  = 
P 

0.141584 

L'  = 
P 

-1.23522 

N' 

P 

-0.027433 

Y'  = 
r 

-0.994562 

L’  = 
r 

1.40799 

N'  = 
r 

-0.414331 

Y '  = 
5DF 

0.015633 

L5DF  = 

-8.76403 

N '  = 
6DF 

-0.890051 

Y5DT  = 

0.033025 

L5DT  = 

-10.5014 

N '  = 

<$DT 

-1.16925 

Y'  = 

5  r 

0.032491 

4r  ' 

1.86570 

N'  = 
5r 

-1.88884 

y;  = 

5c 

0.009350 

L5c 

-0.13479 

= 

5c 

0.748726 

Table  016 


Healthy  Aircraft  Data  File  for  use  with  MULTI  Option  #9 
for  0.3  Mach  at  30  Feet 


l 

14  6  6 

00010000000000 

-32.8627  -.01699  -36.9184  -48.3011  0000  -3.5238  -8.0039  000  22.291 
-.01388  -.000326  -.62996  .995634  0000  -.110589  -.090118  0000 
-.0000333  -.000364  .83371  -.950374  0000  -5.10864  .480707  0  0  -.1637  0 
00000010000000 

0000  .095141  -.215833  .145841  -.994562  0  0  .023889  .032491  .0093504  0 
00000  -23.1367  -1.23522  1.40799  0  0  -11.38938  1.86570  -.13479  0 
00000  2.537  -.027433  -.41433  0  0  -1.182364  -1.8888  .748726  0 
00000000  -20  00000 
000000000  -20  0000 
0000000000  -20  000 
00000000000  -20  0  0 
000000000000  -20  0 
0000000000000-1 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
20  0  0  0  0  0 
0  20  0  0  0  0 
0  0  20  0  0  0 
0  0  0  20  0  0 
0  0  0  0  20  0 
0  0  0  0  0  1 

01000000000000 

-.000014  .0032362  6.91315  -.3663  0000  -1.06302  1.4553  0  0  -.0709  0 
00010000000000 

00000  -1.1456  -.001706  -.12294  0  0  -.2638  -.48038  .42162  0 
00000010000000 
00000001000000 
N 
Y 

111111 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-1000 

1000 

1 


1 

14  6  6 

00010000000000 

-32.8627  -.01699  -36.9184  -48.301  0000  -1.7619  -4.0019  -4.0019  0  0  22.291 
-.01388  -.000326  -.62996  .99563  0000  -.05529  -.045059  -.045059  000 
-.0000333  -.000364  .83371  -.950374  0000  -2.5543  .24035  .24035  0  -.16371  0 
00000010000000 

0000  .09514  -.2158  .14584  -.995  -.004128  .007816  -.007816  .03249  .00935  0 
00000  -23.1367  -1.23522  1.40799  1.3127  -4.382  4.382  1.866  -.13479  0 
00000  2.537  -.027433  -.41433  .14616  -.44503  .44503  -1.8888  .748726  0 
00000000  -20  00000 
000000000  -20  0000 
0000000000  -20  000 
00000000000  -20  00 
000000000000  -20  0 
0000000000000-1 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
20  0  0  0  0  0 
0  20  0  0  0  0 
0  0  20  0  0  0 
0  0  0  20  0  0 
0  0  0  0  20  0 
0  0  0  0  0  1 

01000000000000 

-.000014  .0032362  6.91315  -.3663  0000  -.53151  .57277  .57277  0  -.070922  0 
00010000000000 

00000  -1.1456  -.001706  -.12294  .020385  -.111503  .111503  -.48038  .42162  0 
00000010000000 
00000001000000 
N 
Y 

111111 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-72 

72 

1 

-1000 

1000 

1 

N 
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Appendix  D:  Continued  Design  Responses 

Chapter  II  presents  the  simulation  response  plots  for 
the  transonic  flight  condition  (0.9  Mach  at  20,000  feet). 
This  appendix  gives  the  response  plots  for  the  other  three 
flight  conditions  considered  in  this  thesis,  which  are 
subsonic  flight  (0.6  Mach  at  30,000  feet),  supersonic  flight 
(1.6  Mach  at  30,000  feet),  and  landing  approach  (0.3  Mach  at 
30  feet) .  The  order  of  presentation  is  similar  to  that  of 
Chapter  II. 

Subsonic  Flight  (0.6  Mach  at  30,000  Feet) 

The  subsonic  flight  condition  is  chosen  to  demonstrate 
maneuverabiltiy  at  a  cruising  speed  for  the  aircraft.  The 
dynamic  pressure  is  relatively  low  because  of  the  high 
altitude;  therefore,  the  control  surfaces  are  less  effective 
at  this  flight  condition  than  at  the  transonic  flight 
condition. 

G-Command  Maneuver.  Tables  D-l  and  D-2  list  the  design 
parameters  for  the  g-command  maneuver  for  the  healthy  and 
failed  right  elevator  models.  As  shown  in  the  tables,  the 
same  design  parameters,  a,  and  e  ,  are  used  for  both 
aircraft  models. 

Figures  D-l  through  D-7  show  the  simulation  responses 


for  the  g-command  maneuver.  For  this  maneuver  the  healthy 
aircraft  longitudinal  responses  change  very  little  when  the 
right  horizontal  tail  fails.  Also,  the  healthy  aircraft 


lateral  responses  are  zero  and  are  very  small  with  the 
failed  model.  For  these  reasons  Figures  D-l  through  D-3 
give  simulation  responses  for  both  the  healthy  and  failed 
aircraft  models. 

The  maximum  normal  acceleration  command  of  3.1g  is 
limited  by  the  flaperon  deflection.  For  the  failure  case, 
the  left  elevator  deflects  approximately  twice  the  amount 
required  when  both  tail  halves  are  working.  The  flaperon 
deflection  is  slightly  asymmetric  for  the  failure  case  in 
order  to  compensate  for  the  rolling  moment  created  by  the 
left  elevator  deflection.  The  rudder  and  canard  compensate 
for  the  adverse  yaw  caused  by  a  difference  in  drag  between 
the  left  and  right  halves  of  the  aircraft. 


Table  D-l 


G-Command:  Healthy  Model,  0.6  Mach  at  30,000  Feet 
Sampling  Time:  T  =  0.02  second 


a  =  1.2 


e  =  1.0 


I  =  diag  {  2.0,  0.04,  1.0,  1.0,  1.0,  1.0} 


0.0 

-0.3473E-14 

0.0 

0.0 

0.0 

0.1895E*01 


-0.3565E-04 

0.9869E-03 

0.0 

0.0 

0.0 

0.4139E-02 


-0.8330E-01 

-0.5527E-01 

0.0 

0.0 

0.0 

0.2806E+00 


-0.5402E-02 
0.3198E-02 
0 . 1806E-01 
0 . 8364E-01 
0.1166E+00 
0.4517E-01 


0 . 5  400E-04 
-0.3198E-04 
-0.9219E-02 
0.2127E-02 
-0.1166E-02 
-0.4516E-03 


0.515BE-02 

-0.3054E-02 

-0.3528E-01 

-0.1756E*00 

-0.1114E*00 

-0.4313E-01 


*1 


0.0 

-0 . 4278E-04 

-0.9996E-01 

-0 . 6482E-02 

0 . 6480E-04 

0.6189E-02 

*0 . 4167E-14 

0.X184E-02 

-0.6633E-01 

0.3838E-Q2 

-0.3B37E-04 

-0 . 3665E-02 

0.0 

0.0 

0.0 

0.2168E-01 

-0.U06E-01 

-0.4233E-01 

0.0 

0.0 

0.0 

0.1004E*00 

0.2553E-02 

-0.2108E*00 

0.0 

0.0 

0.0 

0 . 1400E+00 

-0.1399E-02 

-0 . 1336E+00 

0.2274E*01 

0.4966E-02 

0.3367E*00 

0.5420E-01 

-0.5419E-03 

-0.5175E-01 

Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =  0.0 

An  =  3.1  g  (step) 

P 

q  =  0.1672  radian/second  (step) 

A  =  0.0 

P 

p  =  0.0 

r  =0.0 


D-3 


Table  D-2 


G-Command:  Failed  Model,  0.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


o  =  1.2 


e  =  1.0 


Z  =  diag  {2.0,  0.04,  1.0,  1.0,  1.0,  1.0} 


0.4482E-14  -0.7129E-04  -0.1665E*00  -0.1081E-01  0.1129E-03  0.1033E-01 

-0.2244E-14  0.9786E-03  -0.7458E-01  0.2420E-01  -0.1137E-01  -0.4531E-01 

-0.3025E-14  0.9952E-03  -0.3599E-01  -0.1779E-01  0.1131E-01  0.3920E-01 

-0.7669E-16  -0.5852E-05  -0.1367E-01  0.B571E-01  0.6232E-03  -0.1B06E*00 

-0.4928E-15  -0.2051E-06  -0.4791E-03  0.1167E*00  -0.1218E-02  -0.1115E*00 

0.1895E.01  0.4138E-02  0.2804E»00  0.4520E-01  -0.4718E-03  -0.4319E-01 


0.5379E-14  -0.8S54E-04  -0.1999E-00  -0 . 1297E-01  0.1354E-03  0.1240E-01 

-0.2693E-14  0.U74E-02  -0.8950E-01  0.2903E-01  -0.1365E-0!  -0.5437E-01 

-0. 3630E-14  0.1194E-02  -0.4319E-0X  -0.2135E-01  0.1357E-01  0.4703E-01 

-0.9202E-16  -0.7022E-05  -0.1641E-0X  0.1029E»00  0.7479E-03  -0.2167E»00 

-0.5913E-15  -0.2461E-06  -0.5749E-03  0.1400E»00  -0.1462E-02  -0.133BE*00 

0.2274E*01  0.4966E-02  0.3365E*00  0.5424E-01  -0.5662E-03  -0.51B3E-01 


Input  Ramp  Time:  0.4  second 
Command  Vector: 


u  =0.0 


=  3.1  g  (step) 


q  =  0.1672  radian/second  (step) 


l§] 


T  'tIE .  ‘SECONDS 


Figure  D-2:  0.6M  Healthy/Failed  G-Command  —  Angle  of 

Attack,  Roll  Angle,  Sideslip  Angle,  and  Yaw  Rate 


Figure  D-6:  0.6M  Healthy  G-Command  —  Aileron,  Rudder,  and 

Canard  Deflections  and  Thrust 


Figure  D-7:  0.6M  Failed  G-Command  — 

Deflections  and  Thrust 


Rudder  and  Canard 


Pitch  Pointing  Maneuver.  Tables  D-3  and  D-4  give  the 
design  parameters  for  the  pitch  pointing  maneuver.  As  with 
the  g-command  maneuver  the  design  parameters  a  ,  Z_,  and  e  are 
the  same  for  the  healthy  and  failed  models. 

The  simulation  responses  for  the  pitch  pointing  maneu¬ 
ver  are  given  in  Figures  D-8  through  D-14.  The  maximum 
pointing  angle  commanded  is  limited  by  the  flaperon  deflect¬ 
ion.  For  the  same  reasons  discussed  in  the  g-command 
maneuver  section,  Figures  D-8  through  D-10  give  simulation 
responses  for  both  the  healthy  and  failed  aircraft  models. 
For  the  failure  case  the  left  elevator  must  deflect  approx¬ 
imately  twice  the  angle  of  the  healthy  elevator  deflection. 
Also,  for  the  failure  model,  the  flaperons  deflect  asymme¬ 
trically  to  counter  the  roll  induced  by  the  failed  right 
horizontal  tail,  and  the  rudder  and  canard  deflect  to 
compensate  for  adverse  yaw. 


Table  D-3 


Pitch  Pointing:  Healthy  Model,  0.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  2.0 
e  =  1.0 

Z  =  diag  (  2.0,  1.5,  1.0,  1.0,  1.0,  1.0> 


0.0 

-0.1337E-02 

-0 . 8330E-01 

-0.5402E-02 

0.5400E-04 

0.5158E-02 

-0.3473E-14 

0.3701E-01 

-0 . 5527E-01 

0.3198E-02 

-0.3198E-04 

-0.3054E-02 

0.0 

0.0 

0.0 

0.1806E-01 

-0 . 9219E-02 

-0.3528E-01 

0.0 

0.0 

0.0 

O.B364E-01 

0.2127E-02 

-0.1756E*00 

0.0 

0.0 

0.0 

0.1166E+00 

-0.1166E-02 

-0 . 1114E+00 

0. 1895E+01 

0. 1552E+00 

0. 2806E+00 

0.4517E-01 

-0.4516E-03 

-0 . 4313E-01 

0.0 

-0.2674E-02 

-0.1666E+00 

-0.1080E-01 

0.1080E-03 

0.1032E-01 

-0.6946E-14 

0 . 7402E-01 

-0 . 1105E+00 

0.6397E-02 

-0.6395E-04 

-0.6108E-02 

0.0 

0.0 

0.0 

0.3613E-01 

-0.1844E-01 

-0.7055E-01 

0.0 

o 

o 

o 

o 

0.1673E*00 

0.4254E-02 

-0 . 3513E*00 

0.0 

0.0 

0.0 

0.2333E+00 

-0.2332E-02 

-0 . 2227E+00 

0 . 3790E+01 

0 . 3104E+00 

0.5612E*00 

0.9034E-01 

-0.9031E-03 

-0 . 8625E-01 

Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =0.0 

An  =  0.0 

P 

q  =  0.05934  radian/second  (1  second  pulse) 


Table  D-4 


Pitch  Pointing:  Failed  Model,  0.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  2.0 

e  =  1.0 

Z  =  diag  {  2.0,  1.5,  1.0,  1.0,  1.0,  1.0} 


0 .  4482E-14 

-0 . 2673E-Q2 

-0 . 1665E+00 

-0 . 1081E-01 

0.1129E-03 

0 . 1 03  3E-01 

-0 . 2244E-14 

0 . 3670E-01 

-0.7458E-01 

0. 2420E-01 

-0.1137E-01 

-0.4531E-01 

-0. 3025E-14 

0.3732E-01 

-0 . 3599E-01 

-0 . 1779E-01 

0.1131E-01 

0 . 3920E-01 

-0 . 7669E-16 

-0.2194E-03 

-0.1367E-01 

0 . 8571E-01 

0.6232E-03 

-0 . 1806E+00 

-0.4928E-15 

-0.7690E-05 

-0. 4791E-03 

0 . 1167E+00 

-0.1218E-C2 

-0 . 11 15E+00 

0 . 1 895E+01 

0. 1552E+00 

0. 2804E+00 

0 . 4520E-01 

-0.4718E-03 

-0.4319E-01 

0.8965E-14 

-0.5346E-02 

-0.3331E-00 

-0.2162E-01 

0 . 2257E-03 

0.2066E-01 

-0 . 4488E-14 

0 . 7340E-01 

-0.1492E*00 

0.4B39E-01 

-0.22 75E-01 

-0.9062E-01 

-0.6051E-14 

0.7464E-01 

-0.7199E-01 

-0.3559E-01 

0.2262E-01 

0 . 7839E-01 

-0.1534E-15 

-0.4389E-03 

-0.2734E-01 

0 . 1714E+00 

0.1246E-02 

-0 . 3612E-00 

-0.9855E-15 

-0.1538E-04 

-0.9582E-03 

0. 2334E+00 

-0.2437E-02 

-0.2231E-00 

0. 3790E+01 

0.3104E-00 

0.5608E+00 

0.9039E-01 

-0.9436E-03 

-0.8639E-01 

Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =  0.05934  radian/ second  (1  second  pulse) 


Figure  D-10:  0.6M  Healthy/Failed  Pitch  Pointing 

Acceleration  and  Lateral  Acceleration 


Figure  D-ll:  0.6M  Healthy  Pitch  Pointing  —  Elevator  and 
Flap  Deflections 


Figure  D-12:  0.6M  Failed  Pitch  Pointing  —  Left  Elevator 

and  Flaperon  Deflections 
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0-13 :  °*6M  Healthy  Pitch  Pointing  —  Aileron, 

Rudder,  and  Canard  Deflections  and  Thrust 


Tables  D-5  and  D-6 


Longitudinal  Translation  Maneuver, 
give  the  design  parameters  for  the  longitudinal  translation 
maneuver.  Once  again  the  design  parameters  a,  z,  and  e  are 
the  same  for  the  healthy  and  failed  models. 

The  simulation  responses  for  the  longitudinal  transla¬ 
tion  maneuver  are  given  in  Figures  D-15  through  D-21.  As 
discussed  in  the  g-command  maneuver  section  Figures  D-15 
through  D-17  give  simulation  responses  for  both  the  healthy 
and  failed  aircraft  models.  The  maximum  longitudinal 
translation  command  of  0.5g  for  one  second,  which  yields  a 
normal  velocity  of  approximately  22.5  feet/second,  is 
limited  by  the  flaperon  deflection. 

As  is  usually  the  case  with  longitudinal  maneuvers,  the 
failure  model  requires  that  the  left  elevator  must  deflect 
approximately  twice  the  angle  of  the  healthy  elevator  de¬ 
flection.  Also,  for  the  failure  model,  the  flaperons 
deflect  asymmetrically  to  counter  the  roll  induced  by  the 
failed  right  horizontal  tail,  and  the  rudder  and  canard 
deflect  to  compensate  for  adverse  yaw. 


D-16 


Table  D-5 


Longitudinal  Translation:  Healthy  Model 
0.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


o  =  3.5 

e  =  1.0 


£  =  diag  (2.0,  1.0,  0.8,  1.0,  1.0,  1 . 0  ) 


«0 


0.0 

-0 . 8913E-03 

-0.6664E-01 

-0.5402E-02 

0 . 5400E-04 

0 

0.3473E-15 

0.2467E-01 

-0 . 4422E-01 

0.3198E-02 

-0 . 3198E-04 

-0 

0.0 

0.0 

0.0 

0 . 1 806E-01 

-0 . 9219E-02 

-0 

0.0 

0.0 

0.0 

0 . 8364E-01 

0.2127E-02 

-0 

0.0 

0.0 

0.0 

0 . 1166E*00 

-0 . 1166E-02 

-0 

0 . 1895E+00 

0. 1035E+00 

0.2245E+00 

0 . 4517E-01 

-0 . 45 16E-0  3 

-0 

*1 


0.0 

-0 . 3120E-02 

-0 . 2332E*00 

-0.1B91E-01 

0.1890E-03 

0.1215E-14 

0.8635E-01 

-0.1548E*00 

0 . 1 119E-01 

-0.U19E-O3 

0.0 

0.0 

0.0 

0 . 6322E-01 

-0.3226E-01 

0.0 

0.0 

0.0 

0 . 2927E+00 

0.7445E-02 

0.0 

0.0 

0.0 

0 . 4082E*00 

-0.40B1E-02 

0.6632E*00 

0 . 362XE*00 

0.7B57E»0Q 

0 . 158 1E+00 

-0.1580E-02 

Input  Ramp  Time:  0.2  second 
Command  Vector: 


u  =0.0 


An  =  0.5  g  (1  second  pulse) 
P 

q  =0.0 


p  =0.0 

r  =0.0 
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5158E-0 2 
3054E-Q2 
3528E-01 
1756E-00 
11 1 4E  +  00 
4313E-01 


1805E-01 
1069E-01 
1235E*00 
6148E-00 
3  898E  +  00 
1509E+00 


Table  D-6 


Longitudinal  Translation:  Failed  Model, 
0.6  Mach  at  30,000  Feet 


Sampling  Time: 

T  =  0. 

02  second 

a 

=  3 

i.5 

e 

=  1 

.0 

1  M 

=  diag  (2.0, 

• 

o 

* 

o 

• 

H 

03 

• 

O 

• 

o 

• 

O 

__ 

0.4482E-15 

-0.1782E-02 

-0 . 1332E+00 

-0 . 1081E-01 

0. 1129E-03 

0.1033E-01 

-0 . 2244E-15 

0 . 2447E-01 

-0.5967E-01 

0 . 2420E-01 

-0.1137E-01 

-0.4531E-01 

-0.3025E-15 

0.2488E-01 

-0.2879E-01 

-0 . 1779E-0 1 

0.1131E-01 

0 . 3920E-01 

=  ! 

— 0 

-0.7669E-17 

-0.1463E-03 

-0.1094E-01 

0 . 8571E-01 

0.6232E-03 

-0 . 1806E+00 

-0.4928E-16 

-0.5127E-05 

-0.3833E-03 

0 . 1 167E+00 

-0.1218E-02 

-0.1115E*00 

0.1895E*00 

0.1035E*00 

0.2243E*00 

0 . 4520E-01 

-0.4718E-03 

-0.4319E-01 

0.1569E-14 

-0.6238E-02 

-0 . 4663E*00 

-0 . 37B4E-01 

0 . 3950E-03 

0 . 3616E-01 

-0 . 7854E-15 

0.8563E-01 

-0 . 2088E+Q0 

0 . 8468E-01 

-0.3981E-01 

-O.1586E*0O 

-0. 1059E-14 

0.8708E-01 

-0 . 1008E+00 

-0 . 6228E-01 

0.3958E-01 

0.1372E*00 

K. 

s 

-1 

-0.2684E-16 

-0.5120E-03 

-0.3828E-01 

0 . 3000E+00 

0.2181E-02 

-0 . 6321E*00 

-0.1725E-15 

-0.1794E-04 

-0.1342E-02 

0 . 4085E+00 

-0.4264E-02 

-0 . 3904E*00 

0.6632E+00 

0 . 3621E*00 

0 . 7851E+QO 

0 . 1582E+00 

-0 . 165XE-02 

-0 . 1512E+00 

Input  Ramp  Time:  0.2  second 
Command  Vector: 
u  =0.0 

A_  =  0.5  g  (1  second  pulse) 


Lflr.  stfires 


Figure  D-15:  0.6M  Healthy/Failed  Longitudinal  Translation 

—  Pitch  Angle,  Angle  of  Attack,  and  Forward  Velocity 


Figure  D-16:  0.6M  Healthy/Failed  Longitudinal  Translation 

—  Flight  Path  Angle,  Roll  Angle,  Sideslip  Angle,  and  Yaw 


0EG/5EC ) 


Figure  D-17:  0.6M  Healthy/Failed  Longitudinal  Translation 

—  Pitch  Rate,  Normal  Acceleration,  and  Lateral 
Acceleration 
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Roll  About  Velocity  Vector  Maneuver.  As  discussed  in 
Chapter  IV  the  roll  maneuver  requires  some  modification  to 
the  healthy  and  failure  models.  Tables  D-7  and  D-8,  give 
the  design  parameters  a,  j,  and  e  for  the  roll  maneuver. 

The  parameters  are  the  same  for  the  healthy  and  failed 
models . 

The  simulation  responses  for  the  roll  maneuver  are 
given  m  Figures  D-22  through  D-28.  For  this  maneuver  the 
longitudinal  responses  are  very  small  for  both  the  healthy 
and  failure  models.  Also,  the  healthy  aircraft  lateral 
responses  change  very  little  when  the  right  horizontal  tail 
fails.  For  these  reasons  Figures  D-22  through  D-24  give 
simulation  responses  for  both  the  healthy  and  failed 
aircraft  models. 

The  roll  rate  commanded  is  limited  to  200  degrees  per 
second  for  three  seconds  in  order  to  avoid  generating  too 
larce  a  sideslip  angle. 

The  aileron  deflection  for  the  healthy  aircraft  is  les 
than  the  flaperon  deflections  for  the  failure  case  because 
the  aileron  deflection  includes  a  contribution  from  the 
horizontal  tail.  When  the  right  horizontal  tail  fails  the 
f laperons  must  deflect  more  to  achieve  the  commanded  roll 
rate.  The  rudder  and  canard  deflections  are  unaffected  by 
the  failure  of  the  right  horizontal  tail. 


Table  D-7 


Roll  About  Velocity  Vector:  Healthy  Model, 
0.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  2.0 
e  =  1.0 

1  =  diag  (  2.0,  1.0,  0.8,  0.25,  0.18,  0.3  > 


0.0 

-0 . 8913E-03 

-0.6664E-01 

-0 . 1 350E-02 

0 . 9721E-0S 

0. 1547E-02 

-0.3473E-14 

0 . 2467E-01 

-0 . 4422E-01 

0 . 7996E-0  3 

-0.5756E-05 

-0 . 9161E-03 

0.0 

0.0 

0.0 

0.4516E-02 

-0 . 1659E-02 

-*0 . 1058E-01 

o 

o 

0.0 

0.0 

0.209 1E-0 1 

0 . 3829E-0  3 

-0 . 5269E-0 1 

o 

o 

o 

o 

o 

o 

0 . 29 16E-01 

-0 .  2099E-03 

-0 . 3  341E-01 

0 . 1895E+01 

0.1035E-00 

0 . 2245E+00 

0 . 1129E-01 

-0 . 8128E-04 

-0. 1294E-01 

0.0 

-0.1783E-02 

-0 . 1333E+00 

-0.2701E-02 

0. 1944E-04 

0 . 3095E-02 

0 . 6946E-1 4 

0. 4935E-01 

-0.8844E-01 

0 . 1599E-02 

-0 . 1151E-04 

-0.1832E-02 

0.0 

0.0 

0.0 

0 . 9031E-02 

-0.3319E-02 

-0.2117E-01 

0.0 

0.0 

0.0 

0.41B2E-01 

0 . 7658E-0  3 

-0 . X054E+00 

0.0 

0.0 

0.0 

0.5831E-01 

-0 . 4198E-03 

-0.6681E-01 

0 . 3  790E  +  01 

0 . 2069E*00 

0 . 4489E+00 

0. 225BE-01 

-0 . 1626E-03 

-0.2588E-01 

Input  Ramp  Time:  0.3  second 
Command  Vector: 


u  =0.0 

A  =0.0 

np 

q  =0.0 


p  =  3.491  radians/ second  (3  second  pulse) 

r-aT  =  0.0 
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Table  D-8 

Roll  About  Velocity  Vector:  Failed  Model, 
0.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  2.0 


e  =  1.0 

I  =  diag  (  2.0,  1.0,  0.8,  0.25,  0.18,  0.3  > 


0 . 4482E-14  -0.1782E-02  -0. 1332E*00  -0.2703E-02  0.2031E-04  G.3099E-02 

-0.2244E-14  0.2447E-01  -0.5967E-01  0.6049E-02  -0.2047E-02  -0.13S9E-01 

-0.3025E-14  0.2488E-01  -0.2879E-01  -0.4449E-02  0.2035E-02  0.1176E-01 

-0 .  ''669E-16  -0.1463E-03  -0.1094E-01  0.2143E-01  0.1122E-03  -0.5418E-01 

-0. 4928E-1S  -0.5127E-05  -0.3833E-03  0.2918E-01  -0.2193E-03  -0.3346E-01 

0.1895E-01  0.1035E*00  0.2243E«00  0.1130E-01  -0.8493E-04  -0.1296E-01 


0.8965E-14  -0 . 3  564E-02  -0.2665E*00  -0.5405E-02  0.4063E-04  0.6199E-02 

-0.4488E-14  0 . 4893E-01  -0.1193E*00  0.1210E-01  -0.4095E-02  -0.2719E-01 

-0.6051E-14  0.4976E-01  -0.5759E-01  -0.8897E-02  0.4071E-02  0.2352E-01 

-0.1534E-15  -0.2926E-03  -0.2188E-01  0.4286E-01  0.2244E-03  -0.1084E-00 

-0.9855E-15  -0.1025E-04  -0.7666E-03  0.5835E-01  -0.4386E-03  -0.6692E-01 

0 . 3790E+01  0 . 2069E+00  0.4487E«00  0.2260E-01  -0.1699E-03  -0.2592E-01 


Input  Ramp  Time:  0.3  second 


Command  Vector: 


u  =0.0 


A  =0.0 

np 

q  =0.0 

A  =0.0 
YP 

p  =  3.491  radians/ second  (3  second  pulse) 


r-aT  =0.0 
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Figure  D-26: 
Deflections 


0.6M  Failed  Roll  —  Left  Elevator  and  Flaperon 
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Sidef orce  (Flat  Turn)  Maneuver.  Tables  D-9  and  D-10 
give  the  design  parameters  for  both  the  healthy  and  failed 
cases  of  the  sidef orce  maneuver.  The  same  parameters  are 
used  for  both  cases. 

Figures  D-29  through  D-35  present  the  simulation 
responses  for  a  0.2g  sideforce  command.  The  maximum  value 
commanded  is  limited  by  the  canard  deflection.  As  with  the 
roll  maneuver  Figures  D-29  through  D-31  give  simulation 
responses  for  both  the  healthy  and  failed  models. 

For  this  maneuver  the  horizontal  tail  is  mainly  used  to 
counter  the  longitudinal  pitching  moment  produced  by  the  ca¬ 
nard  deflection  although  it  also  contributes  to  roll  control 
for  the  healthy  aircraft  case.  When  the  right  horizontal 
tail  fails  the  left  horizontal  tail  deflects  approximately 
twice  the  deflection  of  the  healthy  elevators.  The  flap- 
erons  deflect  slightly  asymmetrically  to  compensate  for  the 
roll  created  by  the  failed  right  horizontal  tail.  The 
rudder  and  canard  deflections  change  only  slightly  for  the 
failure  case. 


Table  D-9 


Sideforce  Maneuver:  Healthy  Model, 
0.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  = 

1 

.2 

e  = 

1 

.0 

Z  = 

diag  {  2.0 

o 

o 

• 

O 

h-* 

• 

o 

1.0,  1.0  } 

0.0 

-0 . 3565E-04  -0.8330E-01 

-0.5402E-02  0.5400E-04 

0 . 5 1 58E-02 

-0 . 3473E-14 

0 . 9869E-03  -0.5527E-01 

0.3198E-02  -0.3198E-04 

-0.3054E-02 

o 

o 

0.0  0.0 

0.1806E-01  -0.9219E-02 

-0. 3528E-C1 

«0 

- 

0.0 

0.0  0.0 

0.8364E-01  0.2127E-02 

-0 . 1756E+00 

0.0 

0.0  0.0 

0 . 1166E+00  -0.1166E-02 

-0 . 1 1 1 4E*00 

0 . 1895E+01 

0 . 4139E-02  0 . 2B06E+00 

0.4517E-01  -0.4516E-03 

-0.43 1 3E-01 

■ _ 

r~ 

— 

o 

o 

-0 . 4278E-04 

-0.9996E-01 

-0 . 6482E-02 

0.6480E-04 

0.6189E-02 

-0 . 4167E-14 

0 . 1184E-02 

-0 . 6633E-01 

0.3838E-02 

-0.3837E-04 

-0.3665E-02 

0.0 

0.0 

0.0 

0.2168E-01 

-0.1106E-01 

-0.4233E-01 

*1  = 

0.0 

0.0 

0.0 

0 . 1004E+00 

0.2553E-02 

-0 . 2108E*00 

0.0 

0.0 

0.0 

0 . 1400E*00 

-0.1399E-02 

-0.1336E»00 

0.2274E+01 

0 . 4966E-02 

0 . 3367E+00 

0.5420E-01 

-0.5419E-03 

-0.5175E-01 

Input  Ramp  Time:  0.5  second 
Command  Vector: 
u  =0.0 

An  ‘ 

P 

q  =0.0 

A  =  0.2  g  (step) 
yP 

p  =0.0 

r  =  0.01079  radian/second  (step) 


Table  D-10 


Sideforce  Maneuver:  Failed  Model, 
0.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  1-2 

e  =  1.0 

£  =  diag  {  2.0,  0.04,  1.0,  1.0,  1.0,  1.0} 


0. 4482E-14 

-0.7129E-04 

-0.1665E.00 

-0. 1081E-01 

0.  U29E-03 

0.1G33E-0: 

-0.2244E-14 

0 . 9786E-03 

-0.7458E-01 

0.2420E-01 

-0 . 1137E-01 

-0. 4531E-01 

-0.3025E-14 

Q.99S2E-03 

-0.3599E-01 

-0.1779E-01 

0 . 1131E-01 

0. 3920E-01 

-0.7669E-16 

-0.S852E-05 

-0.1367E-01 

0.8571E-01 

0.6232E-03 

-0 . 1806E+00 

-0.4928E-15 

-0 . 2051E-06 

-0.4791E-03 

0.1167E»00 

-0.121BE-02 

-0.1U5E.00 

0.1895E*01 

0.4138E-02 

0 . 2804E.00 

0.4520E-01 

-0.4718E-03 

-0 . 4319E-01 

- 

0.5379E-14 

-0.8554E-04 

-0.1999E.00 

-0.1297E-01 

0.1354E-03 

0.1 2  40E-01 

-0.2693E-14 

0. 1174E-02 

-0.8950E-01 

0.2903E-01 

-0.1365E-01 

-0 . 54  37E-0 1 

-0.3630E-14 

0. 1194E-02 

-0.4319E-01 

-0.2135E-01 

0.1357E-0I 

0 . 4703E-01 

-0 . 9202E-16 

-0.7022E-05 

-0.1641E-01 

0 . 1029E.00 

0 . 7  4  79E-0  3 

-0.2167E.00 

-0.5913E-15 

-0 . 2461E-06 

-0.5749E-03 

0. 1400E*00 

-0.1462E-02 

-0.1338E.00 

0.2274E*01 

0.4966E-02 

0.3365E»00 

0.5424E-01 

-0.5662E-03 

-0 . 5183E-01 

Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =0.0 

A  =  0.2  g  (step) 

YP 

p  =0.0 


r 


0.01079  radian/second  (step) 


Figure  D-29:  0.6M  Healthy/Failed  Sideforce  —  Sideslip 

Angle,  Pitch  Angle,  Angle  of  Attack,  and  Forward  Velocity 


Figure  D-30:  0.6M  Healthy/Failed  Sideforce  —  Roll  Angle 

and  Roll  Rate 


Figure  D-31:  0.6M  Healthy/Failed  Sideforce  —  Normal 

Acceleration,  Lateral  Acceleration,  and  Yaw  Rate 
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Yaw  Pointing  Maneuver .  The  design  parameters  for  the 
yaw  pointing  maneuver  are  given  in  Tables  D-ll  and  D-12  for 
both  the  healthy  and  failed  models.  The  parameters  are  the 
same  for  both  cases. 

Figures  D-36  through  D-42  present  the  simulation 
responses  for  a  four  degree  yaw  pointing  command.  The 
maximum  angle  commanded  for  this  maneuver  is  limited  by  the 
ability  to  hold  the  roll  angle  to  a  small  value.  As  with 
the  roll  maneuver  Figures  D-36  through  D-38  give  simulation 
responses  for  both  the  healthy  and  failed  aircraft  models. 

For  the  yaw  pointing  maneuver  the  horizontal  tail  is 
used  for  both  pitch  and  roll  control.  When  the  right  hori¬ 
zontal  tail  fails  the  flaperons  take  over  the  roll  control 
and  deflect  asymmetrically  to  counter  the  rolling  moment  of 
the  left  horizontal  tail,  which  deflects  for  pitch  control. 
The  rudder  and  canard  responses  are  not  affected  much  by  the 
failure  for  this  maneuver. 


Table  D-ll 


Yaw  Pointing:  Healthy  Model,  0.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

o  =  4.0 
e  =  1.0 

E  =  diag  {  2.0,  0.04,  1.0,  0.7,  2.5,  1.5} 


o 

o 

-0.3565E-04 

-0 . 8330E-01 

-0.3781E-02 

0.1350E-03 

C. 7736E-02 

-0 . 3473E-14 

0 . 9869E-03 

-0.5527E-01 

0 . 2239E-02 

-0 . 7994E-04 

-0.4581E-02 

0.0 

o 

o 

0.0 

0 . 1264E-01 

-0.2305E-01 

-0.5292E-01 

0.0 

0.0 

0.0 

0.S855E-01 

0.5318E-02 

-0.2635E-00 

0.0 

o 

o 

o 

o 

0 . 8164E-01 

-0 . 2915E-02 

-0. 1670E-00 

0.189  5E  +  01 

0 . 4139E-02 

0.2806E+00 

0 . 3162E-01 

-0.U29E-02 

-0 . 6469E-01 

o.o 

-0.1426E-03 

-0 . 3332E+00 

-0.1512E-01 

0 . 5400E-03 

0.3095E-02 

-0.1389E-13 

0 . 3948E-02 

-0.2211E*00 

0.8956E-02 

-0.3198E-03 

-0.1832E-01 

0.0 

o 

o 

0.0 

0.5058E-01 

-0.9219E-01 

•0 . 2117E  +  00 

0.0 

o 

o 

0.0 

0 . 2342E+00 

0 . 2127E-01 

-0.1054E-01 

0.0 

0.0 

0.0 

0 . 3266E+00 

-0.1166E-01 

-0.6681E-C0 

0. 7579E*01 

0.1655E-01 

0. 1122E+01 

0. 1265E*00 

-0.4516E-02 

-0.2588E*00 

Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =  0.0 

A  =  0.0 

n 

P 

q  =  0.0 

A  =  0.0 

P 

p  =  0.0 

r  =  -0.06981  radian/second  (1  second  pulse) 


Table  D-12 


Yaw  Pointing:  Failed  Model,  0.6  Mach  at  30,000  Feet 
Sampling  Time:  T  =  0.02  second 

a  =  4.0 

e  =  1.0 

£  =  diag  (  2.0,  0.04,  1.0,  0.7,  2.5,  1.5> 


0.4482E-14  -0.7129E-04  -0.1665E*00  -0.7567E-02  0.2B21E-03  G.1550E-01 

-0.2244E-14  0.9786E-03  -0.7456E-01  0.1694E-01  -0.2844E-01  -0.6797E-01 

-0.3025E-14  0 . 9952E-03  -0.3599E-01  -0.1246E-01  0.2827E-01  G.5679E-01 

-0.7669E-16  -0.5852E-05  -0.1367E-01  0.6000E-01  0.15S8E-C2  -0.2709E*GC 

-0.4928E-15  -0.2051E-06  -0.4791E-03  0.8169E-01  -0.3046E-C2  -0.1673E-C0 

0. 189SE-01.  0.4138E-02  0.2804E-00  0.3164E-01  -0.1160E-02  -0.6479E-01 


0.1793E-13 

-0.8976E-14 

-0.1210E-13 

-0.3067E-15 

-0.1971E-14 


-0.28S1E-03 

0.3915E-02 

0.3981E-02 

-0.2341E-04 

-0.8203E-06 

0.1655E-01 


-0.6662E*00 
-0.2983E*00 
-0. 1440E*00 
-0 . 5469E-01 
-0.1916E-02 
0.1122E-01 


-0.3027E-01 
0 . 6775E-0 1 
-0 . 4983E-01 
0.2400E-00 
0.3268E-00 


0.1129E-02 
-0.1137E-00 
0.1131E-00 
0 . 6232E-02 
-0.1218E-01 


0.6199E-01 
-0 . 2  719E*0C 
0.2342E*00 
-0. 1084E-C1 
-0.6692E-00 


0 . X265E*00  -0.4718E-02  -0.2592E-00 


Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =  0.0 

A  =  0.0 

np 

q  =  0.0 


p  =  0.0 

r  =  -0.06981  radian/second  (1  second  pulse) 
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Figure  D-38:  0.6M  Healthy/Failed  Yaw  Pointing  —  Norma 

Acceleration,  Lateral  Acceleration,  and  Yaw  Rate 


Figure  D-39:  0.6M  Healthy  Yaw  Pointing  —  Elevator  and  Flap 

Deflections  and  Thrust 


Figure  D-40:  0.6M  Failed  Yaw  Pointing  —  Left  Elevator  and 

Flaperon  Deflections 


Figure  D-42:  0.6M  Failed  Yaw  Pointing  --  Rudder  and  Canard 

Deflections  and  Thrust 
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Lateral  Translation  Maneuver .  Tables  D-13  and  D-14 
give  the  design  parameters  for  the  lateral  translation 
maneuver.  As  indicated  in  the  tables  the  same  design  param 
eters  are  used  for  both  the  healthy  and  failed  aircraft 
models . 

The  simulation  responses  for  a  0.17g  lateral  acceler¬ 
ation  command,  which  corresponds  to  a  steady-state  lateral 
velocity  value  of  8.2  feet/second,  are  given  in  Figures  D-4 
through  D-49.  The  maximum  lateral  acceleration  commanded 
for  this  maneuver  is  limited  by  the  canard  deflection.  As 
explained  for  the  roll  maneuver.  Figures  D-43  through  D-45 
give  simulation  responses  for  both  the  healthy  and  failed 
aircraft  models. 

For  the  lateral  translation  maneuver  the  horizontal 
tail  is  used  mainly  as  an  elevator  to  counter  the  small 
pitch  moment  generated  by  the  canard  deflection.  When  the 
right  horizontal  tail  fails  the  left  horizontal  tail 
deflects  to  approximately  twice  the  elevator  deflection  of 
the  healthy  aircraft  simulation,  and  the  flaperons  deflect 
asymmetrically  to  counter  the  rolling  moment  of  the  left 
horizontal  tail.  The  rudder  and  canard  responses  are  not 
affected  much  by  the  failure  for  this  maneuver. 


Table  D-13 


Lateral  Translation:  Healthy  Model, 
0.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 


a  =  2.0 


e  =  1.0 


£  =  diag  {  2.0,  0.04,  1.0,  3.0,  0.7,  0.7) 


0.0 

-0.3473E-14 


-0.3565E-04 

0.9869E-03 


-0.8330E-01 

-0.SS27E-01 


0.0 


0.0 


0.0  0.0 

0.0  0.0 


0.0 


0.0 


0.0 


0 . 1895E+01 


0.4139E-02  0 . 2806E+00 


-0.1621E-01 
0.9595E-02 
0 . 5419E-01 
0 . 2509E+00 
0. 3499E+00 
0.1355E*00 


0 . 3780E-04 
-0 . 2238E-04 
-0.6453E-02 
0.14B9E-02 
-0.8162E-03 
-0.3161E-03 


0.3610E-02 
-0 . 2138E-02 
-0.2469E-01 
-0.1230E*00 
-0.7795E-01 
-0.3019E-01 


*1 


0.0 

-0. 7130E-04 

-0 . 1666E+00 

-0.3241E-01 

0.7561E-04 

0.7221E-02 

0.6946E-14 

0.1974E-02 

-0 . 1105E*00 

0.1919E-01 

-0.4477E-04 

-0.4275E-02 

0.0 

0.0 

0.0 

0.1084E+00 

-0 . 1291E-01 

-0.4939E-01 

0.0 

0.0 

0.0 

0. 5018E+00 

0.2978E-02 

-0.24S9E»00 

0.0 

0.0 

0.0 

0 . 6998E+00 

-0.1632E-02 

-0.1559E*00 

0.3790E+01 

0 . 8277E-02 

0. 5612E*00 

0.2710E‘00 

-0.6322E-03 

-0.6038E-01 

Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

An  =  0.0 

P 

q  =0.0 

A  =  0.17  g  (1  second  pulse) 
YP 

p  =0.0 

r  =0.0 
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Table  D-14 


Lateral  Translation:  Failed  Model, 
0.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  2.0 
e  =  1.0 

£  =  diag  t  2 . 0 ,  0.04,  1.0,  3.0,  0.7,  0.7} 


0 . 4482E-14 

-0.7129E-04 

-0. 1665E*00 

-0.3243E-01 

0.7900E-04 

0.7232E-02 

-0 . 2244E-14 

0. 9786E-03 

-0.7458E-01 

0.7259E-01 

-0.7962E-02 

-0.3172E-01 

-0.3025E-14 

0.9952E-03 

-0. 3599E-01 

-0 . 5338E-01 

0 . 7915E-02 

0.2744E-01 

-0.7669E-16 

-0.5852E-05 

-0.1367E-01 

0 . 2571E+00 

0.4363E-03 

-0. 1264E-00 

-0.4928E-15 

-0.2051E-06 

-0.4791E-03 

0 . 3501E+00 

-0.8528E-03 

-0.7807E-01 

0.1895E-01 

0.4138E-02 

0. 2804E  +  00 

0. 1356E+00 

-0.3303E-03 

-0.3023E-01 

— 

_ 

— 

__ 

0.8965E-14 

-0.1426E-03 

“0 . 3331E+00 

-0.6486E-01 

0.1580E-03 

0.1446E-01 

-0.4488E-14 

0.1957E-02 

-0.1492E*00 

0.1452E*00 

-0.1592E-01 

-0.6344E-01 

-0.6051E-14 

0.1990E-02 

-0.7199E-01 

-0.1068E*00 

0.1583E-01 

0.5487E-01 

-0.1534E-15 

-0.1170E-04 

-0.2734E-01 

0 . 5143E+00 

0.872SE-03 

-0.252BE*00 

-0.9855E-15 

-0.4101E-06 

-0.9582E-03 

0.7002E*00 

-0.1706E-02 

-0.1561E»00 

0. 3790E*01 

0.8277E-02 

0 . 560BE+00 

0 . 2712E+00 

-0.6605E-03 

-0.6047E-01 

Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

An  =  0.0 

P 

q  =0.0 

A  =  0.17  g  (1  second  pulse) 
YP 

p  =0.0 

r  =0.0 
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Figure  D-43:  0 . 6M  Healthy/Failed  Lateral  Translation  — 

Sideslip  Angle,  Pitch  Angle,  Angle  of  Attack,  and  Forward 
Velocity 


Figure  D-44:  0.6M  Healthy/Failed  Lateral  Translation 

Roll  Angle,  Roll  Rate,  and  Yaw  Rate 


Supersonic  Flight  (1.6  Mach  at  30 , 0Q0  Feet) 

The  supersonic  flight  condition  is  chosen  to  demon¬ 
strate  maneuverabiltiy  of  the  aircraft  at  high  speed.  The 
dynamic  pressure  has  a  large  value  at  this  flight  condition 
because  of  the  high  velocity;  therefore,  designs  are  fairly 
easy  for  the  healthy  aircraft  model  because  the  control 
surfaces  are  very  effective  in  maneuvering  the  aircraft. 
Designs  for  the  failed  right  horizontal  tail  model  are 
slightly  more  difficult  because  the  horizontal  tail  makes  a 
large  contribution  to  roll  control.  When  the  right  horizon¬ 
tal  tail  fails,  the  flaperons  are  required  to  take  over 
complete  roll  control  which  may  lead  to  large  flaperon 
deflections . 

G-Command  Maneuver .  Tables  D-15  and  D-16  list  the 
design  parameters  for  the  g-command  maneuver  for  the  healthy 
and  failed  right  elevator  models.  As  shown  in  the  tables, 
the  same  design  parameters  are  used  for  both  aircraft  models. 

Figures  D-50  through  D-56  show  the  simulation  responses 
for  the  g-command  maneuver.  For  the  same  reasons  given  for 
the  longitudinal  maneuvers  of  the  subsonic  flight  condition, 
Figures  D-50  through  D-52  give  simulation  responses  for  both 
the  healthy  and  failed  aircraft  models. 

The  maximum  normal  acceleration  command  of  8.5  g  is 
limited  by  the  ability  to  hold  the  angle  of  attack  to  a 
relatively  low  value.  For  the  failure  case,  the  left 
elevator  deflects  approximately  twice  the  amount  required 


when  both  tail  halves  are  operational.  The  flaperon 
deflection  is  asymmetric  for  the  failure  case  in  order  to 
compensate  for  the  rolling  moment  created  by  the  left 
elevator  deflection.  The  rudder  and  canard  compensate  for 
the  adverse  yaw  caused  by  a  difference  in  drag  between  the 
left  and  right  halves  of  the  aircraft. 


Table  D-16 


G-Command:  Failed  Model,  1.6  Mach  at  30,000 


Sampling  Time:  T  =  0.02  second 


a  =  3.5 

e  =  1.0 


£  =  diag  f  2.0,  0.021,  1.8,  1.0,  1.0,  1.0* 


^0 


-0.5389E-16 

-0 . 1384E-03 

-0.2388E-01 

0. 3020E-03 

-0 . 2 4 54E-0 5 

-0 

0.0 

0 . 3403E-03 

-0.9S1SE-01 

0.2601E-G1 

-0.U46E-01 

-0 

-0.4012E-15 

0 . 4381E-03 

-0. 7828E-01 

-0 . 2647E-0 1 

0. 1146E-01 

0 

0.2597E-15 

0 . 2030E-04 

0 . 3502E-02 

0 . 3156E-01 

0.7917E-02 

-0 

0.6572E-16 

0 . 1059E-05 

0.1827E-03 

0 . 1767E-0 1 

-0.1436E-03 

-0 

0.2766E+01 

0.161 2E-01 

-0 . 2880E+01 

-0 . 1246E-0 1 

0.101 3E-0  3 

0 

*1 


•0.1886E-1S 

-0.4845E-03 

-0.8357E-01 

0.1057E-02 

-0.8S88E-05 

-0 

0.0 

0.1191E-02 

-0.3330E-00 

0.9103E-01 

-0.4011E-01 

-0 

■0 . 1404E-14 

0.1533E-02 

-0.2740E-00 

-0 . 926  5E-01 

0. 4012E-01 

0 

0.9089E-15 

0 . 7106E-04 

0.1226E-01 

0.1104E*00 

0. 2771E-01 

-0 

0.2300E-15 

0.3707E-05 

0 . 6394E-03 

0.6185E-01 

-0.5025E-03 

-0 

0.9682E»01 

0.5643E-01 

“0 . 1008E*02 

-0. 4363E-01 

0 . 3  544E-03 

0 

Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =0.0 

A  =  8.5  g  (step) 

np 

q  =  0.1719  radian/ second  (step) 

A  =0.0 
YP 

p  =0.0 

r  =0.0 


Feet 


338SE-03 
4627E-0I 
4679E-01 
8413E-0I 
1 9  82E-0 1 
1398E-01 


.  1 1 8  6E-0  2 
.  1619E-00 
.  1638E*00 
.  2944E  +  00 
.  69 3 8E-01 
.  4894E-01 
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Figure  D-51:  1.6M  Healthy/Failed  G-Command  —  Angle  of 

Attack,  Roll  Angle,  Sideslip  Angle,  and  Yaw  Rate 
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Figure  D-52:  1.6M  Healthy/Failed  G-Command  —  P 

Normal  Acceleration,  and  Lateral  Acceleration 


Figure  D-53:  1.6M  Healthy  G-Command  —  Elevator  and  Flap 

Deflections 


Figure  D-54:  1.6m  Failed  G-Command  —  Left  Elevator  and 

Flaperon  Deflections 


Figure  D-55:  1.6M  Healthy  G-Command  —  Aileron,  Rudder 

Canard  Deflections  and  Thrust 


Figure  D-56:  1.6M  Failed  G-Command  —  Rudder  and  Canard 

Deflections  and  Thrust 


Pitch  Pointing  Maneuver .  Tables  D-17  and  D-18  give  the 
design  parameters  for  the  pitch  pointing  maneuver.  Unlike 
the  designs  presented  up  to  this  point,  the  design  parame¬ 
ters  are  not  the  same  for  the  healthy  and  failed  models. 

The  parameters  are  changed  for  the  failure  case  in  order  to 
avoid  overshooting  the  maximum  flaperon  deflection.  Accept¬ 
able  responses  for  the  failure  case  are  achieved  without 
lowering  the  commanded  pitch  pointing  angle. 

The  simulation  responses  for  the  pitch  pointing  maneu¬ 
ver  with  the  healthy  aircraft  model  are  given  in  Figures 
D-57  through  D-61.  Figures  D-62  through  D-66  give  the 
responses  for  the  failure  model.  The  maximum  pointing  angle 
commanded  is  limited  by  the  flaperon  deflection  in  both 
cases.  For  the  failure  model,  the  flaperons  deflect  asymme¬ 
trically  to  counter  the  roll  induced  by  the  left  horizontal 
tail  deflection,  and  the  rudder  and  canard  deflect  to 
compensate  for  adverse  yaw. 

As  shown  in  Figures  D-61  and  D-66  the  thrust  response 
is  negative,  which  is  not  realistic.  The  negative  response 
is  due  to  commanding  large  inputs  to  a  perturbation  model. 
For  the  perturbation  model  a  small  flaperon  deflection  actu¬ 
ally  decreases  drag.  For  large  flaperon  deflections  the 
drag  should  increase,  but  the  model  still  indicates  a  de¬ 
crease  in  drag.  This  decrease  in  drag  is  manifested  by  an 
increase  in  forward  velocity.  The  thrust  is  therefore 
negative  to  maintain  zero  change  in  velocity.  Although  a 
negative  thrust  is  disturbing,  it  does  not  affect  any  of  the 


Table  D-17 


Pitch  Pointing:  Healthy  Model,  1.6 
Sampling  Time:  T  =  0.02  second 
“  =  4.5 

e  =  1.0 

Z  =  diag  f  2.0,  1.5,  2.5,  1.0,  1.0,  1.0 


0.0 

-0 . 4945E-02 

-0. 1658E-01 

0.1530E-03 

0.0 

0 . 2780E-01 

-0. 1204E*00 

-0 . 2361E-03 

0.0 

0.0 

0.0 

0.1531E-01 

0.0 

0.0 

0.0 

0. 3609E-01 

0.0 

0.0 

0.0 

0. 1791E-01 

0.2766E-01 

0 . 1152E  +  01 

-0.4000E*01 

-0 . 1263E-01 

*— 

0.0 

-0.2225E-01 

-0 . 7463E-01 

0.6887E-03 

0.0 

0.1251E-00 

-0 . 5419E+00 

-0.1058E-02 

o 

o 

0.0 

0.0 

0.6892E-01 

0.0 

0.0 

0.0 

0 . 1624E+00 

0.0 

0.0 

0.0 

0.8058E-01 

0 . 1245E+02 

0.5183E+01 

-0  •  1800E+02 

-0.5684E-01 

Mach  at  30,000  Feet 


} 


0.2106E-05 
0 . 3235E-05 
0.6  7 16E-02 
0. 5947E-02 
0 . 2464E-0  3 
0.173BE-03 


-0.9476E-05 
0 . 1456E-04 
-0.3022E-01 
0.2676E-01 
-0 . 1 109E-02 
0.7B21E-03 


-0. 1730E-03 
0 . 2657E-03 
-0. 2720E-01 
-0.921 5E-01 
-0.2024E-01 
0.142BE-01 


-0.7785E-03 
0 . 1196E-02 
-0.1224E»00 
-0.4147E*00 
-0.9109E-01 
0 . 6426E-01 


Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =0.0 

A  =0.0 
nP 

q  =  0.01745  radian/second  (1  second  pulse) 


Table  D-18 


Pitch  Pointing:  Failed  Model,  1.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

ci  =  3.5 
e  =  1.0 

Z  =  diag  {  2.0,  1.5,  1.9,  1.0,  1.0,  1.0} 


-0 . 5389E-1 6 

-0 . 9888E-02 

-0.2520E-01 

0 . 3020E-03 

-0.2454E-05 

-0.3386E-03 

o 

o 

0 . 2431E-01 

-0 . 1004E+00 

0 . 260 1E-0 1 

-0.1146E-01 

-0.4627E-01 

-0.4012E-15 

0 • 3129E-01 

-0 . 8263E-01 

-0 . 2647E-0 1 

0.1 1 46E-01 

0 . 4679E-01 

0.2597E-15 

0 . 1450E-02 

0.3696E-02 

0 . 3156E-01 

0.7917E-02 

-0.841 3E-0 1 

0 . 6572E-16 

0 . 7565E-04 

0 . 1928E-03 

0.1767E-01 

-0.1436E-03 

“0. 1982E-01 

0 . 2766E  +  0 1 

0 . 11 52E  +  01 

-0 . 3040E+0 1 

-0.1246E-01 

0.1013E-03 

0. 1398E-01 

F" -0.1886E-15 

-0.3461E-01 

-0.8822E-01 

0.1057E-02 

-0 . 8S88E-05 

-0.1186E-02 

o 

o 

0 . 8508E-0 1 

-0.3515E*00 

0.9103E-01 

-0.4011E-01 

-0.1619E»00 

-0.1404E-14 

0.1095E*00 

-0.2892E*00 

-0.9265E-01 

0.4012E-01 

0 . 1638E+00 

0.9089E-15 

0.5076E-02 

0. 1294E-01 

0.1104E*00 

0 . 2771E-01 

-0 . 2944E+00 

0.2300E-15 

0.2648E-03 

0 . 6749E-03 

0.6185E-01 

-0.5025E-03 

-0.6938E-01 

0.9682E+01 

0.4031E*01 

-0.1064E»02 

-0.4363E-01 

0 . 3544E-03 

0.4894E-01 

Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =  0.01745  radian/ second  (1  second  pulse) 


Figure  D-63:  1.6M  Failed  Pitch  Pointing  —  Angle  of 

Attack,  Roll  Angle,  Sideslip  Angle,  and  Yaw  Rate 


Figure  D-64:  1.6M  Failed  Pitch  Pointing  —  Normal 

Acceleration  and  Lateral  Acceleration 


Longitudinal  Translation  Maneuver.  Tables  D-19  and 
D-20  give  the  design  parameters  for  the  longitudinal  trans¬ 
lation  maneuver.  Although  the  design  parameters  are  the 
same  for  the  healthy  and  failed  models,  the  commanded  normal 
acceleration  is  reduced  for  the  failed  case  in  order  not  to 
overshoot  the  maximum  flaperon  deflection.  A  change  in  the 
design  parameters  for  the  failure  case  leads  to  unacceptable 
simulation  responses. 

The  simulation  responses  for  the  healthy  model  longitu¬ 
dinal  translation  maneuver  are  given  in  Figures  D-67  through 
D-71.  The  responses  for  the  failure  model  are  given  in 
Figures  D-72  through  D-76.  For  the  failure  model,  the  flap- 
erons  deflect  asymmetrical ly  to  counter  the  roll  induced  by 
the  failed  right  horizontal  tail,  and  the  rudder  and  canard 
deflect  to  compensate  for  adverse  yaw. 


Table  D-20 


Longitudinal  Translation:  Failed  Model, 
1.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  3.5 
e  =  1.0 

Z  =  diag  {  2.0,  1.0,  3.0,  1.0,  1.0,  1.0  J 


-0.5389E-16  -0.6592E-02  -0.3980E-01  0.3020E-03  -0.2454E-05  -0.3388E-03 

0.0  0.1621E-01  -0.1586E.00  0.2601E-01  -0.U46E-01  -0.4627E-01 

-0.4012E-15  0.2086E-01  -0 . 1305E.00  -0.2647E-01  0.1146E-01  0.4679E-01 

0.2597E-15  0.9668E-03  0.5836E-02  0.3156E-01  0.7917E-02  -0.8413E-01 

0 . 6572E-16  0.5044E-04  0.304SE-03  0.1767E-01  -0.1436E-03  -0.19B2E-01 

0 . 2766E+01  0 . 7678E+00  -0 . 4800E»0 1  -0.1246E-01  0.1013E-03  0.1398E-01 


K 


1 


-0.2307E-01 
0.5672E-01 
0.7302E-01 
0.3384E-02 
0 . 1765E-03 
0.2687E»01 


-0.1393E-00 
-0.5550E*00 
-0 , 4566E+00 
0. 2043E-01 
0 . 1066E-Q2 
-0.1680E.Q2 


0.1057E-02 
0.9103E-01 
-0.9265E-01 
0 . 1104E+00 
0.6185E-01 
-0.4363E-03 


-0.8588E-05 
-0.4011E-01 
0.4012E-01 
0. 2771E-01 
-0. 5025E-03 
0. 3544E-03 


Input  Ramp  Time:  0.3  second 
Command  Vector: 
u  =0.0 

=  0.45  g  (1  second  pulse) 
P 

q  =0.0 


p  =0.0 
r  =0.0 
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rqrt,‘ 


Figure  D-70:  1.6M  Healthy  Longitudinal  Translation 

Elevator  and  Flap  Deflections 
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Roil  About  Velocity  Vector  Maneuver .  As  discussed  in 
Chapter  IV  the  roll  maneuver  requires  some  modification  to 
the  healthy  and  failure  models.  Tables  D-21  and  D-22,  give 
the  design  parameters  for  the  roll  maneuver.  The  parameters 
are  the  same  for  the  healthy  and  failed  models;  however,  the 
roll  rate  commanded  for  the  failure  case  is  lowered  from 
that  of  the  healthy  model  in  order  not  to  overshoot  the 
flaperon  deflection  limit. 

The  simulation  responses  for  the  healthy  model  roll 
maneuver  are  given  in  Figures  D-77  through  D-81,  and  the  re¬ 
sponses  for  the  failed  model  are  given  in  Figures  D-82 
through  D-86.  For  the  right  horizontal  tail  case  the  flap- 
erons  deflect  to  achieve  the  commanded  roll  rate. 


Table  D-21 


Roll  About  Velocity  Vector:  Healthy  Model, 
1.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  1.0 

e  =  1.0 

Z  =  diag  2.0,  0.021,  1.8,  1.0,  1.0,  1.0 


0.0 

-0. 6923E-04 

-0.1194E-01 

0. 1530E-03 

-0. 2106E-05 

-0.1730E-03 

o 

o 

0.3892E-03 

-0.8671E-01 

-0.2351E-03 

0 . 3235E-05 

0.2657E-03 

o 

o 

o 

o 

0.0 

0 . 1531E-01 

-0 . 6716E-02 

-0.2720E-01 

0.0 

0.0 

o 

o 

0 . 3609E-01 

0.5947E-02 

-0.9215E-01 

0.0 

o 

o 

o 

o 

0 . 1791E-01 

-0.2464E-03 

-0.2024E-0X 

0. 2766E+01 

0.1612E-01 

-0 . 2880E+01 

-0 . 1263E-01 

0 . 1738E-03 

0.1428E-01 

*1  =  ^0 


Input  Ramp  Time:  0.3  second 
Command  Vector: 
u  =0.0 


P 


0.0 

3.491  radians/second  (3  second  pulse) 


Table  D-22 

Roll  About  Velocity  Vector:  Failed  Model, 
1.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  1.0 


e  =  1.0 

E  =  diag  (  2.0,  0.021,  1.8,  1.0,  1.0,  1.0} 


-0.5389E-16  -0 . 1384E-03  -0.2388E-01  0.3020E-03  -0.2454E-05  -0.3388E-03 

0.3403E-03  -0.9515E-01  0.2601E-01  -0.1146E-01  -0.4627E-01 

-0.4012E-1S  0.4381E-03  -0.7828E-01  -0.2647E-01  0.1146E-01  0.4679E-01 

0.2597E-15  0.2039E-04  0.3502E-02  0.3156E-01  0.7917E-02  -0.8413E-01 

0.6572E-16  0.1059E-05  0.1827E-03  0.1767E-01  -0.1436E-03  -O.1902E-O1 

0.2766E*01  0.1612E-01  -0.2880E-»02  -0.1246E-01  0.1013E-03  0.1398E-01 


h  -  ^0 


Input  Ramp  Time:  0.3  second 


Command  Vector: 


u  =0.0 


A  =0.0 

np 

q  =0.0 


A  =0.0 
YP 

p  =  2.182  radians/second  (3  second  pulse) 


r-oT  =0.0 
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Sidef orce  (Flat  Turn)  Maneuver .  Tables  D-23  and  D-24 
give  the  design  parameters  for  both  the  healthy  and  failed 
cases  of  the  sideforce  maneuver.  The  same  parameters  are 
used  for  both  cases. 

Figures  D-87  through  D-93  present  the  simulation 
responses  for  a  0.6  g  sideforce  command.  The  maximum  value 
commanded  is  limited  by  the  rudder  deflection.  For  this 
maneuver  the  longitudinal  responses  are  very  small  for  both 
the  healthy  and  failure  models.  Also,  the  healthy  aircraft 
lateral  responses  change  very  little  when  the  right  horizon¬ 
tal  tail  fails.  for  these  reasons  Figures  D-87  through  D-89 
give  simulation  responses  for  both  the  healthy  and  failed 
aircraft  models. 

For  this  maneuver  the  horizontal  tail  is  mainly  used  to 
counter  the  longitudinal  pitching  moment  produced  by  the  ca¬ 
nard  deflection  although  it  also  contributes  to  roll  control 
for  the  healthy  aircraft  case.  When  the  right  horizontal 
tail  fails  the  left  horizontal  tail  deflects  approximately 
twice  the  deflection  of  the  healthy  elevators.  The  flap- 
erons  deflect  slightly  asymmetrically  to  compensate  for  the 
roll  created  by  the  failed  right  horizontal  tail.  The 
rudder  and  canard  deflections  change  only  slightly  for  the 
failure  case.  For  the  same  reasons  given  in  the  discussion 
of  the  supersonic  pitch  pointing  maneuver,  a  negative  thrust 
is  generated.  The  negative  thrust  does  not  affect  any  of 
the  other  aircraft  responses  and  is  included  in  the  model 
only  to  hold  forward  velocity  constant. 
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Table  D-23 


n 


Sideforce  Maneuver:  Healthy  Model, 
1.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  3.5 
e  =  1.0 


Z 


diag  (2.0,  0.021,  1.8,  1.0,  1.0,  1.0  } 


y 


0.0 

-0 . 6923E-04 

-0.1194E-01 

0.1530E-03 

-0.2106E-05 

-0.1730E-C3 

0.0 

0 . 3892E-03 

-0. 8671E-01 

-0.2351E-03 

0.3235E-05 

0 . 2657E-03 

0.0 

0.0 

0.0 

0.1531E-01 

-0.6716E-02 

-0 . 2720E-01 

0.0 

0.0 

0.0 

0. 3609E-01 

0 . 5947E-02 

-0.9215E-01 

o 

o 

0.0 

0.0 

0.1791E-01 

-0 . 2464E-03 

-0.2024E-01 

| _ 0 . 2766E  +  01 

0 . 1612E-01 

-0 . 2880E+01 

-0.1263E-01 

0.1738E-03 

0.1428E-01 

0.0 

-0 . 2423E-03 

-0.4179E-01 

0.5356E-03 

-0.7370E-05 

-0.6055E-03 

O 

o 

0.1362E-02 

-0.3035E+00 

-0.8228E-03 

0.U32E-04 

0.9301E-03 

o 

o 

o 

o 

o 

o 

0 . 5360E-01 

-0.23S1E-01 

-0.9518E-01 

o 

o 

0.0 

o 

o 

0 . 1263E+00 

0.2081E-01 

-0 . 3225E+00 

o 

o 

O 

o 

o 

o 

0 . 6268E-01 

-0 . 8624E-03 

-0.7085E-01 

0 . 9682E+01 

0.S643E-01 

-0 . 1008E*02 

-0.4421E-01 

0 . 6083E-03 

0.4998E-01 

Input  Ramp  Time:  0.5  second 
Command  Vector: 
u  =0.0 


q  =0.0 

A  =  0.6  g  (step) 
yP 

p  =0.0 

r  =  0.01214  radian/second  (step) 


*■-  ■'s 

r  v  •.a;^\V.iV«V.>ViW.,i  .ViVi>  II II 


>  --- 
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Table  D-24 


Sideforce  Maneuver:  Failed  Model, 
1.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

o  =  3.5 
e  =  1.0 

Z  =  diag  {  2.0,  0.021,  1.8,  1.0,  1.0,  1.0> 


0 . 5389E-16  -0 . 1384E-03  -0.2388E-01  0.3020E-03  -0.24S4E-05  -0.3388E-03 
0.0  0 . 3403E-03  -0.9515E-Q1  0.2601E-01  -0.1146E-01  -0.4627E-01 
0.4012E-15  0.4381E-03  -0.7828E-01  -0.2647E-01  0.1146E-01  0.4679E-01 
0.2597E-15  0.2030E-04  0.3502E-02  0.3156E-01  0.7917E-02  -0.8413E-01 
0 . 6572E-16  0.1059E-05  0.1827E-03  0.1767E-01  -0.1436E-03  -0.1982E-01 
0 . 2766E+01  0.1612E-01  -0.2880E-01  -0.1246E-01  0.1013E-03  0.1398E-01 


-0.1886E-15  -0.4845E-03  -0.8357E-01  0.1057E-02  -0.8588E-05  -0.11B6E-02 

0.0  0.1191E-02  -0. 3330E+00  0.9103E-01  -0.4011E-01  -0.16X9E*00 

-0.1404E-14  0 . 1533E-02  -0.2740E+00  -0.9265E-01  0.4012E-01  0.1638E*00 

0.9089E-15  0. 7106E-04  0.1226E-01  0.1104E*00  0.2771E-01  -0.2944E-00 

0 . 2300E-15  0 . 3707E-05  0.6394E-03  0.618SE-01  -0.5025E-03  -0.6938E-01 

0 . 9682E+01  0.5643E-01  -0.1008E*02  -0.4363E-01  0.3544E-03  0.4894E-01 


Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =0.0 

A  =  0.6  g  (step) 

YP 

p  =0.0 

r  =  0.01214  radian/ second  (step) 


Yaw  Pointing  Maneuver.  The  design  parameters  for  the 
yaw  pointing  maneuver  are  given  in  Tables  D-25  and  D-26  for 
both  the  healthy  and  failed  models.  The  parameters  are  the 
same  for  both  cases. 

Figures  D-94  through  D-100  present  the  simulation 
responses  for  a  one  degree  yaw  pointing  command.  The 
maximum  angle  commanded  for  this  maneuver  is  limited  by  the 
ability  to  hold  the  roll  angle  to  a  small  value.  As  with 
the  sideforce  maneuver  Figures  D-94  through  D-96  give  simu¬ 
lation  responses  for  both  the  healthy  and  failed  aircraft 
models . 

For  the  yaw  pointing  maneuver  the  horizontal  tail  is 
used  mainly  for  roll  control.  When  the  right  horizontal 
tail  fails  the  flaperons  take  over  the  roll  control.  The 
rudder  and  canard  responses  are  not  affected  much  by  the 
failure  for  this  maneuver.  As  explained  in  the  discussion 
of  the  supersonic  pitch  pointing  maneuver,  the  negative 
thrust  does  not  affect  the  simulation  responses  of  interest. 
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Table  D-25 


Yaw  Pointing:  Healthy  Model,  1.6  Mach  at  30,000  Feet 

Sampling  Time:  T  =  0.02  second 

o  =  5.5 

c  =  1.0 

I  =  diag  {  2.0,  0.021,  1.8,  1.0,  2.0,  3.1} 


- - z)  v  w  / 

0.0 

-0 . 6923E-04 

-0.1194E-01 

0 . 1530E-03 

-0.4212E-05 

-0 . 5363E-03 

0.0 

0 . 3892E-03 

-0 . 8671E-01 

-0 . 2351E-03 

0 . 6470E-05 

0 . 8238E-03 

0.0 

0.0 

0.0 

0 . 1531E-01 

-0 . 1343E-01 

-0 . 8430E-01 

0.0 

0.0 

0.0 

0 . 3609E-01 

0.1189E-01 

“0 . 2857E  +  00 

0.0 

0.0 

0.0 

0 . 1791E-01 

-0.4928E-03 

-0 . 6275E-01 

0 . 2766E+01 

0 . 1612E-01 

-0 . 2880E+01 

-0 . 1263E-01 

0 . 3476E-03 

0 . 4427E-0 1 

0.0 

-0 . 3808E-03 

-0 . 6567E-01 

0 . 8417E-03 

-0.2316E-04 

-0.2950E-02 

0.0 

0 . 2141E-02 

-0 . 4769E+00 

-0 . 1293E-02 

0.3558E-04 

0 . 4531E-02 

0.0 

o 

o 

o 

o 

0 . 8423E-01 

-0.7388E-01 

-0 . 463  7E  +  00 

0.0 

o 

o 

o 

o 

0 . 1985E+00 

0.6542E-01 

-0 . 1571E+01 

0.0 

0.0 

0.0 

0 . 9849E-01 

-0.2710E-02 

-0 . 3451E+00 

0 . 1521E+02 

0 . 8868E-01 

-0 . 1584E+02 

-0 . 6948E-Q1 

0.1912E-02 

0 . 2435E+00 

_ 

— 

Input  Ramp  Time:  0.2  second 

Command  Vector: 

u  =  0.0 

A  =  0.0 

np 

q  =  0.0 

A  =  0.0 

P 


r 


0.0 

0.01745  radian/second  (1  second  pulse) 


Table  D-26 


Yaw  Pointing:  Failed  Model,  1.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  5.5 

e  =  1.0 

Z  =  diag  {  2.0,  0.021,  1.8,  1.0,  2.0,  3.1} 


-0.2964E-15  -0 . 7614E-03  -0.1313E+00  0.1661E-02  -0.2699E-04  -0.5777E-02 

0.0  0.1872E-02  -0 . 5233E+00  0.1430E*00  -0 . 1261E«-00  -0.7888E*00 

-0.2207E-14  0.2410E-02  -0.4305E+00  -0.1456E.00  0.1261E+00  0.7977E+00 

0.1428E-14  0 . 1117E-03  0.1926E-01  0.1736E*00  0.B709E-01  -0.1434E*01 

0.3615E-15  0.5825E-05  0.1005E-02  0.9719E-01  -0.1579E-02  -0.33B0E*00 

0.1521E*02  0.8868E-01  -0.1584E»02  -0.6856E-01  0.1114E-02  0.2384E-»00 


-0.5389E-16  -0.1384E-03  -0.2388E-01  0.3020E-03  -0.490BE-05  -0.1D50E-02 

0.0  0.3403E-03  -0.9515E-01  0.2601E-01  -0.2292E-01  -0.1434E*00 

-0.4012E-15  0.4381E-03  -0.7828E-01  -0.2647E-01  0.2293E-01  0.1450E*00 

0.2597E-15  0.2030E-04  0.3502E-02  0.3156E-01  0.1583E-01  -0.2608E*00 

0.6572E-16  0 . 1059E-05  0.1827E-03  0.1767E-01  -0.2B71E-03  -0.6145E-01 

0.2766E*01  0.1612E-01  -0.2880E*01  -0.1246E-01  0.2025E-03  0.4335E-01 


Input  Ramp  Time:  0.2  second 


Command  Vector: 


r  =  -0.01745  radian/second  (1  second  pulse) 


SO!  1  RMPI  t  'Ctrl  ^  P  (DtP/Stn 
:  c  -c.*c  c.cc  c.*c 


Figure  D-94:  1.6M  Healthy/Failed  Yaw  Pointing  —  Sideslip 

Angle,  Pitch  Angle,  Angle  of  Attack,  and  Forward  Velocity 


Figure  D-95:  1.6M  Healthy/Failed  Yaw  Pointing  —  Roll  Angl 

and  Roll  Rate 


Figure  D-96:  1.6M  Healthy/Failed  Yaw  Pointing  —  Norma 

Acceleration,  Lateral  Acceleration,  and  Yaw  Rate 


Figure  D-98:  1.6M  Failed  Yaw  Pointing  --  Left  Elevator  and 

Flaperon  Deflections 
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Lateral  Translation  Maneuver.  Tables  D-27  and  D-28 
give  the  design  parameters  for  t.he  lateral  translation 
maneuver.  As  indicated  in  the  tables  the  same  design  paranv 
eters  are  used  for  both  the  healthy  and  failed  aircraft 
mode 1 s . 

The  simulation  responses  for  a  0.45  g  lateral  acceler¬ 
ation  command,  which  corresponds  to  a  steady-state  lateral 
velocity  value  of  21.7  feet/ second,  are  given  in  Figures 
D-101  through  D-107.  The  maximum  lateral  acceleration 
commanded  for  this  maneuver  is  limited  by  the  rudder  de¬ 
flection.  As  explained  for  the  sideforce  maneuver,  Figures 
D-101  through  D-103  give  simulation  responses  for  both  the 
healthy  and  failed  aircraft  models. 

For  the  lateral  translation  maneuver  the  horizontal 
tail  is  used  mainly  for  roll  control.  When  the  right 
horizontal  tail  fails  the  flaperons  take  over  complete  roll 
control  The  rudder  and  canard  responses  are  not  affected 
much  by  the  failure  for  this  maneuver. 


Table  D-27 


Lateral  Translation:  Healthy  Model, 
1.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


o  =  3.5 
e  =  1.0 

1  =  diag  {  2.0,  0.021,  1.8,  0.9,  1.8,  2.9  } 


«o 


o.o 

0.0 

0.0 

0.0 

0.0 

0 . 2766E+01 


-0.6923E-04 

0.3892E-03 

0.0 

0.0 

0.0 

0 . 1612E-01 


-0 . 1194E-01 
•0.8671E-01 
0.0 
0.0 
0.0 

■0.2880E+01 


0.1377E-03 
■D.2116E-03 
0.1378E-01 
0 . 3248E-01 
0 . 1612E-01 
0.1137E-01 


-0.3790E-05 

0.5823E-05 

-0.1209E-01 

0.1070E-01 

-0.4435E-03 

0.3129E-03 


■0. 5017E-03 
0.7707E-03 
0.7887E-01 
•0. 2672E  +  00 
0.5B70E-01 
0 . 4141E-01 


*1 


0.0 

0.0 

0.0 

0.0 

0.0 

0. 9682E+01 


-0.2423E-03 

0.1362E-02 

0.0 

0.0 

0.0 

0.5643E-01 


-0.4179E-01 

-0.3035E*00 

0.0 

0.0 

0.0 

-0.1008E+02 


0. 4821E-03 
-0.7406B-03 
0 . 4824E-01 
0.1137E*00 
0.5643E-01 
-0.3979E-0X 


-0.1327E-04 
0.2038E-04 
-0.4231E-01 
0. 3747E-01 
-0.1552E-02 
0. 1095E-02 


-0.1756E-02 
0.2697E-02 
-0 . 2760E*00 
-0.9354E*00 
-0.2055E*00 
0.1449E*00 


Input  Ramp  Time:  0.5  second 
Command  Vector: 
u  =0.0 


q  =  0.0 

A  =  0.45  g  (1  second  pulse) 

yP 

p  =0.0 
r  =  0.0 
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Table  D-28 

Lateral  Translation:  Failed  Model, 
1.6  Mach  at  30,000  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  3.5 


e  =  l.O 


Z  =  diag  (2.0,  0.021,  1.8,  0.9,  1.8,  2.9  } 


-0.5389E-16  -0.1384E-03  -0.2388E-01  0.271BE-03  -0.4417B-05  -0.9826E-03 

0.0  0.3403E-03  -0.9515E-01  0.2341E-01  -0.2063E-01  -0.1342E-00 

-0.4012E-15  0.4381E-03  -0.782BE-01  -0.2383E-01  0.2063E-01  0.1357E*00 

0.2597E-15  0. 2030E-04  0.3502E-02  0.2840E-01  0.1425E-01  -0.2440E*00 

0.6572E-16  0.1059E-05  0.1827E-03  0.1590E-01  -0.25B4E-03  -0.5749E-01 

0.2766E+01  0.1612E-01  -0.2880E*01  -0.1122E-01  0.1823E-03  0.4055E-01 


-0.1886E-1S  -0.4845E-03  -0.83S7E-01  0.9514E-03  -0.1546E-04  -0.3439E-02 

0.0  0.U91E-02  -0.3330E*00  0.B193E-01  -0.7220E-01  -0.46a6E*00 

-0.1404E-14  0.1533E-02  -0.2740E*00  -0.8339E-01  0.7222E-01  0.4749E*00 

0.90B9E-15  0.7106E-04  0.1226E-01  0.9940E-01  0.4988E-01  -0.8539E-00 

0.2300E-15  0.3707E-05  0.6394E-03  0.5566E-01  -0.9045E-03  -0.2012E*00 

0 . 9682E+01  0.5643E-01  -0.1008E*02  -0.3926E-01  0.6380E-03  0.1419E»00 


Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

A  =  0.0 
nP 

q  =0.0 

A  =  0.45  g  (1  second  pulse) 
yP 

p  =0.0 


r  =0.0 
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Figure  D-101:  1.6M  Healthy/Failed  Lateral  Translation  — 

Sideslip  Angle,  Pitch  Angle,  Angle  of  Attack,  and  Forward 
Velocity 
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Figure  D-103:  1.6M  Healthy/Failed  Lateral  Translation 

Normal  Acceleration  and  Lateral  Acceleration 
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Landing  Approach  (0.3  Mach  at  30  Feet) 

The  landing  approach  flight  condition  is  chosen  to 
demonstrate  the  aircraft's  maneuverabiltiy  at  a  low  speed 
and  low  altitude.  Designs  at  this  flight  condition  are  dif¬ 
ficult  due  to  the  low  dynamic  pressure  and  high  trim  angle 
of  attack.  In  reality  the  AFTI/F-16  lands  with  the  flap- 
erons  deflected  fully  downward,  and  the  design  of  control 
laws  requires  much  information  from  pilots  "flying"  in 
simulators.  Because  this  thesis  is  concerned  with  evaluat¬ 
ing  maximum  performance  at  various  flight  conditions  within 
the  flight  envelope,  designs  are  attempted  with  the  flap- 
erons  trimmed  at  zero  degrees,  and  pilot  interfacing  is 
postponed  for  future  research. 

G-Command  Maneuver ■  Tables  D-29  and  D-30  list  the 
design  parameters  for  the  g-command  maneuver  for  the  healthy 
and  failed  right  elevator  models.  As  shown  in  the  tables, 
the  same  design  parameters  are  used  for  both  aircraft  models 

Figures  D-108  through  D-114  show  the  simulation  re¬ 
sponses  for  the  g-command  maneuver.  For  this  maneuver  the 
healthy  aircraft  longitudinal  responses  change  very  little 
when  the  right  horizontal  tail  fails.  Also,  the  healthy 
aircraft  lateral  responses  are  zero  and  are  very  small  with 
the  failed  model.  For  these  reasons  Figures  D-108  through 
D-110  give  simulation  responses  for  both  the  healthy  and 
failed  aircraft  models. 


The  maximum  normal  acceleration  command  of  0.3  g  is 


limited  by  the  flaperon  deflection.  For  the  failure  case, 
the  left  elevator  deflects  approximately  twice  the  amount 
required  when  both  tail  halves  are  working.  The  flaperon 
deflection  is  slightly  asymmetric  for  the  failure  case  in 
order  to  compensate  for  the  rolling  moment  created  by  the 
left  elevator  deflection.  The  rudder  and  canard  compensate 
for  the  adverse  yaw  caused  by  a  difference  in  drag  between 
the  left  and  right  halves  of  the  aircraft. 
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Table  D-29 


G-Command:  Healthy  Model,  0.3  Mach  at  30  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  l.O 
e  =  l.C 


diag  (2.0,  1.0,  1 

f— 

.0,  1.0, 

1.0,  1. 

0  } 

0.0 

0 .  4500E-02 

-0 . 1072E+00 

-0 . 616SE-02 

-0.7253E-04 

0 . 6201E-02 

0.0 

0 .  4782E-01 

-0.9951E-01 

0.7566E-02 

0. 8900E-04 

-0.7609E-02 

0.0 

0.0 

0.0 

0.1034E-01 

-0. 1581E-01 

-0.2613E-01 

0.0 

0.0 

0.0 

0.7860E-01 

0.1089E-01 

-0 . 1751E+00 

0.0 

© 

o 

o 

o 

0. 2146E»00 

0.2525E-02 

*0 . 2158E+00 

0 . 8972E+00 

0 . 3577E.00 

-0 . 1054E+01 

0.3484E-01 

0. 4098E-03 

-0.3504E-01 

<• 


Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =  0.0 

An  =  0.3  g  (step) 

P 

q  =  0.02884  radian/second  (step) 

A  =  0.0 

P 

p  =  0.0 

r  =  0.0 


D-107 


Table  D-30 


1 

i 


i* 


G-Command:  Failed  Model,  0.3  Mach  at  30  Feet 


Sampling  Time:  T  =  0.02  second 

a  =  1.0 

e  =  l.o 

I  =  diag  {  2.0,  1.0,  1.0,  1.0,  1.0,  1.0} 


0.4171E-16 

0.8964E-02 

-0.2136E+00 

-0.1236E-01 

-0.1721E-04 

0 , 1256E-01 

0.0 

0 . 4924E-01 

-0 . 1333E+00 

0.1921E-01 

-0. 2071E-01 

-0. 4001E-01 

0.0 

0. 4645E-01 

-0.6682E-01 

-0.4041E-02 

0 . 2073E-01 

0.2459E-01 

0.1275E-15 

0. 2867E-03 

-0.6833E-02 

0.7887E-01 

0.9883E-02 

-0 . 1764E+00 

0 . 6509E-15 

0.6306E-03 

-0.1503E-01 

0 . 2152E+00 

0 . 2996E-03 

-0.2186E-00 

0 . 8972E+00 

0 . 3578E+00 

-0 . 1056E+01 

0.3493E-01 

0.486  3E-04 

-0 . 3  549E-0 1 

*1 


K 


0 


Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =0.0 

An  =  0.3  g  (step) 

P 

q  =  0.02884  radian/second  (step) 

A  =  0.0 

p  =0.0 

r  =0.0 


a 
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Figure  D-110:  0.3M  Healthy/Failed 

Normal  Acceleration,  and  Lateral 


G-Command  —  Pitch  Rate, 
Acceleration 
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Figure  D-lll:  0.3M  Healthy  G-Command  —  Elevator  and  Flap 
Deflections 


Figure  D-112:  0.3M  Failed  G-Command  —  Left  Elevator  and 

Flaperon  Deflections 


Pitch  Pointing  Maneuver .  Tables  D-31  and  D-32  give  the 
design  parameters  for  the  pitch  pointing  maneuver.  As  with 
the  g-command  maneuver  the  design  parameters  are  the  same 
for  the  healthy  and  failed  models. 

The  simulation  responses  for  the  pitch  pointing  maneu¬ 
ver  are  given  in  Figures  D-115  through  D-121.  The  maximum 
pointing  angle  commanded  (3.5  degrees)  is  limited  by  the 
flaperon  deflection.  For  the  same  reasons  discussed  in  the 
g-command  maneuver  section.  Figures  D-115  through  D-117  give 
simulation  responses  for  both  the  healthy  and  failed  air¬ 
craft  models.  For  the  failure  case  the  left  elevator  must 
deflect  approximately  twice  the  angle  of  the  healthy  eleva¬ 
tor  deflection.  Also,  for  the  failure  model,  the  flaperons 
deflect  asymmetrically  to  counter  the  roll  induced  by  the 
failed  right  horizontal  tail,  and  the  rudder  and  canard 
deflect  to  compensate  for  adverse  yaw. 
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Table  D-31 


Pitch  Pointing:  Healthy  Model,  0.3  Mach  at  30  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  2.0 

c  =  1.0 

E  =  diag{  2.0,  2.0,  1.2,  1.0,  1.0,  1.0 } 


0.0 

0.9000E-02 

-0.1287E*00 

-0.6165E-02 

-0.7253E-04 

0.6201E-02 

0.0 

0 . 9565E-01 

-0.1194E*00 

0.7566E-02 

0.8900E-04 

-0.7609E-02 

0.0 

0.0 

0.0 

0.1034E-01 

-0.1581E-01 

-0.2613E-01 

0.0 

0.0 

0.0 

0.7860E-01 

0.1089E-01 

-0.1751E‘00 

0.0 

0.0 

0.0 

0.2146E+00 

0.2S25E-02 

-0.2158E‘00 

0.8972E+00 

0. 7153E+00 

-0.1264E*01 

0.3484E-01 

0.4098E-03 

*0 . 3504E-01 

— 

— 

_ 

- 

0.0 

0. 1800E-01 

-0.2S74E*00 

-0.1233E-01 

-0.1451E-03 

0.1240E-01 

0.0 

0.1913E+00 

-0.2388E*00 

0.1513E-01 

0.1780E-03 

-0.1522E-01 

0.0 

0.0 

0.0 

0.2067E-01 

-0.3161E-01 

-0.5225E-01 

0.0 

0.0 

0.0 

0.1572Z*00 

0. 21 79E-01 

-0. 3502E+00 

0.0 

0.0 

0.0 

0.4292E»00 

0.5049E-02 

-0.4317E*00 

0.1794E+01 

0.1431E+01 

-0 . 2529E  +  01 

0.6968E-01 

0.8196E-03 

-0.7007E-01 

Input  Ramp  Time:  0.3  second 
Command  Vector: 
u  -0.0 

\ '  °'° 

q  =  0.04363  radian/ second  (1  second  pulse) 

A  =  0.0 

p  =0.0 

r  =  0.0 
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Table  D-32 


Pitch  Pointing:  Failed  Model,  0.3  Mach  at  30  Feet 


Sampling  Time:  T  =  0.02  second 


o  =  2.0 
E  =  1.0 

z  =  diag  *2.0,  2.0,  1.2,  1.0,  1.0,  l.o) 


-0.4171E-16  0 . 1793E-01  -0.2563E*0Q  -0.1236E-01  -0.1721E-04  0.1256E-01 

0.0  0.9848E-01  -0.1$99E»00  0.1921E-01  -0.207XB-01  -0.4001E-01 

0.0  0.9290E-0X  -0.8019E-01  -0.4041E-02  0.2073E-01  0.2459E-01 

0.1275E-15  0.5775E-03  -0.8199E-02  0.7BB7E-01  0.9883E-02  -0.1764E*00 

0.6509E-1S  0.1261E-02  -0.1803E-01  0.2152E*00  0.2996E-03  -0.2186E»00 

0.8972E«00  0. 7155E+00  -0.1267E+0X  0.3493E-01  0.4863E-04  -0.3549E-01 


= 


0.8343E-16  0. 3586E-01  -0.5127E+00  -0.2473E-01  -0.3442E-04  0.2512E-01 
0.0  0.1970E*00  -0.3198E+00  0.3843E-01  -0.4143E-01  -0.8001E-01 
0.0  0.1858E«00  -0.1604E+00  -0.8083E-02  0.4147E-01  0.491BE-01 
0.2551E-15  0.1147E-02  -0.1840E-01  0.1577E*00  0.1977E-01  -0.3527E*00 
0.1302E-14  0.2522E-02  -0.3607E-01  0.4305E*00  0.5991E-03  -0.4372E»00 
0.1794E*01  0. 1431E+01  -0.2535E*01  0.6986E-01  0.9725E-04  -0.7098E-01 


Input  Ramp  Time:  0.3  second 
Command  Vector: 


u  =0.0 


A  =0.0 
"p 


q  =  0.04363  radian/second  (1  second  pulse) 


A  =  0.0 

p  =  0.0 


r  =0.0 
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Figure  D-117:  0.3M  Healthy/Failed  Pitch  Pointing  —  Normal 

Acceleration  and  Lateral  Acceleration 
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Figure  D-119:  0.3M  Failed  Pitch  Pointing  —  Left  Elevator 

and  Flaperon  Deflections 
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Figure  D-120:  0.3M  Healthy  Pitch  Pointing 

Rudder,  and  Canard  Deflections  and  Thrust 


Aileron 


Rudder  and 


Longitudinal  Translation  Maneuver.  Tables  D-33  and 


D-34  give  the  design  parameters  for  the  longitudinal  trans¬ 
lation  maneuver.  Once  again  the  design  parameters  are  the 
same  for  the  healthy  and  failed  models. 

The  simulation  responses  for  the  longitudinal  transla¬ 
tion  maneuver  are  given  in  Figures  D-122  through  D-128.  As 
discussed  in  the  g-command  maneuver  section  Figures  D-122 
through  D-124  give  simulation  responses  for  both  the  healthy 
and  failed  aircraft  models.  The  maximum  longitudinal  trans¬ 
lation  command  of  0.2  g  for  one  second,  which  yields  a 
normal  velocity  of  approximately  9.0  feet/second,  is  limited 
by  the  flaperon  deflection. 

As  is  usually  the  case  with  longitudinal  maneuvers,  the 
failure  model  requires  that  the  left  elevator  must  deflect 
approximately  twice  the  angle  of  the  healthy  elevator  de¬ 
flection.  Also,  for  the  failure  model,  the  flaperons 
deflect  asymmetrically  to  counter  the  roll  induced  by  the 
failed  right  horizontal  tail,  and  the  rudder  and  canard 
deflect  to  compensate  for  adverse  yaw. 
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Table  D-33 


Longitudinal  Translation:  Healthy  Model, 
0.3  Mach  at  30  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  2.0 
e  =  1.0 

E  =  di ag  ^  2.0,  3.0,  1.0,  1.0,  1.0,  1.0  > 


0_q  0.1350E-01  -0 . 1072E+00  -0.6165E-02  -0.7253E-04  0.6201E-02 
O.o  0 . 1435E+00  -0.9951E-01  0.7566E-02  0.8900E-04  -0.7609E-02 
O.o  0.0  0.0  0. 2034E-01  -0.1581E-01  -0.2613E-01 
O.o  0.0  0.0  0.7860E-01  0.1089E-01  -0.1751E»00 
O.o  0.0  0.0  0.2146E»00  0.2525E-02  -0.2158E*00 
0 . 8972E+00  0 . 1073E+01  -0.1054E»01  0.3484E-01  0.4098E-03  -0.3504E-01 


0.0  0.2700E-01  -0.2145E*00  -0.1233E-01  -0.1451E-03  0.1240E-01 
O.o  0 . 2869E+00  -0 . 1990E+00  0.1513E-01  0.1780E-03  -0.1522E-01 
0_0  0.0  0.0  0.2067E-01  -0.3161E-01  -0.5225E-01 
00  0.0  0.0  0.1572E+00  0.2179E-01  -0.3502E*00 
O.o  0.0  0.0  0 . 4292E+00  0.5049E-02  -0.4317E*00 
0 . 1794E*01  0 . 2146E+01  -0.2146E«01  0.6968E-01  0.8196E-03  -0.7007E-01 


Input  Ramp  Time:  0.4  second 
Command  Vector: 
u  =0.0 

A  =  0.2  g  (1  second  pulse) 
P 

q  =0.0 

A  =0.0 
YP 

p  =0.0 

r  =0.0 


D-121 


Table  D-34 


Longitudinal  Translation:  Failed  Model 
0.3  Mach  at  30  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  2.0 


e  =  1.0 


E  =  diag  {  2.0,  3.0,  1.0,  1.0,  1.0,  1.0  ) 


«0 


-0 . 4171E-16 

0.2689E-01 

-0.2136E»00 

-0 . 1236E-0 1 

-0.1721E-04 

0.0 

0  - 1477E+00 

-0 .  133  3E  +  00 

0.1921E-01 

-0.2071E-01 

0.0 

0.1394£*00 

-0.6682E-01 

-0. 4041E-02 

0.2073E-01 

0 . 1275E-15 

0.8602E-03 

-0.6833E-02 

0 . 7887E-01 

0.9BB3E-02 

0 . 6509E-15 

0.1892E-02 

-0.1503E-01 

0 . 2152E+00 

0.2996E-03 

0 . 8972E  +  00 

0.1073E.01 

-0.1056E.0l 

0. 3493E-01 

0. 4863E-04 

*1 


0 . B343E-16 

0.5378E-01 

-0 . 4272E.00 

-0 . 2473E-01 

-0.3442E-04 

0 

0.0 

0 . 2954E+00 

-0.2665E.00 

0.3843E-01 

-0.4143E-01 

-0 

0.0 

0.27B7E.00 

-0 . 1336E+00 

-0 . 8083E-02 

0, 4147E-01 

0 

0.2551E-15 

0.1720E-02 

-0. 1367E-01 

0.1577E.00 

0.1977E-01 

-0 

0 . 1302E-14 

0 . 3784E-02 

-0.3005E-01 

0 . 4305E.00 

0.5991E-03 

-0 

0 . 1794E+01 

0 . 2147E.01 

-0 . 2112E.01 

0.6986E-01 

0.9725E-04 

-0 

Input  Ramp  Time:  0.4  second 
Command  Vector: 
u  =  0.0 

=  0.2  g  (1  second  pulse) 
P 

q  =0.0 


p  =0.0 
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•1256E-01 
.4001E-01 
.  2459E-01 
.  1764E  +  00 
.  2186E  +  00 
3549E-D1 


.2512E-01 
.  800XE-01 
.  49 18E-01 
.  3527E+00 
.  4372E  +  00 
.7098E-01 
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Figure  D-122:  0.3M  Healthy/Failed  Longitudinal  Translation 

—  Pitch  Angle,  Angle  of  Attack,  and  Forward  Velocity 


Figure  D-123:  0 . 3M  Healthy/Failed  Longitudinal  Translation 

—  Flight  Path  Angle,  Roll  Angle,  Sideslip  Angle,  and  Yaw 
Rate 
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Figure  D-124:  0.3M  Healthy/Failed  Longitudinal  Translation 

—  Pitch  Rate,  Normal  Acceleration,  and  Lateral 
Acceleration 


Figure  D-125:  0.3M  Healthy  Longitudinal  Translation 

Elevator  and  Flap  Deflections 
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Roll  About  Velocity  Vector  Maneuver .  As  discussed  in 
Chapter  IV  the  roll  maneuver  requires  some  modification  to 
the  healthy  and  failure  models.  Tables  D-35  and  D-36,  give 
the  design  parameters  for  the  roll  maneuver.  The  parameters 
are  the  same  for  the  healthy  and  failed  models;  however,  the 
roll  rate  commanded  for  the  failure  case  is  lowered  from 
that  of  the  healthy  model  in  order  not  to  overshoot  the 
flaperon  deflection  limit. 

The  simulation  responses  for  the  healthy  model  roll 
maneuver  are  given  in  Figures  D-129  through  D-133,  and  the 
responses  for  the  failed  model  are  given  in  Figures  D-134 
through  D-138.  For  the  right  horizontal  tail  case  the  flap- 
erons  deflect  to  achieve  the  commanded  roll  rate. 
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Table  D-35 


Roll  About  Velocity  Vector:  Healthy  Model, 
0.3  Mach  at  30  Feet 


Sampling  Time:  T  =  0.02  second 
a  =  2.5 

e  =  1.0 

£  =  diag  (  2.0,  1.0,  1.0,  1.7,  2.0,  3.3  > 


*o  = 


o.o 

0.0 

0.0 

0.0 

0.0 

0 . 8972E+00 


0 . 4500E-02  -0.1072E*00  -0.1048E-01  -0.1451E-03 


0.4782E-01 

0.0 

0.0 

0.0 


-0.9951E-01 

0.0 

0.0 

0.0 


0.3S77E*00  -0.1054E+01 


0.1286E-01 
0.1757E-01 
0 . 1336E+00 
0 . 3648E+00 
0.5922E-01 


0.1780E-03 

-0.3161E-01 

0.2179E-01 

0.5049E-02 

0.8196E-03 


*1 


0.0 

0.0 

0.0 

0.0 

0.0 

0.2243E*01 


0.U25E-01  -0.2681E*00  -0.2620E-01  -0.3626E-03 
O.U96E.OO  -0.2488E*00  0.3215E-01  0.4450E-03 
0.0  0.0  0.4393E-01  -0.7903E-01 
0.0  0.0  0.3341E*00  0.5447E-01 
0.0  0.0  0 . 9121E.00  0 . 1262E-01 
0 . 8942E+00  -0 . 2634E*01  0.1481E»00  0.2049E-02 


Input  Ramp  Time:  0.3  second 
Command  Vector: 
u  =0.0 


n 


=  0.0 

q  =0.0 

A  =0.0 

P 

p  =  1.222  radians/ second  (3  second  pulse 


r-°T  =  0.0 


0.20 46E-01 
-0.25UE-01 
-0.8622E-01 
-0 . 5778E+00 
-0 . 7122E+00 
-0 . 1 156E+00 


0.5115E-01 
-0.6277E-01 
-0.21S5E*00 
-0 . 1444E*01 
-0.1781E*01 
-0 . 2890E+00 


) 
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Table  D-36 


Roll  About  Velocity  Vector:  Failed  Model, 
0.3  Mach  at  30  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  2.5 


e  =  1.0 

Z  =  diag  {  2.0,  1.0,  1.0,  1.7,  2.0,  3.3} 


«0 


-0.4171E-16 

0.0 

0.0 

0.1275E-15 

0.6509E-15 

0.8972E-00 


0 . 9964E-02 
0.4924E-01 
0 . 464SE-01 
0.2867E-03 
0 . 6306E-03 
0.3578E*00 


-0.2136E*00 
-0 . 1333E+00 
-0.6682E-01 
-0.6833E-02 
-0.1503E-01 
“0. 1056E+01 


-0.2102E-01 
0.3266E-01 
-0.6870E-02 
0 . 1341E+00 
0.36SBE*00 
0.5938E-01 


-0.3442E-04 

-0.4143E-01 

0.4147E-01 

0.1977E-01 

0.5991E-03 

0.9725E-04 


0.4145E-01 
-0.1320E*00 
0.8U5E-01 
-0.5B20E*00 
-0. 7215E+00 
-0.H71E*00 


- 

-0 . 1043E-15 

0.2241E-01 

-0.5340E»00 

-0.5255E-01 

-0.8606E-04 

0. 1036E+00 

0.0 

0.1231E«00 

-0 . 3332E+00 

0. 8166E-01 

-0.1036E‘00 

*0 . 3300E+00 

0.0 

0. 1161E*00 

-0.1671E+00 

-0.1718E-01 

0.1037E-00 

0 . 2029E*00 

0.3188E-15 

0.7168E-03 

-0.1708E-01 

0. 3352E+00 

0.4941E-01 

-0.1455E»01 

0.1627E-14 

0.1576E-02 

-0.3757E-01 

0 . 9146E+00 

0.149BE-02 

-0 . 1B04E*01 

0. 2243E*01 

0 . 8944E+00 

-0. 2640E+01 

0.1485E*00 

0.2431E-03 

-0.2928E*00 

Input  Ramp  Time:  0.3  second 
Command  Vector: 
u  =0.0 


0.0 

0.0 

0.0 

1.047  radians/second  (3  second  pulse) 


ROLi 


Figure  D-133:  0.3M  Healthy  Roll  —  Aileron,  Rudder,  and 

Canard  Deflections 


Figure  D-134:  0.3M  Failed  Roll  —  Sideslip  Angle,  Pitch 

Angle,  Angle  of  Attack,  and  Forward  Velocity 
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Figure  D-136:  0.3M  Failed  Roll  —  Normal  Acceleration  and 

Lateral  Acceleration 
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Figure  D-138:  0.3M  Failed  Roll  —  Rudder  and  Canard 

Deflections  and  Thrust 
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Sidef orce  (Flat  Turn)  Maneuver .  Tables  D-37  and  D-38 
give  the  design  parameters  for  both  the  healthy  and  failed 
cases  of  the  sideforce  maneuver.  The  same  parameters  are 
used  for  both  cases. 

Figures  D-139  through  D-145  present  the  simulation 
responses  for  a  0.1  g  sideforce  command.  The  maximum  value 
commanded  is  limited  by  the  canard  deflection.  For  this 
maneuver  the  longitudinal  responses  are  very  small  for  both 
the  healthy  and  failure  models.  Also,  the  healthy  aircraft 
lateral  responses  change  very  little  when  the  right  horizon¬ 
tal  tail  fails.  For  these  reasons  Figures  D-139  through 
D-141  give  simulation  responses  for  both  the  healthy  and 
failed  aircraft  models. 

For  this  maneuver  the  horizontal  tail  is  mainly  used  to 
counter  the  longitudinal  pitching  moment  produced  by  the  ca¬ 
nard  deflection  although  it  also  contributes  to  roll  control 
for  the  healthy  aircraft  case.  When  the  right  horizontal 
tail  fails  the  left  horizontal  tail  deflects  approximately 
twice  the  deflection  of  the  healthy  elevators.  The  flap- 
erons  deflect  asymmetrically  to  compensate  for  the  roll 
created  by  the  failed  right  horizontal  tail.  The  rudder  and 
canard  deflections  change  only  slightly  for  the  failure 
case. 
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Table  D-37 


Sideforce  Maneuver:  Healthy  Model, 
0.3  Mach  at  30  Feet 

Sampling  Time:  T  =  0.02  second 

o  =  l.o 

e  =  1.0 

I  =  diag  {  2.0,  1.0,  1.0,  1.0,  1.0,  1.0} 


Input  Ramp  Time:  0.5  second 

Command  Vector: 

u  =0.0 

A  =0.0 
P 

q  =0.0 

A  =  0.1  g  (step) 

YP 

p  =0.0 

r  =  0.009641  radian/second  (step) 


Table  D-38 


1 


Sideforce  Maneuver:  Failed  Model, 
0.3  Mach  at  30  Feet 


Sampling  Time: 

T  =  0. 

02  second 

a 

=  1 

..0 

e 

=  1 

..0 

E 

=  diag  t  2.0, 

1.0,  1. 

o 

H* 

• 

o 

H-* 

• 

o 

H* 

• 

O 

f 

-0.4171E-16 

0.8964E-02 

-0.2136E+00 

-0.1236E-01 

-0.1721E-04 

0. 1256E-01 

0.0 

0.4924E-01 

-0 . 1333E+00 

0.1921E-01 

-0.2071E-01 

-0.4001E-01 

*0 

= 

0.0 

0.1275E-15 

0.4645E-01 

0.2867E-03 

-0.6682E-01 

-0.6833E-02 

-0.4041E-02 

0. 7BB7E-01 

0 . 2073E-01 

0.9883E-02 

0 , 2459E-0 1 

-0 . 1 764E+00 

0.6509E-1S 

0.6306E-03 

-0.1503E-01 

0.2152E+00 

0 . 2996E-03 

-0 . 2186E+00 

0. 8972E+00 

0.3578E*00 

-0. 1056E+01 

0.3493E-01 

0.4863E-04 

-0 . 3549E-01 

Input  Ramp  Time:  0.5  second 
Command  Vector: 
u  =0.0 

\  *  °-° 
q  =0.0 

A  =  0.1  g  (step) 
yP 

p  =0.0 

r  =  0.009641  radian/second  (step) 
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Figure  D-141:  0.3M  Healthy/Failed  Sideforce  —  Normal 

Acceleration,  Lateral  Acceleration,  and  Yaw  Rate 


^Deflections  aSd^hrus^  S1<JeEorce  "  and  Canard 
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Yaw  Pointing  Maneuver.  The  design  parameters  for  the 
yaw  pointing  maneuver  are  given  in  Tables  D-39  and  D-40  for 
both  the  healthy  and  failed  models.  The  parameters  are  the 
same  for  both  cases. 

Figures  D-146  through  D-152  present  the  simulation 
responses  for  a  three  degree  yaw  pointing  command.  The 
maximum  angle  commanded  for  this  maneuver  is  limited  by  the 
ability  to  hold  the  roll  angle  to  a  small  value.  As  with 
the  sideforce  maneuver  Figures  D-146  through  D-148  give 
simulation  responses  for  both  the  healthy  and  failed  air¬ 
craft  models. 

For  the  yaw  pointing  maneuver  the  horizontal  tail  is 
used  for  both  pitch  and  roll  control.  When  the  right  hori¬ 
zontal  tail  fails  the  flaperons  take  over  the  roll  control 
and  deflect  slightly  asymmetrically  to  counter  the  rolling 
moment  of  the  left  horizontal  tail,  which  deflects  for  pitch 
control.  The  rudder  and  canard  responses  are  not  affected 
much  by  the  failure  for  this  maneuver. 


Table  D-39 


Yaw  Pointing:  Healthy  Model,  0.3  Mach  at  30  Feet 

Sampling  Time:  T  =  0.02  second 

a  =  3.0 
e  =  1.0 

Z_  =  diag  {  2.0,  1.0,  1.0,  1.0,  3.0,  2.0} 


K 


1 


0.1350E-01 
0 . 1435E+00 
0.0 
0.0 
0.0 

0.1073E*01 


-0.321 7E*00 
-0.29B5E*00 
0.0 
0.0 
0.0 

-0. 3161E-01 


-0.1850E-03 
0.2270E-01 
0 . 3101E-01 
0 . 2358E-00 
0. 6438E+00 
0 . 1045E+00 


-0.6527E-03 

0.8010E-03 

-0.1422E-00 

0.9B05E-01 

0.2272E-01 

0.3688E-02 


Input  Ramp  Time:  0.2  second 
Command  Vector: 
u  =  0.0 

v  °-° 

q  =0.0 

A  =  0.0 

P 

p  =  0.0 


0.05236  radian/second  (1  second  pulse) 
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Table  D-40 


Yaw  Pointing:  Failed  Model,  0.3  Mach  at  30  Feet 


Sampling  Time:  T  =  0.02  second 


a  =  3.0 


e  =  1.0 


E  =  diag  {2.0,  1.0,  1.0,  1.0,  3.0,  2.0} 


*0 


— 

— 

-0 . 41 7  IE-  If) 

0 . 8964E-02 

-0 . 2136E+00 

-0.1236E-01 

-0.S164E-04 

0.2512E-01 

0.0 

0.492  4E-01 

-0.1333E+00 

0.1921E-01 

-0 . 6214E-01 

-0.8001E-01 

0.0 

0 . 4645E-01 

-0 . 66B2E-01 

-0 . 4041E-02 

0.6220E-01 

0.4918E-01 

0.1275E-1S 

0.2867E-03 

-0 . 6833E-02 

0.7887E-01 

0.2965E-01 

-0 . 3527E+00 

0.6509E-15 

0. 6306E-03 

-0.1503E-01 

0.2152E+00 

0.6987E-03 

-0 . 4372E+00 

0 . 8972E+00 

0 . 3578E+00 

-0 . 1056E+01 

0.3493E-01 

0. 1459E-03 

-0.7098E-01 

K 


1 


-0.1251E-1S 

O.2689E-01 

-0.6408E*00 

-0.3709E-01 

-0.1549E-03 

0.7537E-01 

0.0 

0.1477E*00 

-0 . 399BE+00 

0.5764E-01 

-0 . 1864E+00 

-0.2400E+00 

0.0 

0.1394E+00 

-0.200SE*00 

-0.1212E-01 

0 . 1866E+00 

0.1476E+00 

0.3826E-15 

0.8602E-03 

-0.2050E-01 

0 . 2366E+00 

0.8895E-01 

-0.105BE+01 

0.1953E-14 

0.1892E-02 

-0.4508E-01 

0.6456E*00 

0.2696E-02 

-0 . 1312E*01 

0 . 2692E+0 1 

0 . 1073E+0X 

-0 . 3168E+01 

0. 1048E+00 

0.4376E-03 

-0.2129E+00 

— 

_ 

Input  Ramp  Time:  0.2  second 
Command  Vector: 
u  =  0.0 

v  °-° 

q  =  0.0 

A  =  0.0 

p  =  0.0 

r  =  -0.05236  radian/second  (1  second  pulse) 
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Figure  D-148:  0.3M  Healthy/Failed  Yaw  Pointing  — 

Acceleration,  Lateral  Acceleration,  and  Yaw  Rate 


Figure  D-150:  0.3M  Failed  Yaw  Pointing  —  Left  Elevator  and 

Flaperon  Deflections 
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Figure  D-152:  0.3M  Failed  Yaw  Pointing  —  Rudder  and  Canard 

Deflections  and  Thrust 


D-148 


Lateral  Translation  Maneuver.  Tables  D-41  and  D-42 
give  the  design  parameters  for  the  lateral  translation 
maneuver.  As  indicated  in  the  tables  the  same  design  param 
eters  are  used  for  both  the  healthy  and  failed  aircraft 
models . 

The  simulation  responses  for  a  0.1  g  lateral  acceler¬ 
ation  command,  which  corresponds  to  a  steady-state  lateral 
velocity  value  of  4.5  feet/second,  are  given  in  Figures 
D-153  through  D-159.  The  maximum  lateral  acceleration 
commanded  for  this  maneuver  is  limited  by  the  canard  deflec 
tion.  As  explained  for  the  sideforce  maneuver,  Figures 
D-153  through  D-155  give  simulation  responses  for  both  the 
healthy  and  failed  aircraft  models. 

For  the  lateral  translation  maneuver  the  horizontal 
tail  is  used  mainly  as  an  elevator  to  counter  the  small 
pitch  moment  generated  by  the  canard  deflection.  When  the 
right  horizontal  tail  fails  the  left  horizontal  tail 
deflects  to  approximately  twice  the  elevator  deflection  of 
the  healthy  aircraft  simulation,  and  the  flaperons  deflect 
asymmetrically  to  counter  the  rolling  moment  of  the  left 
horizontal  tail.  The  rudder  and  canard  responses  are  not 
affected  much  by  the  failure  for  this  maneuver. 


Table  D-41 


Lateral  Translation:  Healthy  Model, 
0.3  Mach  at  30  Feet 

Sampling  Time:  T  =  0.02  second 


a  =  1.0 


e  =  1.0 


£  =  diag<  2.0,  1.0,  1.0,  2.5,  1.5,  1.0  } 


*0 


0.0 

0.4500E-02 

“0. 1072£*00 

-0.1541E-01 

-0.1D88E-03 

0.6201E-02^ 

0.0 

0 . 4782E-01 

-0.9951E-01 

0.1891E-01 

0.1335E-03 

-0 . 7609E-02 

0.0 

0.0 

0.0 

0.2584E-01 

-0.2371E-01 

-0.2613E-01 

0.0 

0.0 

0.0 

0 . 1965E+00 

0.1634E-01 

-0 . 1751E+00 

0.0 

0.0 

0.0 

0 . 5365E+00 

0. 37B7E-02 

-0 . 2158E+00  ! 

I 

0.8972E*00 

0 . 3577E+00 

*0 . 1054E+01 

0.8709E-01 

0.6147E-03 

-0.3504E-01  | 

K 


1 


K 


0 


Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =  0.0 

An  =  0.0 

P 

q  =  0.0 

Ay  =  0.1  g  (1  second  pulse) 


Table  D-42 


Lateral  Translation:  Failed  Model, 
0.3  Mach  at  30  Feet 

Sampling  Time:  T  =  0.02  second 

£  =  l.O 

E  =  1.0 

l  =  diag  {  2.0,  1.0,  1.0,  2.5,  1.5,  1.0} 


0.8964E-02 
0 . 4924E-0 1 


0.0  0.4645E-01 

0.1275E-15  0.2867E-03 


0.6509E-15  0.6306E-03 


0.3578E+00 


-0. 2136E+00 
•0 . 1333E+00 
-0.6682E-01 
-0.6833E-02 
-0.1503E-01 
-0.1056E*01 


-0.3091E-01 
0.4B03E-01 
-0.1010E-01 
0 . 1972E+00 
0. 5380E+00 
0.8733E-01 


-0.2582E-04 

-0.3107E-01 

0.3110E-01 

0.1482E-01 

0.4493E-03 

0.7294E-04 


0.1256E-01 
-0.4001E-01 
0.2459E-01 
-0.1764E+00 
-0 . 2186E+00 
-0.3549E-01 


Input  Ramp  Time:  0.4  second 

Command  Vector: 

u  =0.0 

A  =0.0 

np 

q  =0.0 

A  =  0.1  g  (1  second  pulse) 
YP 

p  =0.0 

r  =  0.0 
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Q.4C 


>)  *.  qr  io 

.02  C.04  0.06  0.08  O.tO 


Figure  D-155:  0.3M  Healthy/Failed  Lateral  Translation  — 

Normal  Acceleration  and  Lateral  Acceleration 
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control  laws  are  tailored  to  perform  seven  maneuvers  at  four  flight  conditions. 

Maximum  maneuvers  are  commanded  to  yield  maximum  control  surface  deflections. 

The  same  control  law  designs  are  then  applied  to  the  model  with  a  failed  right 
horizontal  tail,  and  the  performance  is  evaluated.  Some  maneuvers  require  modifi¬ 
cations  to  the  designs  or  lowered  maximum  maneuver  requirements  to  avoid  overshooting 
the  deflection  limits  of  the  operational  control  surfaces. 

Simulation  responses  are  presented  for  both  the  healthy  and  failure  aircraft 
models.  Generally,  when  the  right  horizontal  tail  fails,  the  left  horizontal  tail 
assumes  primary  pitch  control  and  the  flaperons  take  over  complete  roll  control.  The 
flaperons,  rudder,  and  canards  deflect  to  counter  the  rolling  and  yawing  moments  pro¬ 
duced  by  the  left  horizontal  tail  deflection. 

The  design  process  used  in  the  study, -is  applicable  to  other  control  surface 
failures.  Recommendations  are  presented  for  future  study. 
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